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Patent Applied for 


Now standardised on many 
new aircraft types 


Ministry of Supply and Air Registration Board approval 
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As so often before, British aircraft designers show 
Once again that aircraft can be built to give still 
greater safety, comfort and economy. From 
Vickers-Armstrongs comes the world’s first pro- 
peller turbine airliner—the Vickers Viscount 700. 
For safety, the Viscount 700 has four Rolls-Royce 
‘Dart’ engines using kerosene fuel; crash-proof 
tanks: twin wheels to each undercarriage unit: 
thermal de-icing; and large doubie-slotted flaps 
making possible the use of small undeveloped 
airfields. For comfort, there is 4 fully pressurized, 
air-conditioned cabin; the noise level has been 
reduced to a figure never attained before and 
vibration almost eliminated. For economy, the 
cruising speed is 325 m.p.h.—with a full load of 
53 passengers—giving a fuel consumption of 1.45 
air miles per gallon. 


VICKERS-ARMSTRONG LIMITED 


AIRCRAFT SECTION: 


Head-on view showing the twin wheels and the high-set tail 
plane which leaves the trim setting unaffected by changes of 


engine power or when lowering the flaps and undercarriage. 


A view showing part of the cabin where passengers can relax 


in quiet comfort never before experienced in air travel. 


VICKERS HOUSE BROADWAY * LONDON swf 
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The right 
approach 
to every 
aeronautical 
problem 


The factual reporting, informed commentaries and ONE SHILLING 
unrivalled technical drawings which have distinguished AN bes nn 
‘The Aeroplane’ for the past 38 years still make £3.1.0 

a careful study of its pages the right approach to POST FREE 


every aeronautical problem. 


TEMPLE PRESS LIMITED 


BOWLING GREEN LANE, LONDON, E.C.1. TERMINUS 3636 


Associated Aeronautical Publication: 
““THE AEROPLANE’ DIRECTORY” incorporating “WHO'S WHO IN BRITISH AVIATION.” 
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Weston Thermometer Test Sets 


Weston Thermometer Test sets are designed to provide an accurate, 
simple and speedy means of checking electrical aircraft thermometers. 
Two models are available—the Model S76 for resistance type 
thermometers and the Model S77 for thermocouple types. Both 
models are portable and self-contained and are direct reading. They 
will test all English type thermometers which conform with A.M. 
specifications and in addition various American types. These Weston 
Test Sets are an invaluable asset to all airports, both for routine 
inspection in the field and for overhaul and service depots. Please 


write for details. 


SANGAMO WESTON LIMITED 


GREAT CAMBRIDGE ROAD, ENFIELD, MIDDLESEX 
Telephone: Enfield 3434 & 1242 
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*FLEXELITE’ 
The new LIGHT WEIGHT 


FLEXIBLE FUEL TANK 
of SYNTHETIC RUBBER 


and NYLON FABRIC. 


MAREX* 


LIGHT ALLOY 
HEAT EXCHANGERS 
COOLERS, INTERCOOLERS 
RADIATORS, ete., ete. 


Fullest co-operation and advice from our Development Departments 


MARSTON 


(A bsidi 


EXCELSIOR LTD 


y of Imperial Chemical Industries Ltd.) 


WOLVERHAMPTON (TEL. FORDHOUSES 2181) AND LEEDS 


MAR.51a 


43 PLACES & 3 BADEN POWELL 


MEMORIAL PRIZES GAINED IN R.AE.S. EXAMINATIONS 


\ 
nd careet g at 


by home-study students of The T.I.G.B. 
furnish satisfactory proof that The T.I.G.B. 
courses are an authoritative means of equip- 
ping men in the aeronautics industry with 
technological knowledge. 


You can start any time on a T.I.G.B. pro- 
fessional course for the A.F.R.Ae.S., A.M.I. 
Mech.E., A.M.I.E.E., or other institutional 
examination in which you are interested; but 
the best time is NOW, because Government 
and industry both require recognised techno- 
logical attainment to be in the possession of 
those who aspire to higher posts. Write to- 
day to The T.I.G.B. for ‘‘ The Engineer’s 
Guide to Success,’’ giving particulars of your 
technological and qualification requirements. 


The Professional Engineering 
and Aeronautics Tutors 


THE TECHNOLOGICAL INSTITUTE 
OF GREAT BRITAIN 


39 Temple Bar House, London, EC4 
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The industrial, transport and 
automotive uses of ESSO have long been 
known in many countries. And ESSO 
Aviation Products are to be found along 
the airways of the world. The operators 

of large or small aircraft, whether commercial 
carriers or private owners, are now looking 


to the famous ESSO oval for high quality facilities please write to Anglo-American Oil Co. Ltd., 
aviation petroleum products. 


Aviation Dept., Artillery House, London, S.W.1. 


AVIATION PRODUCTS 


FUELS e LUBRICANTS e DE-ICING FLUIDS 
RUST PREVENTIVES «+ SPECIAL PRODUCTS 


For contract terms and foreign travel 


ANGLO-AMERICAN us LIMITED 
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Inereased and stabilised yields..,, 
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The helicopter’s role in food production 


Fundamental advantages of the helicopter for crop-spraying, pest-control, and fertiliser-laying duties are (a) its 


great maneeuvrability ; (b) its ability to fly low at slow speeds; (c) simplicity of taking-off and landing within 
confined spaces; (d) its ability to cover large areas in a short time. The “Bristol” Type 171 helicopter has 
been constructed and designed to give a fully versatile performance with an exceptionally high safety factor. 


Agricultural development is but one aspect of many industrial and social services offered by this type of aircraft. 


THE use AEROPLANE COMPANY LIMITED ENGLAND 


Daily Express, May 4, '49 Co! Internatio” 
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; the Home produced 


Elektron” magnesium is abundant in the 
sea-water from which it is obtained. Two-thirds the’ 
weight of aluminium, it has been successfully proved in 


its 

in aircraft for over twenty-five years. New “Elektron” 
as magnesium-zinc-zirconium and magnesium-zirconium-rare-earth 

‘ alloys are readily cast, extruded, forged or welded. 

5. Ample supplies are available! 

HUGHES & COMPANY LIM! 
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Australia’s International Airline 
provides complete modern facilities 
for Air Travel, Air Mail and Air Cargo 


SYDNEY—LONDON 

via Singapore, India, Egypt 

by ‘Kangaroo’ Service. 

By Constellation (with B.O.A.C.) 

@ SYDNEY—NEW GUINEA 
Bird of Paradise Service by DC.3 Airliner. 
Sydney — Northern Queensland Airports — New 
Guinea — Rabaul. 


@ ISLAND SERVICES 
Sydney — Norfolk Island. Sydney — 
Noumea — Suva. Sydney —Lord Howe 


Island. 

SYDNEY—AUCKLAND 
Trans-Tasman Service 
(with T.E.A.L.) 
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QANTAS EMPIRE AIRWAYS 


in parallel with BRITISH OVERSEAS AIRWAYS CORPORATION 
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SPERRY 


GYROSCOPIC INSTRUMENTS FOR NAVIGATION 


AND PRECISE CONTROL OF AIRCRAFT 


<< ELECTRIC GYRO HORIZON Type H.L.3 


A sectioned drawing of the movement 


With the advance in design of British Aircraft has come the require- 
ment for new instrumentation to meet the exacting needs of high-speed 
flight at extremes of altitude and temperature. 


To meet these requirements Sperry have developed new precision-built 
instruments, which will give the necessary accuracy and reliability. 
These are now in production. They have the full approval of the Air 
Registration Board. The Electric Gyro Horizon provides an accurate 
indication of aircraft attitude, is topple-free, and has been specified on 
all new British Transport aircraft. 


GYROSYN COMPASS, C.L.1and €.L.2 GYROPILOTS, A.3, A.L.1, A.12 
The latest design of Gyro-Magnetic Compass . DIRECTIONAL GYRO 

for aircraft. Free from Northerly turning " 

errors, acceleration errors, and Gyro wander. e AIR DRIVEN GYRO HORIZON 
Specified by the British Corporations. 


GRUst WEST ROAD, BRENTFORD, MIDDX. TELEPHONE : EALING 6771 (10 LINES) ee 
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THE SECOND LOUIS BLERIOT LECTURE 


FRENCH PRACTICAL AERODYNAMIC 
METHODS 
by 
JEAN BROCARD and FRANCOIS HUSSENOT 


THE SECOND Louis Bleriot Lecture and the 766th Lecture to be read before 
the Royal Aeronautical Society, was given at the Institution of Civil Engineers, 
Great George Street, London, §.W.1 on Thursday, 10th February 1949. Dr. H. 
Roxbee Cox, F.R.Ae.S., F.LAe.S., the President, presided. 

The President: It was on 25th July 1909 that Louis Bleriot flew from France 
and landed on the slopes near Dover Castle, winning the Daily Mail Prize. He not 
only made one of the great flights of history, but he made history itself. It was not 
surprising, therefore, that when the Association Francaise des Ingénieurs et Tech- 
niciens de |’Aéronautique decided to found an annual lecture, they should name 
it after him. They also had the very kind thought of inviting the Royal Aeronautical 
Society to be associated with that annual lecture and it had been arranged that the 
lecture should be given first in France by an Englishman and then in England by 
a Frenchman and so on. Thus the name of the man who first joined France 
and England by air became the name of the lecture which provided a fraternal link 
between their two societies, a link which they hoped would last as long as their 
two societies. 

It had been his good fortune to be President of the Society when the First Louis 
Bleriot Lecture was given in France in 1948, and it was a privilege to be present on 
that inaugural occasion. He counted it a great privilege also to be in the Chair on 
the present occasion for the Second Louis Bleriot Lecture, and the first one in this 
country. 

A number of their French friends had come over for the meeting and he was 
sure all the members of the Society would wish him to welcome them. 

They also wished to welcome their two lecturers, M. Jean Brocard and M. 
Francois Hussenot, and M. J. Jarry, President of the A.F.I.T.A., who would introduce 
the lecturers. 

M. Brocard would read the first part of the Lecture and Professor A. A. Hall, 
at short notice, had kindly agreed to give a résumé of the rest of the Lecture. 

Monsieur Jarry: On the occasion of the Second Louis Bleriot Lecture it was an 
honour as President of the A.F.I.T.A. to have the privilege of introducing the authors. 

Monsieur Brocard graduated from the French Higher National College of Aero- 
nautics in 1929 and that same year served as a Lieutenant in the French Air Force. 
— to 1940 he served as a Reserve Officer, flying many types of French 
aircraft. 

He joined the Breguet Aircraft Company in 1931 and since then had been in 
charge of the Aerodynamics Department and the Wind Tunnel testing of all Breguet 
types. He had done much research on seaplane floats, powered models, propellers 
and engines. He was at present Chief Engineer of the Louis Breguet Company and 
Director of the Breguet Wind Tunnel. 
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He was a member of Council of the A.F.I.T.A. and President of the Aero- 


dynamics Committee of the A.F.I.T.A. 


Monsieur Hussenot attended the French Polytechnique from 1930-32 and then 
served for a year with the French Air Force testing new aircraft. 

From 1933-35 he attended the French National College of Aeronautics and then 
joined the Flight Testing Centre as an engineer, taking part in the testing of proto- 
types of both landplanes and seaplanes, and in operations on aircraft carriers, 

At present he was Chief of the Methods Department at the Flight Testing 
Centre at Bretigny. M. Hussenot was the inventor of several multiple registering 


Flight Recording Instruments. 


He was a member of the Council of the A.F.I.T.A. and President of the Aero- 
dynamics of Flight Committee of the A.F.I.T.A. 


PPRANCE. paralysed during many years of 

a crue! occupation, has been unable to 
maintain a sustained technical effort, com- 
parable with that of England. 

Nevertheless, the work of French engineers 
was never discontinued. In difficult and 
sometimes dramatic circumstances, more 
often than not with only restricted facilities, 
their work of research continued, stubbornly 
and in silence. 

A survey of this research work has been 
made by the Association Frangaise des 
Ingénieurs et Techniciens (A.F.I.T.A.), and 
a number of aerodynamic test results have 
been brought to light. It was thought that 
these tests might well be of interest to 
aeronautical technicians in view of their 
methodical nature, both in wind tunnel and 
flight test stages. 

This paper is the result of teamwork. The 
team is a large one, consisting of specialists, 
either from Government organisations or 
from national and private aeronautical con- 
struction establishments. 

We are indebted to the Directors of the 
Government departments concerned and the 
Presidents of the various Companies who 
have made it possible to undertake this work, 
and also to the engineers for their contri- 
bution of the best results obtained. 

The paper will not describe experimental 
methods well known to all; our purpose is to 
discuss the ways and means of obtaining 
from experimenial work, results which will 
speak for themselves and which it will be 
possible to put to practical use; results which 
will afford the engineer the best possible 
chance of anticipating correctly the handling 
characteristics and performances of the 
prototype for which he is responsible. 

By attempting to obtain results with only 
restricted facilities at their disposal, the 
French engineers followed to the best of their 
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ability the example shown almost forty years 
ago by their great elder Louis Bleriot. 

This paper, for the main part, is devoted 
to low speed aerodynamic testing. This may 
appear paradoxical at a time when yet higher 
Mach numbers are being called for daily in 
the skies, by the Vampires, Meteors, 
Attackers, Swallows, Shooting Stars and 
many others. 

Nevertheless, although the aircraft which 
carries human lives may have been designed 
to attain speeds up to and even beyond that 
of sound, the start is from the ground and at 
zero speed. First there is the take-off run, 
then the take-off and the gathering of speed 
and altitude. The cycle is then reversed and 
the pilot must land and arrest his machine 
without having, for one single moment, lost 
control of the aircraft. 

It follows that the engineer must provide 
excellent characteristics of controllability, 
manceuvrability and lift, even at the “low” 
speeds for take-off and landing. 

The considerable effort which is being 
made all over the world in the field of high 
speed research must not neglect low speed 
conditions. 


It is our belief that the low-speed tunnel 
will always retain its place side-by-side with 
the more powerful transonic and supersonic 
tunnel. 

We know, through classical mechanics, 
that the motion of a body is governed by two 
sets of equations: the equations of forces and 
the equations of the moments of these forces. 
For the aeronautical technician, the “ per- 
formances ” of the prototype are governed by 
the equations of forces, whereas its handling 
characteristics are governed by the whole of 
the equations. The art of the engineer con- 
sists in anticipating both as accurately as 
possible. 
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Are wind tunnel tests reliable for anti- 
cipations of this nature? ; 

And to what extent shall the flight test 
results bear out wind tunnel anticipations? 
That is the subject which we propose to 
discuss. We do not claim to be supplying 
conclusive answers to these “eternal” 
questions. We merely intend to review 
certain methods which have been used in 
France for several years and the positive 
results obtained. 


|, AERODYNAMIC TESTS 


The knowledge of the already firmly estab- 
shed laws of aerodynamics, augmented by 
the daily contribution of an abundant 
technical literature, together with the experi- 
ence gained through previous achievements, 
will of course allow for “ roughing out” an 
aeroplane project. 

There are a considerable number of points 
which are not covered by calculation how- 
ever. 

Fach new aeroplane offers its own 
problems and these require numerous aero- 
dynamic tests, involving longitudinal and 
lateral stability, control surface effectiveness 
and forces, drag, high lift and side slipping. 


l.l. THE SLIPSTREAM 


One of the leading factors in the develop- 
ment of an aircraft is the slipstream. 

The slipstream has a great influence upon 
the flow pattern and is the cause of local over 
speeding and downwash. It is accompanied 
by a tip vortex sheet which may alter 
considerably the characteristics of the control 
surfaces. The intensity and shape of these 
vortices vary considerably with engine r.p.m., 
forward speed, propeller pitch and direction 
of rotation of the propellers. 

The Froude theory enables the approxi- 
mation of the mean value of the propeller 
slip and of the total downwash rearward of 
the wing, but cannot be used for estimating 
the special affect of a given axial or 
tangential induced velocity, or for assessing 
the effect of the slipstream upon the engine 
nacelles, the wings, the fuselage and the 
control surfaces. Nor is this theory of any 
use for estimating the effect of a vortex 
coming into contact with the tailplane of the 
aircraft. 

Figure | shows how important this vortex 
sheet action can be. 

It will be seen that a vortex appears to be 
crossing the tailplane of the aircraft. This 


phenomenon, which has often been observed 
in wind tunnels, is accounted for by the 
violent rotational motion which is imparted 
to the air. After the air has left the surface 
situated in its path, the rotational motion is 
still sufficiently intense to form another 
vortex further downstream. 

Since the engineer cannot “kill” the 
vortices, he must compromise with them. It 
is essential therefore that they be reproduced 
accurately in the wind tunnel, so that the 
engineer may become thoroughly familiar 
with them. Hence the necessity to equip the 
aerodynamic models with engine/propeller 
combinations capable of simulating the true 
streamline flow of actual propellers. These 
models, when properly used, are an invalu- 
able source of information and are known as 
powered models.” 


1.2. DEVELOPMENT OF THE POWERED 
MODEL IN FRANCE 

As early as 1934, the work of Mr. Jacques 
Valensi on the visualisation of the slipstream 
(published in No. 73 of the French Air 
Ministry’s Scientific and Technical Publi- 
cations) had brought home to the constructors 
the importance of taking into account the 
slipstream structure in wind tunnel tests. It 
was the difficulty of producing engines small 
enough to be fitted in scaled-down models 
and at the same time, powerful enough to 
meet mechanical similarity requirements, 
which delayed the development of powered 
models until the beginning of 1936. 

Then, several constructors, having at last 
found the answer to the problem of electric 
motors with the help of qualified electrical 
engineers, among whom we would mention 
the Seguin brothers and E. Grambin, were 
able, for the first time, to make a number of 
powered-model tests, thanks to the support 
of Chief Engineer Paul Dupont, Head of the 
Aerodynamics Division of the French Air 
Ministry. 

These constructors were: first, the Société 
Morane, which made visualisations of the 
air flow, then Les Ateliers d’Aviation Louis 
Bréguet, which were the first in France to 
make a methodical “ weathercock investi- 
gation ” of the influence of the slipstream on 
the stability of an aeroplane; they were 
followed closely by the Société Dewoitine 
(French Air Ministry Technical and Scientific 
Publications, No. 128, 1938). 

The tests were made in the Air Ministry 
tunnel at Issy-les-Moulineaux, of which M. 
Micheau was then the Director. 
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The results obtained clarified many points 
hitherto unexplained, and the constructors 
were soon able to bring to a conclusion the 
development of the prototypes which they 
had on test at the airfield. 

Powered-model testing soon became wide- 
spread among the constructors and was the 
cause of many new developments of which 
we shall give a few examples. 


1.3. FRENCH POWERED-MODEL TECHNIQUE 
—COMPARISON WITH FLIGHT TESTING 


The finished aircraft is the result of a 
careful compromise between handling 
characteristics and performance. The air- 
craft must have stability in the broad sense 
of the term, and at the same time be 
maneeuvrable; that is, there must be sufficient 
response to the control surfaces; another 
requirement is controllability, whereby the 
controls are easy for the pilot to handle. 
Drag must be brought down to a minimum. 
Those are the problems. 

The following is a sequence of tests for 
finding the answers:— 

First, the aeroplane project is “ roughed 
out” by using a non-powered model, the 
scale of which is substantially smaller than 
that of the powered model; then, the latter 
is made. It is an accurately scaled-down 
version of the prototype; control surfaces, 
high lift devices and internal air circuits are 
included. Electric motors fitted with adjust- 
able pitch airscrews provide mechanical con- 
- similar to those of the actual power 
plant. 

Equipped in this manner, the powered- 
model is ready for giving many indications 
on longitudinal stability, deflection and forces 
on the control surfaces, effect of thrust and 
slipstream on the trim of the machine and on 
the performance. 

These tests are delicate operations and a 
number of precautions must be taken:— 

First, special care must be given to the 
propellers. 


1.3.1. Powered-model propellers 


For a complete defintion of propeller oper- 
ation, it is necessary to know the forward 
Velocity V, the revolutions per second n, the 
power absorbed W, the diameter D, and the 
specific mass of the air p. 

The expression of power absorbed: 


W=ypn*D*, 
together with the expression of the thrust: 


T=7pn’D*, 
will give the efficiency; 
y nD 
where J=V /nD. 


For each value of pitch angle «h*, there 
are the corresponding functions: x=f (J), 

Curves of x, 7, 7 in terms of J for given 
values of pitch make up the complete 
characteristic curves of the propeller. 

We first ask to what extent the model 
propeller characteristic curves, which define 
dimensionless coefficients, are influenced by 
scale effect, since diameters ranging between 
6.2 in. and 15.7 in. are frequently used to 
represent propeller diameters of 13 ft. and 
above. 

A substantial loss may be feared in the 
efficiency of the models or, the efficiency 
remaining unaffected, a sharp falling off of 
power. This would result in the relative 
thrust values being underestimated and in a 
decrease of slipstream effect. In all these 
cases, the tests would give incorrect results. 

It might be hoped to overcome these diffi- 
culties by the use of artifices, such as altering 
the pitch or speed of an unsatisfactory 
propeller, but it is our opinion that these may 
prove dangerous methods, since vortex sheets 
are created which are no longer similar to the 
vortices of the actual size propeller and, if 
the thrust coefficient - is brought up to the 
actual 7 obtained in flight, the effects of slip- 
stream and wake thus created will no longer 
have anything in common with true con- 
ditions, and the conclusions of the test may 
well be erroneous. 

By taking certain precautions in the con- 
struction of the models and during testing, 
the troubles mentioned may be prevented. 

The scaled-down propeller blades should 
be of metal. Profile dimensions should be 
adhered to within a tolerance of 0005 in. and 
carefully inspected. 

Engine torque should be measured by 
means of a spring torque-meter (first made 
by Mr. Silber). There are various kinds of 
torque-meters, and these may be installed 
inside the nacelles of the model and directly 
connected to the motor stator. 


* The pitch angle of a propeller normally refers 
to the angle between the zero lift axis of a given 
section and the plane cf rotation. The section 
normally chosen is at 70 per cent. of the blade 
radius, which is taken from the axis of rotation. 
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GRAMBIN ELECTRICAL 
MOTOR EQUIPPED WITH 
| TORQUE STAND 


GRAMBIN ELECTRICAL MOTOR 
WITH PROPELLER SHAFT 


EXTERNAL TORQUE - MEASURING STAND 


MONNIN PITCH- ADJUSTING APPARATUS 


Fig. 2. 


Equipment for wind tunnel measurements of propeller characteristics. 


If it is desired to prevent a rise in temper- 
ature which might affect the accuracy of 
torque readings, the spring is placed outside 
the tunnel experimental chamber and con- 
nected to the engine by means of an 
extremely thin wire which offers practically 
no aerodynamic resistance. Some of this 
equipment is shown in Fig. 2. 

By way of example, we will review a 
methodical investigation made by the 
Société Nationale du Centre. This is a 
comparison of characteristic curves of a 
number of propellers of various diameters, 
but geometrically similar. 

The diameters of the propellers tested 
were:— 
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0.167 m., 0.180 m., 0.200 m., 0.287 m., 
0.300 m., 0.324 m., 0.400 m.* 

The testing equipment included:— 

A torque-meter with spiral spring 
attached to the propeller. Remote read- 
ing of the deflections was done by means 
of binoculars. ay 

A special “ thrust balance ” with sensitivity 
adjustment. 

A “Stroboret” for measuring rotational 
speeds. 

The pitch settings were adjusted by means 
of a “ Monnin” pitch measuring instru- 
ment. 


* For conversion one metre = 3.281 ft. 


one Km.p.h.=0.6214 m.p.h. 
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Curves and photographs contributed by the S.N.C.A.C. resulting from research work by 


Société Nationale de Constructions Aéronautiques du Centre. 


Fig. 3. 


Monsieur Monnin. 
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[SAINTFCYR WIND-TUNNEL MODEL TESTS ( 27.5” DIAMETER) 
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Comparison of take-off propeller efficiencies, 


Curves contributed by The Société Nationale de Constructions Aéronautique du Sud-Ouest 


Photographs of these propellers and the 
results obtained are shown in Fig. 3. 

The following conclusions were drawn from 
these tests:— 


1. Although the loads given by propellers of 
this size are small, the experimental points 
made a regular curve, due to the precision 
and accuracy of the instruments used. 

2. The efficiencies obtained, which were of 
the order of 82 per cent. agree with the 
efficiencies which may be expected from 
the actual propeller. 

3. The results obtained with similar pro- 
pellers of different diameters are in good 
agreement. 

4. The results obtained with two identical 
propellers are in perfect agreement. This 
proves the possibility of thoroughly 
identical manufacture and adjustment. 

In this same field, the Société Nationale 
du Sud-Ouest has made methodical com- 
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parisons of the efficiencies given by propellers 
of different diameters. 

Efficiencies, measured at l'Institut Aéro- 
technique de Saint-Cyr, of a 0.70 m. 
diameter propeller were compared with 
efficiencies, measured at the Hispano-Suiza 
wind tunnel, of an actual size propeller. 

The variation of efficiency with J for both 
these propellers is shown in Fig. 4. 

These curves are scaled in values of pitch 
angle and are very close to one another. 

It appears, therefore, that small diameter 
model propellers can have efficiencies very 
close to those of the actual size propellers 
they represent. ~ 

The agreement between the characteristics 
of the actual size propellers and those of the 
scaled-down models is not unduly surprising. 

In fact, the Reynolds number for the 0.70 
tip radius of a 3.85 m. diameter propeller 
running at 1,300 r.p.m. and having a forward 
velocity of 400 km.p.h. is of the order 
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of 4.10°; the corresponding Reynolds number 
for the 1/18-scale model of the same propeller 
will be 150,000. 

Therefore, the mean profile drag of the 
model propeller blade sections should be 
increased by 40 per cent. in relation to the 
profile drag of the actual size blade. 

The vortex theory, when applied to the 
propeller, shows that the loss of efficiency due 
to the 40 per cent. profile drag increase will 
be less than 2 per cent. 

Such approximation in a test is satis- 
factory. 

There is, however, one restriction, and this 
is for rotational speeds corresponding to static 
rp.m. or to very low values of J; the maxi- 
mum C;, values of the model propeller blade 
sections will in effect be lower than in the 
case of actual size propellers, with a possi- 
bility of greater losses. But it should be 
pointed out that these values of J do not 
generally correspond to the values of J 
usually obtained in flight. 

We shall now give a few examples of tunnel 
tests with powered models. 


1.4. HANDLING CHARACTERISTICS TESTS 


1.4.1. Study of longitudinal static stability 
in steady straight flight—Com- 
parison with flight tests 


14.1.1 Principle of method used 


The conditions for steady flight of an 
aeroplane are governed by the following 
parameters: (a) of the aircraft: 


Configuration of the aeroplane 
Flight incidence i 
Velocity along the flight path V 
Altitude of flight: 

and (b) of the power plant: 
Horse power 
Revolutions per second of propeller 
Pitch of propeller. 

To study longitudinal stability for a given 
flight condition is to assess the value of the 
ratio dM/di for the equilibrium incidence 
corresponding to this flight condition, M 
being the moment about the centre of gravity 
of the aerodynamic forces acting on the 
aeroplane. 

A definite value of the lift coefficient C, 
of the aircraft corresponds to this incidence 
and there is also at a given altitude a corres- 
ponding velocity which defines, the airscrew 
rotational speed being known, a value of J. 

By a suitable measure made on the model, 
the value of pitch angle which will absorb 


the power required for meeting the require- 
ments of mechanical similarity can be 
deduced. 

The angle of deflection 8. of the elevator 
can then be adjusted so as to eliminate the 
aerodynamical torque about the centre of 
gravity. 

The model can be mounted either as an 
aerodynamical “ weathercock,” as shown in 
Fig. 5, or on a special balance. 

Then a few points of the curve represent- 
ing the pitching moment coefficient about 
the angle of equilibrium are plotted in order 
to determine whether the stability is satis- 
factory. 

Furthermore, an examination of the model 
oscillating freely about its centre of gravity 
is most interesting; it occurs frequently that, 
at a certain incidence, the tailplane is 
immmersed in the propeller tip vortex sheet. 
It can then be seen that the model is oscil- 
lating about its equilibrium position, the 
amplitude of the oscillations being governed 
by the incidences at which the tailplane is 
not subjected to the vortex sheet. 


1.4.1.2. Influence of airscrew vortex sheet. 
—The Bréguet 460 aircraft 


Figure 6 shows a series of methodical wind 
tunnel and flight experiments during which 
this phenomenon was brought out very 
clearly. 

The photograph shows the twin-engined 
Bréguet 460 aircraft, and the powered model. 

Longitudinal stability curves are drawn for 
various directions of rotation of the pro- 
pellers, all other things being equal. 

With the propellers rotating inwardly, the 
tailplane is, at the test incidence, almost 
completely immersed in the tip vortex sheet. 

The stability curve which gives the vari- 
ation in aerodynamic moment coefficient 
about the centre of gravity, against incidence, 
shows certain anomalies. Certain regions of 
the curve indicate model oscillations due to 
tip vortices. 

With the propellers rotating outwardly, 
there is quite a noticeable alteration in the 
aspect of the vortex sheets; in this instance, 
the vortices are, to a certain extent, away 
from the tailplane. The regions of oscillation 
are eliminated and the stability gradient 
dC,,,/di is substantially increased. 

These wind tunnel observations have been 
completely confirmed by flight tests. In this 


connection, Fig. 6 shows curves of control 
column force plotted against velocity. 
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Fig. 5. 


Vertically supported model in the Breguet wind tunnel for weathercock instability measures. 
Ateliers d’Aviation Louis Breguet. 


These indicate a reversal due to longi- 
tudinal instability experienced in flight with 
inwardly rotating propellers, whereas stability 
was excellent with outwardly rotating pro- 
pellers, and the curves of control column 
force plotted against velocity were normal in 
this latter instance. 


1.4.1.3. One example of thorough wind 
tunnel investigation._Bréguet 730 
flying-boat. 


Thorough experimentation on models of 
different scales was done in various wind 
tunnels prior to the testing of the Bréguet 730 
long-range flying-boat. 

The main concern was to obtain the best 
possible handling characteristics and at the 
same time the performances required for the 
aircraft. 


The tests were made:— 
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to 


scale model which has a span of 2.01 m., of 
the 1/4-scale model which has a span of 
10.09 m., and of the prototype (span 
= 40.36 m.). 


stability curve obtained in the same con- 
ditions with 1/30, 1/20 and 1/4-scale models. 


. At the Eiffel laboratory with a non- 


. At the large size tunnel of Issy-les- 


. At the large size tunnel of Chalais-Meudon 


. At Le Havre and at Cherbourg, a series 


powered 1/30-scale model. 


Moulineaux with a 1/20-scale powered 
model. 


with a 1/4-scale powered model. 


of flight tests was made with the proto- 
type in order to check tunnel test results. 
It will be seen that the flight tests con- 
firmed the tunnel test results in a most 
definite manner. 

Figure 7 shows photographs of the 1/20- 


The figure also shows a_ longitudinal 


PysH 
T 


= 
| 
Pi: 


Longitudinal stability test. 


Influence of different directions of rotation of the propellers all other things being equal. 
CONVERGENT SMOKE VORTEX 


DIVERGENT SMOKE VORTEX 


Breguet 460. 


100Cm 
| | | | | | 
| 
OSCILLATIONS OF MODEL } , ELIMINATION OF OSCILLATIONS 
| | | 
| | | | 
| | 
= | | | & 
| | 
| 
WINO-TUNNEL TESTS ON MOTORIZED 76, SCALE MODEL 
PITCHING MOMENT COEFFICIENTS PLOT TED 
Kes | 
z4 + x4 | | 3 
IRECTION OF ROTATI = + 
| | | 
WO 160 100200 220 240 260 SPEED IN KmH 140 «160 180 200 220 240 280 SPEED IN KmH 
Fig. 6. 


Control column forces. Flight tests. Ateliers d’Aviation Louis Breguet. 
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Breguet 730 No. 1 moored in Powered model mounted on Powered Model (scale: { /20), 
Havre. the Chalais Meudon Wind- 


Tunnel Balance (scale: 4). 


BREGUET 730 FLYING-BOAT BREGUET 730 FLYING-BOAT 
le Distance between Control Column and Instrument Pannel, Control column position plotted against speed 
150! | 110, ] | 
1 190 + 
| boat 
200 + 
| | —o— Flight test with C.G. ot 36.25% of reference chord (7-1v-%) 
—o— Flight test with C.G. at 36.25% of reference chord (28-1v-8) 
100120. 40. 60. 180 200. 220 ™O 260 ~ - Scale model wind-tunnel test (C6. at 38% of reference chord) 
s Flight tests : Flaps down + 25°. Atlerons neutral 0° - ~~ ~Y4 Scale model wind-tunnel test (C.6 at 38% of reference chord) 


4 scale model Flaps down +15°. Ailerons 
Fright tests: Flaps neutral Ailerons neutral O° 
—o— 14 scale model tests. 0.85. Flaps down Ailerons +10° 


BREGUET 750 FLYING-BOAT 


COMPARATIVE RESEARCH COVERING REACTIONS 


16 


| | } | 
ON CONTROL COLUMN MEASURED IN FLIGHT AND ) —— scale model | 
DURING WIND-TUNNEL TESTS —o— ‘50 scale mode] | 


[Yeo scale model] | 


All models centered at 38% 


| with | Flaps neutral L 

8 
t 
of 

Speedo 
| | | | 
-6 + 

| 


-- YAScALe MODEL TESTS WITH 50 % AUTOMATIC COMPENRATOR 
FLIONT TEST WITH 85 AUTOMATIC COMPENSATOR 

- - SCALE MODEL TESTS WITH 100% AUTOMATIC compens® 
FLIGHT TEST WiTH 50% auTomaTic comPeNsaToR 


LONGITUDINAL STABILITY, 


Pitching moments are considered as 


Fig. 7. 
Ateliers d’Aviation Louis Breguet. 


Most satisfactory agreement between tests for each velocity V, the angle of deflection 
obtained in the various tunnels is shown by 8, of the elevator for given engine conditions 


these tests. (for example at rated r.p.m. and full throttle, 
or at low boost). 
Comparisons with flight tests For each velocity V, there is a corres- 


One of the means of checking the flying ponding flight incidence i. 
qualities of the aircraft consists in measuring The sign of the expression d./di, obtained 
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positive, according to French standards. 
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fom the tests, defines the longitudinal 
ability of the aircraft. 

On Fig. 7 curves of control column move- 
ment—{the control column position relative 
o the instrument panel defines the angle f. 
of elevator deflection)—are plotted against 
velocity V measured on a 1/4-scale model at 
the Chalais-Meudon wind tunnel and in 
fight on the prototype. 

The two curves nearly coincide. 

Later, a methodical investigation was made 
of the effect of camber flap deflection on 
these curves of control column movement, 
and tunnel tests were carried out with the 
faps deflected. Curves were also established 
in flight tests and the comparison of the 
results obtained is shown in Fig. 7. 

It will be seen that there is a satisfactory 
agreement between actual conditions and 
that which had been anticipated. Further, 
at a constant velocity, deflection of the 
camber flaps causes a forward movement of 
the control column in order to maintain the 
equilibrium of the aircraft. Therefore, flap 
deflection tends to pull the aeroplane’s nose 
up, although nose-heaviness of the wing is 
set up. 

This experiment proves that in this parti- 
cular case, it is the aerodynamic deflection 
on the tailplane causing a tail-heaviness 
moment which has most influence on the 
general equilibrium. 


1.4.1.4. Influence of engine conditions on 
stability 


The influence of engine conditions on the 
equilibrium of the aircraft has been studied 
on the 1/20-scale model, the 1/4-scale model 
and the prototype. 

An aeroplane is said to be stable in flight 
against engine conditions at a given speed 
and therefore a given incidence when, by 
suddenly cutting the engines at this incidence, 
a pitching moment is set up, whereby the 
nose of the aeroplane is depressed. 

If the reverse occurs, the aeroplane will be 

unstable against engine conditions. This 
quality is readily measured on a powered 
model. Tests were made on the 1/20-scale 
model of the Bréguet 730 flying-boat and 
these enabled us to determine accurately the 
Optimum horizontal tail unit position. 
_ An interesting point is that the variations 
In tm after cutting the engines were found 
to be the same for both the prototype and its 
1/20-scale model. 


control 


1.4.2. Pre-determination of 
column forces 


Correct pre-determination of the 
air forces on the control column 
against velocity is of paramount 
importance, particularly with 
large aircraft 


With a large size horizontal tail unit as on 
the Bréguet 730 flying-boat (area=over 33 
sq. m.*), an error in the estimation of the 
hinge moment may be the cause of loads on 
the control column greater than the pilot’s 
physical strength. There is no need to stress 
the seriousness of this. 

The elevator of the 1/4-scale flying-boat 
model, for this reason, was equipped with 
automatic tabs. 

eflection 
elevator deflection 


“automaticity ” is adjustable at will. 

Careful determination of “automaticity” 
adjustment was made in the wind tunnel. 

The sequence of operations was as 
follows: 

For each longitudinal stability test point, 
a measure was made of the elevator hinge 
moment corresponding to the deflection £, 
which neutralised the general pitching 
moment. Then the full-scale loads on the 
control column were deduced by calculations 
based on the control rod layout of the actual 
size aircraft. 

As the next step, the forces acting on the 
control column were plotted against velocity, 
which was easily deduced, the incidence and 
polar of the aircraft being known. 

The wind tunnel tests enabled us to draw 
a curve of the forces on the control column 
for tab “automaticities” of 100 per cent. and 
50 per cent. These curves are shown in 
Fig. 7 together with curves drawn from flight 
test results. It will be seen that for 50 per 
cent. tab “automaticity,” the curves of tunnel 
tests and those of flight tests are in close 
agreement. An 85 per cent. “automaticity” 
obtained by interpolation, which gave accept- 
able loads for the pilot, was finally adopted 
for the aircraft (see Fig. 7). 

The whole of the aerodynamic tests 
amounted to some 2,000 hours of wind 
tunnel work. This may appear excessive, 
but it should be pointed out that the flight 
tests were kept to a minimum since these 
entailed practically no extra development 
work. 


1.4.2.1. 


* For conversion one sq. m.=10.76 sq. ft. 
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TESFS MADE THE ETABLISSEMENT DE CHE AERONAUTIQUE (TouLoube) 
—o-— TESTS MADE AT|THE WIND-TRrINEL | 

L 
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| | } \ | 
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Actual tailplane of the Breguet 890 mounted in the Chalais-Meudon tunnel. 


COMPARED ELEVATOR HINGE MOMENTS OF THE BREGUET 890 


AND OF iTS % SCALE MODEL 


Fig. 9. 


Ateliers d’Aviation Louis Breguet. 
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1.4.2.2 


When it is not possible to make a very 
large size powered model for tests at Chalais- 
Meudon, a smaller scale model will be used 
(span=2 m. approximately) which will give 
useful indications on stability and on the 
polar, but will not be sufficiently reliable for 
tailplane hinge moment measurements. 
These latter will be made on a separate tail- 
unit model, much larger than the aeroplane 
model (for example, 5 times larger). 

These measurements will be the basis for 
calculating the forces on the contro] column 
and the rudder bar. The local velocity will 
be measured on the powered model. 

An example of the application of this 
method is given in Fig. 8. The method was 
used for calculating the stick forces on the 
four-engined Bréguet 482 aircraft. Com- 
parisons with flight tests are also shown in 
the figure and these denote a good agreement 
between results as regards order of magni- 
tude. 

The figure shows a comparison of the 
curves of hinge moments of an actual tail 
unit tested in the large tunnel at Chalais- 
Meudon with those of the 1/5-scale model of 

’ the same tail unit (see Fig. 9). 

It will be seen that the various results are 

very close to each other. 


1.4.3. Application of operational calculus 
to the study of disturbed motions 
of aircraft 


Le Centre d’Etudes de Mécanique du 
Vol, founded in 1947 by Le Service Tech- 
nique Aeronautique of the French Air 
Ministry, set themselves the task of applying 
the modern methods of operational calculus 
and of frequency response to various 
problems of the dynamics of aircraft, par- 
ticularly those relative to stability. 

One of the many interesting results was a 
comparison of a stick fixed dynamic stability 
test with the result of the calculation by the 
operational method. 

The calculation of stability and control- 
lability of aircraft involves linear equation 
systems which may or may not contain a 
right-hand side. The solving of these 
equations by means of the classical analysis 
methods would be most tedious. Particularly, 
research of motion for given initial conditions 
would necessitate a large number of 
operations, and the time factor involved 
would not be acceptable for an Industrial 
Research Office. The amount of work is 
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lessened considerably by applying the 
operational method. 

There are various ways of Presenting 
operational calculus which may be looked 
upon as an aspect of symbolic calculus. We 
shall introduce it in the form of Carson's 
transformation which is identical but for one 
factor to that of Laplace. 

Taking the function f(r) of the variable ; 
identically zero for <0, Carson’s trans. 
formation gives the corresponding function 
defined by 


It is said that ¢(p) is the operational 
equivalent of f (7). 

The advantage of this transformation lies 
in the fact that operational equivalents of 
complicated functions can be very simple 
and, above all, that involved algebraic 
calculations may be replaced by simple 
calculations on the operational equivalents. 
In many cases, particularly in the case of 
differential linear equations, there will be 
advantage in calculating on the operational 
equivalents and then returning to the initial 
function from the results obtained in the 
shape of operational equivalents, by means 
of a table of operational equivalents, or 
the Wagner-Bromwich formula : — 


1 4-106 ert (p) 
| “Pap 
In order to solve linear systems, knowledge 
of the following is sufficient : — 
Let @(p) be the operational equivalent of 
f(r). Then we have 


Function Operational equivalent 


f(r) 
(7) (p)- pf (0) 
f’ (7) (p) - p*f (0) - pf (0) 

unit step 1 

—p|(p-2) 

It is readily seen that a system of linear 
differential equations changes into a system 
of linear algebraic equations. The initial 
conditions enter the equations automatically 
and therefore the solution for given initial 
conditions is arrived at directly. 

Another important simplification 1s 
obtained by applying the principle of super- 
imposition, which leads to integrations with 
the actual solutions, whereas it leads 
to a simple multiplication in the case of 
operational equivalents, using the solution 
corresponding to the “ unit step.” 
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| || 
| 
CALCULATED CURVE 
| | | |-+-+--FtleHT TESTS | 
| | 
| 
100 {20 140 160 180 TIME] (sec) 
| | 
| 
Fig. 10. 


Compared flight-recorded and calculated trim and speed oscillations. : 
Experimental sequence: The control column is moved approximately 4 in. and then fixed in 
its new position. 
Leo 451 aeroplane. 
Société Nationale de Constructions Aéronautiques du Sud-Est (S.N.C.A\S.E.). | 


Further advantages, with which electrical 
engineers have been familiar for some time, 
indicate that operational calculus will play a 
leading part in questions relative to stability 
and controllability. We give as an example 
of application, the response of an aeroplane, 
the Leo 45, to a “unit step” elevator 
deflection measured in flight and compared 
with the result obtained through operational 
calculus. 

Figure 10 shows the variations in trim and 
speed calculated and measured against time. 
Agreement is most satisfactory. 


1.4.4. Lateral dynamic stability in steady 
symmetrical flight 


Let us consider an aeroplane travelling with 
a uniform rectilinear motion, at constant 
incidence and trim, and without side slip. If 
a slight outside disturbance sets up a rolling 
or yawing moment or lateral force, an 
unsteady disturbed motion will be added to 
the initial motion. It is said that the 
aeroplane is stable when the disturbance 
amplitudes tend to disappear, that is to say 


that, after a sufficient lapse of time, the 
amplitudes become as small as is desired. 
The quicker the disturbed motion disappears, 
the greater the degree of stability. 

In the study of dynamic stability, the 
normal practice is to check the stability of 
the aircraft without bothering about the 
degree of stability. 

When using automatic pilots, or for guided 
missiles and pilotless aircraft, the investiga- 
tion must be carried further, that is, we must 
find out whether the machine is stable or not, 
but we must also assess the degree of stab‘lity, 
and often the disturbed motion subsequent to 
an outside disturbance, as well. 

The analysis of the equations of lateral 
motion of an aeroplane (equations of motion 
about the axis of roll, of yaw, and along the 
lateral axis in the case of stick fixed stability) 
shows that there are two main conditions of 
stability indicating that the real parts of the 
roots of an equation of the fourth order are 
negative. The remaining conditions for 
obtaining negative real parts of the roots are 
always verified in practical cases. If R and 
E are certain expressions which are functions 
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S.E.2300 light aeroplane. 


INSTABILITY ZONE 
TT Engine shut-off. 


ZONE 


INSTABILITY LIMIT INCIDENCE 


Ditedral Angle 
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Influence of wing dihedral angle. 
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ig 
Model wing dihedral 6" 
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4 
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SE. 2310. SHORTENED TAIL FIN 


Influence of tail fin surface. 


Fig. 


Société Nationale de Constructions Aéronautiques du Sud-Est (S.N.C.A.S.E.). 


of the aircraft characteristics and flying con- 
ditions, we can write:— 


R>0O (R=Routhian discriminant) 

E>0 (E=condition of static stability). 

In studying the nature of the instabilities 
introduced by the negative sign of the 
previous quantities we conclude that when R 
is negative the aeroplane will “rock” with 
increasing amplitude, and this is called oscil- 
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latory instability. When E is negative, the 
aeroplane will start to turn, and the turn will 
become increasingly narrow; this is called 
spiral instability. 
The lateral stability derivatives for a given 
aircraft are determined by the general con- 
figuration of the aircraft, with the following 
exceptions: the rolling moment rate of 
change with side slip which depends mainly 
on wing dihedral and the yawing moment 
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rate of change with side slip, which, besides 
the fuselage, depends on the size of the fin. 
in order to determine the dihedral and fin 
ize which will provide lateral stability for 
the main flying conditions, stability limit 
curves can be drawn for each incidence and 
therefore each velocity : 


R (x, y)=0 and E(x, yj=0 


For a given value of x, there are two values 
of y between which there is stability. By 
varying the speed, two stability limit curves 
are obtained which represent the limiting 
Cn3 
values of Cly 


with x= 
fi 


against speed, between which 


there is stability. This ratio is measured in 
the wind tunnel and a flight test is made to 


ascertain that when the point (v one 


situated between the two afore-mentioned 
limiting curves, the aeroplane is stable. 

This method was applied to the D.520 S. 
aircraft. 

A correct aircraft was expected fron the 


=f(V) curves obtained with the modei, 
and this was confirmed by the flight tests. 

The Société Nationale de Construction 
Acronautique du Sud-Est (S.N.C.A.S.E.) 
have investigated the influence of wing 
dihedral and tail-unit area on lateral stability. 

This study was made on the four-seater 
SE.2300 touring aeroplane (see Fig. 11). 

Ina first series of tests on dihedral effect, 
the aeroplane’s model was at a given 


_ incidence and the tests were made with the 


engines at full throttle and at small throttle 
Openings. 

It was noted that the slipstream increased 
the value mt and at the same time brought 
the spiral instability limit down. 

Further, dihedral effect was very marked 
up t06—7°, after which it was less apparent. 
Measurements in flight were made in a 
different manner and therefore results cannot 
be compared directly. 

The spiral instability velocity limit was 
measured in flight and the corresponding 
dihedral limit was deduced therefrom. The 


incidence limit was plotted on a graph 
against dihedral. 


The conclusions were identical. 
The slipstream decreases spiral stability; 
improvements due to dihedral are substantial 
when passing from 4° to 6° and are less 
marked when passing from 6° to 8°. 
On the other hand, the quantitative assess- 
ment of the spiral instability limit, based on 
the tunnel test, leads to a definitely optimistic 
estimation which is not confirmed by the 
flight tests. 
Finally, two aeroplanes, the SE.2300 and 
SE.2310, which differ only in regard to 
undercarriage and fin area, were compared 
in a second series of tests on fin area effect. 
The fin area, in the case of the SE.2310, is 
obtained by doing away with the upper 
portion of the usual fin. 
The wind tunnel tests have shown that by 
reducing the fin area the value Cn;/Cl; was 
decreased, but at the same time spiral 
instability limits were also decreased. 
In short, the result was a same degree of 
stability in both cases. 
Flight tests confirmed that the speeds at 
which lateral instability occurs are the same 
in both cases. 

It is concluded from this “Sud-Est” 
investigation that the problem of lateral 
stability has not been solved to the same 
extent as that of longitudinal stability. 

But it should be possible to select from 
several possibilities the one which will give 
the best results. 


1.4.5. Measurements of rudder deflection 
required for compensation of the 
yawing moment set up when two 
engines are “cut” 


The wind tunnel tests were carried out as 
follows : — 

The 1/4-scale model was set at an inci- 
dence of 7°5 corresponding to the maximum 
climbing angle. Two propellers on the one 
side were windmilling and the two propellers 
on the other side were driven by their 
respective engines. These latter were run at 
conditions representing the rated r.p.m. and 
rated torque of the prototype’s engines scaled 
down to conditions of mechanical similarity 
for the model. 

A variation of J was then carried out. 

To each value of J at these engine con- 
ditions, corresponded a given value of 
airscrew thrust and slipstream effect and 
therefore a yawing moment. This yawing 


moment was compensated by appropriate 
rudder deflection. 
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The curve of 4 plotted against J, thus 
obtained, is shown in Fig. 12. 

A point corresponding to the tests of the 
1/20-scale model and a point corresponding 
to the flight tests are shown in the figure, and 
it will be noted that the agreement between 
all these tests is most satisfactory. 

We have given only a few examples of the 
laboratory development work which has been 
done in connection with this aeroplane. 
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SCALE POWERED MODEL 


IN THE CHALAIS MEUDON WIND - TUNNEL 


Fig. 


THIS FIGURE iS DRAWN 
TO SHOW THE SCALE 
OF THE MODEL. 


ATELIERS D'AVIATION 
LOUIS BREGUET 


1.4.6. A method for the measurement in 
flight of aileron effectiveness 


L’Arsenal de  JlAéronautique have 
developed an interesting method for the 
measurement in flight of aileron effectiveness. 
We mention it because of its originality. 

The ailerons on which the measurements 
were made were fitted on the wings of the 
“Emouchet” glider. 
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Aileron efficiency measurements. in flight, on 


The glider was fitted with two water- 
ballasts, one on each wing tip. During take- 
off and climb the water-ballasts were kept 
full When the measurements were to be 
done, the pilot drained one water-ballast by 
means of a remote control, thus applying to 
the aircraft a pre-determined rolling moment. 

This moment was neutralised by an appro- 
priate aileron deflection such as to bring the 
bank indicator ball bearing back to the zero 
mark. The induced yaw was compensated 
by means of the rudder bar so as to cause the 
turn indicator pointer to coincide once again 
with zero. The one factor that cannot be 
eliminated is the roll induced by rudder 
deflection, but the influence of this is com- 
paratively small because of the slight 
deflections involved. 

Before proceeding with the actual measure- 
ments, the pilot flew his machine once or 


an “Emouchet” (Sparrowhawk) Sailplane. 


twice in a straight line of stabilised flight and 
recorded the various control positions, which 
were to be taken as the “zeros” from which 
control movement would be counted during 
the test. 

It was checked that speed had no influence 
on these “zeros.” 

While the measurement work was in pro- 
gress, the pilot kept the machine flying 
straight. The Hussenot-Beaudoin recorder 
was connected by the pilot and he noted his 
control column and rudder bar positions. 

From 10 to 15 such measurements could 
be recorded during a descent. 

Before landing, the pilot drained the 
second water-ballast. 

The following curves were obtained from 
these measurements (see Fig. 13):— 

Ci=f («) and C\=f (7) 
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Reynolds effect on lift. 
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the rudder bar deflection. 


progress. 


(1) Wing of classical form. 
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Fig. 14. 
Société Nationale de Constructions Aéronautiques du Sud-Est— 
Ateliers d’'Aviation Louis Breguet. 
TI 
where C, is the rolling moment coefficient. Although there is a certain dispersion, the 2 
f, mean curve C\=f (z) is a straight line. The 
a slope of this curve defines aileron effective: 
p m 
V*Sb ness. 
x represents control column deflection and + Comparison tests with models are now I 7 
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iRFOIL MODEL. WING $.U.c.10 AIRPLANE RIGGED WITH 
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Fig. 15. 


Calculated pressure distribution tests by Société Breguet and Centre d’Essais En Vol. 


1.5. A FEW TESTS FOR THE PREDETERMIN- 
ATION OF PERFORMANCE 
15.1. Measurement of maximum C,. 


The correct prediction of the maximum C, 
Value for the prototype is most important. 
This value has a direct bearing on the total 
weight of the aircraft and therefore on its 
payload and range. 

Methodical tunnel and flight tests with 
models of different scales have been made by 
Various constructors. Fig. 14 shows Cy, 
Increases obtained by the Société Bréguet 
on the Bréguet 730 flying-boat with 1/30, 


1/20 and 1/4-scale models. The test results 
obtained were in good agreement with the 
N.A.C.A. tests. 

Similar measurements were effected by the 
S.N.C.A.S.E. on a swept-back wing. The 
aircraft tested was a light aircraft of the 
flying wing type. With these aircraft, the 
effect of Reynolds number is less apparent 
than it is with conventional aircraft. This 
was predicted by certain German tests. 


1.5.2. Profile drag measurements 


The Société Bréguet is now proceeding 
with a thorough study of laminar flow 
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23 AIRFOU | 
| 


Fig. 16. 

Drag measurements in flight by pitot traverse method on Breguet L.23 laminar aerofoil. 
Top left: Total multi-pressure head. Top right: Multi-pressure recording gauge. 
Sleeve rigging of an L.23 aerofoil on the wings of a S.U.C.10 light aeroplane. 
Compared drags of conventional and laminar aerofoils (wind tunnel tests). 

Tests by Société Breguet and Centre d’Essais En Vol. 


profiles. Further to a thorough theoretical The Cy of this profile was the object of a + 
study, a laminar flow profile was made, of methodical wind tunnel investigation. f 
which the transition point is very far back, It was measured by direct weighing and by | | 

corresponding to high C, values of the order the Pitot traverse method* in the wind tunnel \ 


of 0.5; this is the Bréguet L 23 profile. at l'Institut Aérotechnique de Saint-Cyr. N 
* The Pitot traverse method:— In practice, only the wake is affected by the fr) 
A test method known as the Figg — integration. 4 
method,” consists in measuring the drag of a 
body by determining the loss of momentum of 
the air current influenced by the said body. Po= upstream infinity static pressure. WF | 
A steady flow is assumed. un : 
The quantities measured are: the total pressure &= total saaiianiaiiel “9 wake. 
and the static pressure in a plane downstream p= static pressure in wake. ; a 
of the profile perpendicular to the general current y= transversal distance to an axis parallel with | 
and also the static and total pressures far up- speed. 


stream of the body. The method was established 
by Betz for incompressible flows and the use of 
the results was notably simplified by Melvill 


1= aerofoil chord. 
It should be noted that this formula is ®. 


B. Jones. applicable to profiles where the flow 1s 
is gi i incompressible. ~ 
sie — ee With a few modifications, it can be 
C,=2) / ( extended to compressible flows and bodies 
Neo—po\ of any form. 
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The results were conclusive, as is shown 
by the curves in Figs. 15 and 16. 

But the extent of laminar flow decreases 
with an increase of the Reynolds number or, 
in other words, the transition point moves 
closer to the leading edge. 

The Bréguet L 23 profile has a negative 

dient pressure distribution, whereby it 
should retain its laminar character even at 
high Reynolds numbers. In order to check 
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recording multi-pressure gauge 
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Profile drag flight recording equipment. 


VI = 450 


this properly, large profile models have been 
made which fit on the wings of a test aircraft, 
the SUC 10 “Courlis.” With this aircraft, 
which has a pusher type propeller, the wings 
are not subjected to the disturbance effect of 
the slipstream. 

Figures 15 and 16 contain photographs of 
the wind tunnel model, the profile “ model ” 
tested on the aircraft, and the measuring 
equipment. The measurement of wake 
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pressures is done by means of a SFIM comb 
placed in the vertical plane parallel to the 
aircraft’s plane of symmetry, and containing 
the profile under study. This comb includes 
37 total pressure heads distributed along the 
same vertical plane and nine static pressure 
heads. A photograph of the comb is shown. 

The pressures detected by the comb are 
recorded on the Hussenot-Beaudoin multi- 
recorder by means of a Hussenot-Beaudoin 
multi-pressure recording apparatus whereby 
21 pressures may be measured. 

A photograph of this apparatus is shown 
in Fig. 17. The main parts are three 
diaphragms which measure the difference 
between a common reference pressure and 
the pressures which arrive from the comb and 
are transmitted by a distributor (one distri- 
butor to each diaphragm). The three 
distributors have respectively 6, 7 and 8 
channels and rotate at the same speed. This 
enables the recording and identification of 
6, 7 and 8 pressures. The photograph in Fig. 
17 represents the multi-pressure recorder, the 
distributor mechanism and a record obtained 
during a preliminary test; the marks corres- 
ponding to each diaphragm can be seen. 

Measurements in flight are now being 
made. 


1.5.2.1. Measurements of drag at super- 
sonic speeds 


At the same time as these profile tests were 
being done, the Société Bréguet carried out 
measurements of drag with supersonic flow 
on a body of revolution, by the Pitot traverse 
method. Measurements were effected at 
Mach numbers ranging from 0.87 to 1.4 and 
the wave drag found for M=1.4 is in good 
agreement with that given by Karmdan’s 
theory (see Fig. 18). 


1.5.3. More comparisons of polars 
measured in flight with rolars 
obtained in the wind tunnel 


These tests involve three aircraft of sub- 
stantially different tonnage:— 
The Bréguet 500—total weight = 14 tons 
The Bréguet 730 =30 tons 
The Latécoére 631 =66 tons 
Figure 19 shows photographs of these 
aircraft and also the polars obtained for each 
machine in the wind tunnel and in flight. 
The polars were obtained as follows:— 
The “ wind tunnel polars” were measured 
with propellerless models. The “flight polars” 
were calculated from the speeds measured 
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in flight, the known engine power, and 
propeller efficiencies measured in the wing 
tunnel. It will be seen that there is satis. 
factory agreement between the various polars 
of the models. 

It will be noted that the drag with the 
Bréguet 500 aircraft is slightly higher than it 
is with the corresponding model. In the case 
of the Bréguet 730 flying-boat, drag is slightly 
lower than with the model. With the 
Latécotre 631 flying-boat, drag is notably 
lower than it is with the model. This 
evolution is explained by the fact that 
Reynolds number effect grows as aircraft size 
is increased. 

Scale effect does not quite compensate for 
additional drag due to certain parts which are 
not represented in the wind tunnel in the case 
of the Bréguet 500 (such as rivet heads, 
joints, fittings and so on). With the Bréguet 
730, this additional drag is substantially 
compensated by scale effect. In the case of 
the Latécoére 631, the Reynolds number is 
clearly the most important factor. 


1.5.4. Helicopter tests 


Thorough wind tunnel experiments were 
made prior to the design and production of 
the Bréguet “ Gyroplane G.11.E.” 

Figure 20 shows photographs of this 
machine and of its 1/4.6-scale model. 


Fig. 18. 
Drag measurements on a body of revolution by the 
pitot traverse method ry subsonic and supersonic 
ow. 
Tests made at le Centre d’Essais de Mécanique 
des Fluides. | 
Ateliers d’Aviation Louis Breguet. 
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G.11 E. Breguet Gyroplane, bench tests. ; 
1/4.6 scale powered model mounted in the Breguet wind tunnel. 
Ateliers d’Aviation Louis Breguet. 
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A model of the “ Gyroplane G.11.E” was 
equipped with a 4 h.p. electric motor. The 
rotor control system was an exact replica of 
that of the actual machine. Particularly, it 
was possible to carry out cyclic blade pitch 
variation and to reproduce the actual flying 
conditions in the wind tunnel. Mechanical 
similarity requirements were met substantially 
by the blades. In the course of numerous 
wind tunnel tests, M. Louis Bréguet was able 
to check the correctness of the theoretical 
research work on rotating wings he has been 
doing for a great number of years, and to 
make interesting comparisons with the 1936 
“Gyroplane Laboratoire Bréguet.” 

During the flight tests with the Bréguet 
“Gyroplane G.11.E,” which were begun 
while this paper was being prepared, 
methodical comparisons with laboratory tests 
will be made. 

These few concrete examples of industrial 
aerodynamic test methods show that the 
engineer is now provided with the means for 
solving beforehand a large number of the 
problems which accompany the initial design 
of an aircraft. 

But it is most important to be familiar 
with the main aspects of the technique of 
these tests, if full advantage is to be taken 
of the possibilities they offer. 


Most important of all, aerodynamic tests 
call for thorough care and extreme precision, 
since the quantities to be measured, forces, 
moments, incidences, are on the whole very 
small. The slightest neglect when carrying 
out measuring operations, the slightest error 
in the design of the equipment, may be the 
cause of secondary phenomena, which may 
well be important enough to conceal or even 
teverse the phenomenon under study. 

Furthermore, industrial aerodynamic tests 
are characterised by their very urgency. The 
test programme is closely linked with that of 
the Research Offices where results are 
awaited in order to decide upon constructive 
solutions. 

The results must arrive “ on time.” 

Many are the instances which could be 
recalled of defects observed on aircraft during 
flight testing, which it is too late to cure by 
means of a logical solution, because it would 
take too long or be too costly to embody a 
major modification to the structure; it is a 
tail unit which cannot be made larger with- 
out a risk of endangering the good behaviour 
of the fuselage. or without taking back the 
CG. position in a prohibitive manner—it is 


an engine nacelle, which causes harmful 
aerodynamic interference, but the shape or 
position of which it is too late to alter. 

There is then no alternative but to adopt 
half-measures, solutions which multiply the 
flight tests with their attendant delays, risks 
and also the considerable expenses they 
involve. 

The wind tunnel must not play the part of 
a “ Defect witnessing office” for prototypes; 
it should, on the contrary, serve as a 
“ Prediction Office ” for measurements to be 
made in order to prevent these defects. 

Aerodynamic test results should be care- 
fully filed. A clear and accurate report 
should be made for each test and should give 
without any exceptions, all details which are 
likely to bring out its purpose and usefulness. 
The report should be made out immediately 
following the test, while impressions are still 
tresh in the experimenter’s mind. 

A well-written report does not merely aim 
at supplying the special item of information 
for which the test was made. It is also an 
invaluable source of information for other 
cases of a similar nature. The information 
will be an efficient contribution to future 
developments and possibly will be used for 
establishing laws and regulations. 

It will prevent the undertaking of tests 
uselessly and also be an important factor for 
reducing research costs, which should be a 
permanent concern to the engineer. 

Lastly, it should be remembered that the 
ultimate aim of any research work and of all 
aerodynamic tests is to obtain, with a mini- 
mum of development work, an aircraft with 
performances and handling characteristics 
true to the predictions which will have been 
made. 

It is essential therefore to have perfect 
knowledge of what is being measured in 
flight and of the qualities required by the 
operators. 

It follows that the engineer and test pilot 
must maintain a close and constant co-oper- 
ation, be thoroughly familiar with the 
respective possibilities of wind tunnel and 
airfield, and analyse objectively the results 
obtained. 

The vigilant observation of these main 
rules which this lecture is endeavouring to 
symbolise and which will enable all those 
interested to obtain the utmost from this 
delicate experimentation, could be termed 
“The Art of Aerodynamic Tests,” if Air 
Commodore Banks’ expression were to be 
applied to this field. 
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Since 1941, and despite the obstacles 
resulting from the Occupation, a new tech- 
nique of flight testing has been developed in 
France. This technique was first applied 
in North Africa and heavy flying-boats 
(Latécoére 631, SE 200 and CAMS 161) 
which have benefited by this equipment have, 
at the same time, been used for its develop- 
ment. 

Co-operation was the basis of this develop- 
ment work. In this connection, we would 
like to mention the “ Station Service ” which 
has been set up for the purpose of centralising 
the ordering and development of new type 
measuring instruments. The “ Station Ser- 
vice” is also responsible for supplying the 
instruments at short notice to the manu- 
facturers who require them for prototype 
tests. 

“Station Service ” has a special factory in 
the suburbs of Paris, where equipment in 
short supply on the French market is designed 
and made. 

The equipment which is finally adopted 
results from our fairly long experience in the 
field of flight measurements, in the course 
of which we were compelled to become 
thoroughly familiar with even the humblest 
aspects of flight testing. 

About 1939, it was realised that, more 
often than not, a flight test outside the normal 
routine of climbs and full throttle level flights 
was deemed a special problem necessitating 
the design of new equipment. We felt that 
this policy was in opposition to modern 
practice, whereby there is no hesitation to 
allot large sums for fooling in order to mini- 
mise production costs. It is indeed true that, 
the fewer the resources available, the greater 
the care must be taken on the tooling to 
avoid failure. 

The “tooling” for flight measurements 
should have the following characteristics:— 


(a) Simple equipment so as to reduce weight 
and cost and facilitate servicing and 
operation. It should be possible to 
start up the whole of the equipment 
either automatically or by means of a 
simple push-button. 


(b) Recording should be on photographic 
paper by projection on to a slot; the 
paper tape used should be wide enough, 
so that photographic enlargement should 
not be necessary before checking 
results. If an additional tape is required, 
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synchronisation must be both automatic 
and accurate. 

(c) True indications should be the main con- 
cern; frequent and accurate calibrations 
as close as possible to flying conditions 
should be effected. 

(d) Standardisation of equipment, spare 
parts, and so forth. 


As an example of the above, we would 
mention that the overall dimensions of the 
recording apparatus have not been altered in 
the past eight years and that the prototype 
is still in service, all production modifications 
having been embodied. We feel that we have 
achieved our objective, which was not merely 
to produce equipment for an aircraft, but a 
real tool to be developed with great care 
before using it on the prototypes of to- 
morrow. 

There was the danger that this “tool” 
would limit the scope of flight tests and 
prevent the development of new measuring 
methods. We feel that these difficulties have 
been avoided, thanks to the “rear framework” 
which will fit on the back of any recorder. 
Instruments of all kinds, such as pendulums, 
gyroscopes, galvanometers, etc., for which 
there is no room on the normal recorder, can 
be installed on this “ rear framework.” 

We will now give an example of the 
utilisation of the equipment which will show 
the possibilities of a continuous recording of 
numerous parameters during evolutions; this 
continuous recording avoids tedious waiting 
for stabilised readings. We refer to the 
measuring of the polar in flight; this measure- 
ment when effected with engines cut and 
propellers feathered, supplies terms for direct 
comparison with the wind tunnel tests. 

The test method, thought out and 
developed at the Centre d’Essais en Vol de 
Bretigny, has been applied many times and 
in various conditions to an aircraft, the polar 
of which had also been determined in the 
wind tunnel with a non-powered model 
(scale: 1/20). The reliability of the results 
obtained is at least equivalent to wind tunnel 
possibilities. 

Figure 21 shows the aircraft which was 
the object of the tests: a twin-engined 
NC 702 Martinet powered by two Renault 
12 S 00 engines, developing 609 CV at take- 
off. Its total weight is 5,500 kilos. The mast 
on the fuselage and other externally mounted 
accessories correspond to the aircraft's special 
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Test equipment. 


Test equipment. 


Fig. 21. 
N.C.702 Martinet and its test equipment. 
N.C.702 No. 136. 


use as a flying laboratory. The polar has 
also been measured on an aircraft of the same 
type without any special external equipment. 

The evolution required from the pilot is 
shown diagrammatically in Fig. 22. The 
evolution is in two periods :— 

One period of acceleration started at the 
lowest possible speed, during which constant 
longitudinal trim is maintained in a steep 


dive and checked by means of a collimator 
trained on the visible horizon. 

A period of deceleration corresponding to 
the same flying methods at constant trim and 
applied to a second collimator setting giving 
a longitudinal trim with slight tail-heaviness. 

Thus, the polar is described twice in 
approximately three minutes. Figure 23 
shows examples of recording on photo- 
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Fig. 22 
N.C.702 
graphic paper tape 89 mm. wide, i.e. 3.5 in., Cy. — 
which can be set to run at the rate of 0.5 to c =tan (= +i) 


170 mm. per second, i.e. 0.02 to 15 in. per 
second. 

The main parameters recorded are: time, 
dynamic pressure, altitude, incidence and 
direction of the resultant acceleration, to 
which generally added normal 
acceleration, outside temperature, manifold 
pressures, engine r.p.m. and also, elevator 
position and forces thereon. The last few 
parameters will enable curves of stick forces 
and movement to be drawn identical to those 
which would be obtained in the stabilised 
condition, since there are neither velocities 
nor angular accelerations. 

The three quantities requiring extremely 
accurate measurement are: dynamic pressure, 
incidence and direction of the resultant 
acceleration. For these latter two angles, 
errors must be kept down to a twentieth of a 
degree. 

This accuracy is achieved by means of 
“pressure intakes” bored in the fuselage nose 
or in a leading edge and by using a multi- 
mirror, heavily damped, pendular level 
indicator giving a 13 mm. deflection, i.e. 
0.51 in. per degree. These two parameters 
give us the Cp by means of the formula: 
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where Cp is the drag coefficient, C, the lift 
coefficient, i the incidence and = the pen- 
dulum angle with the perpendicular to the 
incidence reference. 

Figure 23 shows an example of recording 
on which we find: the time, the dynamic 
pressure, the static pressure, the longitudinal 
trim, the incidence and pendulum angle, and 
below, points from checks of various flights, 
aimed at showing the order of magnitude and 
dispersion of results. 

Figure 24 enables comparison to be made 
of the polars obtained for the Martinet 136 
with externally mounted equipment, the 
normal Martinet 282 and the wind tunnel 
model. 

We have applied the same flying methods 
to a large variety of cases, for example: to 
assess performances at most of the flying 
and engine conditions, with the essential 
advantage that the results obtained remain 
independent of vertical motions of the 
atmosphere; it follows that this method gives 
an accuracy which cannot be hoped for when 
using classical methods. 

When it is desired primarily to find the 
vertical velocity v that will be observed at 
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the same configuration in stabilised flight, the 
following equivalent relation is used : — 
v=V sin (i+7z) 

The polars in flight which were left until 
now on one side, because of the magnitude of 
possible errors and high cost, are now com- 
non practice in flight testing due to standard 
measuring equipment. 

We have just seen a case where the 
measurements effected in slowly changing 
fying conditions have been utilised. The 
recording equipment now in use was not 
designed for the study of phenomena which 
have a fundamental frequency greater than 
100 cycles per second. These will probably 
call for cathode oscillographs. However, the 
equipment available enables us to study the 
fastest aircraft manoeuvres, such as deck- 


landing on an aircraft-carrier. The example 


which follows is related to this latter 
manceuvre. 

It is not possible as yet to carry out with 
prototype models tests which are entirely 
similar to actual deck-landing. In this field, 
it appears that the flight test can be used for 
directing and preparing wind tunnel testing. 
Such tests are already frequently made on 
powered models, with rapid engine-thrust 


Fig. 25. 
The S.B.D.5 Dauntless and its test equipment. 


variation. This was one of the aims of the 
simulated, or rather “idealised” deck-landing 
manceuvre experimented with at the Centre 
d’Essais en Vol. 

The manceuvre was suggested by the 
Service Technique Aéronautique to ascer- 
tain certain qualities of a deck-landing 
aircraft. It consists in taking the aircraft in 
level flight at the speed for final approach. 
The engine is cut and then, immediately, the 
controls are locked. Continuous recording 
of the various parameters which determine 
the flight path and the altitude of the aircraft, 
against time, while the controls remain 
locked, gives numbers which will be the basis 
for comparisons between various deck- 
landing aircraft. 

The first tests, which we present here, were 
made in 1947 with the $.B.D. 5 Dauntless, a 
two-seater deck-landing aircraft well known 
for its qualities. 

Figure 25 shows this aircraft and the 
recording equipment, which is located facing 
the observer; Fig. 26 gives an example of 
recording and some of the curves deduced 
from six tests made in succession; these are 
mean curves which show only a very slight 
dispersion. 
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Fig. 26. on 
Curves from records during idealised deck landing. é 


Fig. 27. Angu 


Spin case selected in view of a comparison between tunnel and flight tests. Loss 
C.G. position—forward. tud 
Elevator control—in the middle. Spl 
Rudder control—fully starboard. luti 
Ailerons control—fully port. 
Flaps up and cockpit closed. ; 
The two-seater aircraft Fokker §.11 Instructor on the aerodrome of Bretigny. 
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When the aircraft on test has longitudinal 
instability when deck-landing, as is the case 
with the Seafire, the method, obviously, will 
have to be modified. Elevator control is then 
ket to the pilot who has to maintain constant 
jongitudinal trim. 

There is no doubt that direct recording in 
the aircraft itself, which does not necessitate 
the presence of observers on the ground, 
renders the tests far more economical for, 
even with radio communication, the fewer the 
people required for preparing the test, the 
hetter the tests are likely to proceed. 

In spite of these setbacks, measurements 
carried out from the ground or from an 
accompanying aircraft are resorted to when- 
ever they offer greater precision, or when 
they are essential in conjunction with classical 
measurements. As an example of these 


| measurements, we will give the reproduction 


of a spin with a Fokker, type II aircraft, 
made in order to determine the aircraft’s 
geometrical characteristics, particularly the 
instantaneous axis of rotation of the aircraft. 
At the same time, we were able to compare 
the results obtained with the tunnel tests, 


which we summarise. 


The model of the two-seater Fokker S.II 
Instructor was built and tested at 1’Institut 
de Mécanique des Fluides de Lille 
(ON.E.R.A.). It had a span of 0.80 metres 
and a weight of 0.500 kilos which was 
determined by the rules of similarity for 
altitudes ranging from 500 to 2,000 metres. 

Spin and spin exit characteristics were 
assessed in the two metre diameter vertical 
wind tunnel through photographs which were 
taken at the rate of 60 per second on a 
35 mm. film. The parameters were measured 
on the pictures where the aircraft is shown 


Tunnel Flight test Flight test 


Characteristics IMFL. Bretigny NLL. 
values Units No. 24 No. 2D No. 6 
Average 
altitude metres 2000 1200 2000 
Longitudinal 
attitude degrees —-46 -46 
Transversal 
attitude degrees —6 —2 —4 
Rate of metres per 
descent second —27 -28 -33 
Angular revolution 
speed persecond 0.53 0.45 
Loss of alti- 
tude per 
Spinrevo- metres per 
lution revolution 51 67 73 
Fig. 29. 


Principal characteristics of the spin. 


with its longitudinal axis parallel or perpen- 
dicular to the axis of the camera. The model 
is kept at the same height as the camera axis 
by adjustment of the air current velocity. 

The actual aircraft, which has a span of 
11 metres and weighs 1,100 kilos,* made four 
spins at Bretigny, flown by the Fokker chief 
pilot, Mr. Sonderman. 

The characteristics of these actual spins 
were assessed by means of films obtained 
with two Askania cinetheodolite sets 1,300 
metres apart. Pictures were taken at the rate 
of two per second and the checking methods, 
based on descriptive geometry, provided 
figures accurate within a metre for heights 
and within a degree for angles of trim. 

Earlier, detailed tests had been made in 
Amsterdam by the National Lucht Vaart 
Laboratium with the recording of the various 
parameters by photographing an instrument 
panel. The satisfactory agreement between 
the tests made in Lille and those in Amster- 
dam has been confirmed by the Bretigny tests. 
We illustrate this agreement with an example 
chosen from the nine spins recorded by the 
N.L.L., the four tests at Bretigny and the 
nine tests filmed at Lille. Each spin was of 
ten revolutions in flight and of four in the 
tunnel. 

Figure 27 shows the Fokker aircraft and 
indicates the test conditions. Fig. 28 shows 
the model in the Lille wind tunnel and the 
aircraft in flight over Bretigny; the latter 
photograph was taken with one of the cine- 
theodolites and indicates height and azimuth. 

Figure 29 summarises the mean values 
obtained for the most important parameters, 
in flight and in the tunnel. A few curves 
representing a spin filmed at Bretigny are 
shown on Fig. 30. 

In this lecture we have been able to give 
only very briefly the lines of general policy, 
illustrated by a few examples. We hope, 
however, that you will have seen our 
endeavour to conduct these experiments in a 
scientific manner. 

In the field of aeronautics measurement 
work whether it be carried out on the ground 
test installations or in the actual aircraft, 
requires a thorough knowledge of general 
physics, electricity, electronics and even 
metallurgy. Fruitful co-operation in these 
fields already exists between British tech- 
nicians and French specialists. As examples 
of this co-operation, we will only quote the 
adoption in France of the British Measuring 


*For conversion, one kilo=2.205 lb. 
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Fig. 28. 


The model of Fokker S.11 aircraft in the wind tunnel of Lille. 
An example of picture taken by the Askania cinetheodolite. 
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Methods of Control Column Forces per “g,” 
or the use of British magnets in our recorder 
galvanometers. 

It is in the name of all French technicians 
that we assure our British colleagues of our 
desire to help them to the best of our ability. 
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THE VOTE OF THANKS 


Sir Ben Lockspeiser: He was glad to have the opportunity of joining in the 
welcome to their French colleagues, and in particular to thank on behalf of the 
Society M. Brocard and M. Hussenot for the Second Louis Bleriot Lecture which they 
had just heard. 

He had an added interest in listening to that lecture because, as the authors 
quite rightly and properly pointed out in their opening remarks, France was paralysed 
during the Occupation, and their achievements had to be judged against the back- 
ground in which the French scientists had been working. They had worked under 
very difficult conditions indeed. He thought they would all agree that that 
achievement under those circumstances was a great thing, to be admired and for 
which they were thankful. 

It so happened that he had been in France at the period when science and other 
activities were at their lowest ebb for many years in French history. He had been 
in Paris on that memorable day of liberation, and saw the people, every man, woman 
and child in Paris it seemed, come out into the light and sunshine. That was a 
memorable event in the history of France and anyone who witnessed it, as he had 
done, must have been deeply moved. He did not say that he actually wept—no 
Englishman dared do that, least of all in France—that would be letting the side 
down, but it made an unforgettable impression on his mind. He remembered passing 
the newly gilded statue of Joan of Arc, resplendent in the sunshine. As often 
happened on occasions like that some trivial thing passed through the mind and he 
remembered recalling that St. John Irvine had said that the French were the only 
nation which could be trusted with a pot of gold paint. : 

He was in Paris on the invitation of the French authorities to give them a 
résumé of the achievements which had been made during the war in aeronautics, 
and to talk matters over with them. Having done this to the best of his ability they 
then showed him the condition of their wind tunnels and laboratories. The wind 
tunnels were shells; they were useless because the power plants had been removed 
and in the laboratories much of the equipment had gone. What was left convinced 
him that they would have to start from scratch. But the spirit was there and he 
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saw people in the laboratories making sketches and drawings of those inimitable 
pieces of small mechanism that they all recognised could only come from France. 

There was one ray of sunshine; it so happened that the French Zone in Germany 
contained a very fine wind tunnel, so they had a chance of getting their own back 
and something more. That wind tunnel, half completed, was intended to be 
harnessed to a water-fall. The fans were driven direct, through Pelton wheels and 
were capable of maintaining a 24 ft. diameter stream of air at a Mach number of one. 

When he heard of it he said “That is a thing you must really go out for because 
a tunnel like this would be unique” and he was glad to be able to say that in the 
ensuing years they were able to give their French colleagues considerable technical 
assistance and help in procuring the machinery, part of which was in the French 
Zone and part in the British Zone. In those days it was not easy to get passage of 
material from one zone to another. He believed that tunnel had made extremely 
good progress. The French had found their own water-fall in the Rhone valley and 
there was good hope that the fans would begin ticking over in about 1950. It would 
be a very great accession of strength to the equipment for French aeronautical science 
and he thought they might hope to share in some measure the benefits of the tunnel 
that had come to their French colleagues and that thereby their collaboration in 
technical matters would be enhanced. 

They were proud to be associated with their French colleagues on the occasion 
of the Louis Bleriot Lecture, and in particular to record their appreciation that as 
a Society they were associated with the name of Bleriot. Bleriot was the kind of 
man of whom there never could be many in any one generation. He was one of 
a small band who really brought home to the consciousness of the people in this 
country that flying was a serious business. The Wright Brothers had flown in 
America, but America in those days was a long way away. He well remembered 
the event of Bleriot flying the Channel. It was one of those dramatic events which 
made people realise that an entirely new force had entered their World which was 
going to affect the country and Europe, economically, politically and militarily. 
It was the time, he believed, when the phrase was born, “ We are no longer an island.” 

They should, also, recognise that Bleriot was a great engineer. He never saw 
the model of his aeroplane, or a picture of it without admiring the simplicity of it— 
the design for strength with lightness—and he thought it right to say that with the 
engines and materials available in those days no better design could have been made. 

Bleriot was not only a great engineer, he was also a man of great courage and 
he thought they recognised in the work which their French colleagues had been doing 
since the liberation, and in particular in their lecturers MM. Brocard and Hussenot, 
the same quality of courage that Bleriot showed in his flight. ; 

He proposed, he was sure with their hearty concurrence, a vote of thanks to 
M. Brocard and M. Hussenot for the lecture, which they valued not only for itself, 
but for the spirit behind it. They were heartened by it and he was sure their 
lecturers must also be heartened by the achievements which they had made under 
such great difficulties. 


Dr. H. C. Watts: He would have taken it as a great honour to be asked to second 
a vote of thanks at any of those epoch-marking lectures, but to do so on that occasion, 
at the first of the Louis Bleriot lectures in Great Britain was a great personal pleasure. 

The lecture was not open to discussion and Sir Ben had expressed their thanks 
fully. As he watched the slides and saw and read the names of the French experi- 
mental establishments there was one feeling uppermost in his mind. That was a 
feeling of thankfulness that in spite of everything that France had suffered the French 
genius for aerodynamical testing still lived and flourished as it did in the years 
immediately following the exploit of Louis Bleriot. As he listened to the lecture 
his thoughts went back nearly 40 years to a day when, as a raw youngster in his 
first job at what was now the Bristol Aeroplane Company, he was told to design 
a propeller. Without any data and with no text books, it was to the proceedings of 
Les Academie des Sciences that he turned to find the Dregweicki Theory and to 
the work of M. Eifel, published in the same year as the Cross-Channel flight, that he 
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found the aerofoil data with which to apply it. In those days he was most grateful 
to the early French experimenters and he was grateful on that occasion to MM. 
Brocard and Hussenot for a paper which had worthily upheld the traditions those 
early experimenters had established. Seconding the vote of thanks he not only 
expressed their sincere thanks, but acknowledged his own debt to “ French practical 
aerodynamical methods.” 


Following the Second Louis Bleriot Lecture a Council Dinner was held at which 
the following were present : — 


M. Jacques Allez, President, Aero Club de France; L. S. Armandias, Esq., Correspondent 
in Great Britain for Association Frangais des Ingénieurs et Techniciens de l’Aéronautique. 

Air Commodore F. R. Banks, C.B., O.B.E., F.R.Ae.S., First Louis Bleriot Lecturer 1948; 
M. P. Baudet, Minister Plenipotentiary and First Counsellor of the French Embassy; Lord 
Brabazon of Tara, M.C., F.R.Ae.S., President, Royal Aero Club, Past President, Royal Aero- 
nautical Society; M. Louis Breguet, Membre d'honneur de |’ Association Frangais des Ingénieurs 
et Techniciens de |’Aéronautique; M. J. M. R. Brocard, Co-lecturer; Sir John S. Buchanan, 
C.B.E., A.M.I.Mech.E., F.R.Ae.S., Vice-President, Royal Aeronautical Society; Major G. P. 
Bulman, C.B.E., B.Sc., F.R.Ae.S., Director for Construction of Research Facilities, Ministry of 
Supply (Air), Member of Council. : 

S. Camm, Esq., C.B.E., F.R.Ae.S., Director and Chief Designer, Hawker Aircraft Ltd., 
Member of Council; A. V. Cleaver, Esq., A.R.Ae.S., Special Projects Engineer, de Havilland 
Engine Company, Member of Council; Air Marshal Sir W. Alec Coryton, K.B.E., C.B., M.V.O., 
D.F.C., Controller of Supplies (Air), Ministry of Supply; Dr. H. Roxbee Cox, D.LC., B.Sc. 
F.I.Ae.S., F.R.Ae.S., President, Royal Aeronautical Society, Chief Scientist, Ministry of Fuel and 
Power. 

Ingénieur en Chef N. Daum, Representing Le Directeur Technique et Industriel; Sir Charles 
Darwin, K.B.E., M.C., D.Sc., F.R.S., Director, National Physical Laboratory; Dr. G. P. Douglas, 
O.B.E., M.C., D.Sc., F.R.Ae.S., Head of Aerodynamics Department, Royal Aircraft Establish- 
ment, Member of Council. 

W. S. Farren, Esq., C.B.E., F.R.S,. F.R.Ae.S., Technical Director, A. V. Roe and Co. Ltd., 
Member of Council; Sir A. H. Roy Fedden, M.B.E., D.Sc., M.I.Mech.E., M.S.A.E., F.R.Ae.S., 
Director of Research and Development, Board of Leyland Motors, Past President, Royal Aero- 
nautical Society. 

Francis Grierson, Esq., Representing Lord Rothermere. 

Professor A. A. Hall, M.A., F.R.Ae.S., Zaharoff Professor of Aviation, University of 
London, Member of Council; Mr. S. Scott-Hall, A.C.G.I., M.Sc., D.I.C., F.R.Ae.S., Principal 
Director of Technical Development (Air), Ministry of Supply, Member of Council; Ingénieur 
en Chef de l’Air F. Hussenot, Co-lecturer. 

M. J. Jarry, A.F.R.Ae.S., President, Association Francais des Ingénieurs et Techniciens 
de l’Aégronautique; E. T. Jones, Esq., O.B.E., M.Eng., F.R.Ae.S., Director of Instrument Research 
and Development, Ministry of Supply, Member of Council. 

Sir Ben Lockspeiser, M.A., F.C.S., F.R.Ae.S., Chief Scientist, Ministry of Supply, Member 
of Council; Captain René Lucien, F.R.Ae.S., l’Association Francais des Ingénieurs et Techniciens 
de 1Aéronautique. 

P. G. Masefield, Esq., M.A., F.R.Ae.S., Assistant to the Chairman and General Manager. 
British European Airways, Member of Council; Sir Francis McClean, A.F.C., F.R.Ae.S., Vice- 
President, Royal Aero Club. 

Sir Arnold Overton, K.C.M.G., M.C., Permanent Secretary to the Ministry of Civil Aviation. 

Sir Frederick Handley Page, C.B.E., F.R.Ae.S., Managing Director, Handley Page Ltd., 
Past President, Royal Aeronautical Society; The Rt. Hon. Lord Pakenham, Minister of Civil 
Aviation; W. G. A. Perring, Esq., F.R.Ae.S., Vice-President, Royal Aeronautical Society: 
Captain J. Laurence Pritchard, Hon.F.I.Ae.S., Hon.F.R.Ae.S., Secretary Royal Aeronautical 
Society. 

Sir Alliott Verdon-Roe, Hon.F.I.Ae.S., Hon.F.R.Ae.S., Chairman, Aviation Developments 
Ltd.; N. E. Rowe, Esq., C.B.E., B.Sc., D.LLC., F.R.Ae.S., Controller of Research and Long Term 
Development, British European Airways, Vice-President, Royal Aeronautical Society; Sir 
Archibald Rowlands, G.C.B., M.B.E., Permanent Secretary to the Ministry of Suoply. 

General P. Stehlin, O.B.E., Air Attaché to the French Embassy. Wy wt: 

M. E. Tendron, Secretary, l'Association Francais des Ingénieurs et Techniciens de l’Aéro- 
nautique. ‘ : 

Captain C. F. Uwins, A.F.C., O.B.E., F.R.Ae.S., Managing Director (Aircraft), Bristol 
Aeroplane Co. Ltd., Honorary Treasurer ,Royal Aeronautical Society. : 

Dr. H. C. Watts, M.B.E., D.Sc., M.Inst.C.E., F.R.Ae.S., Consulting Engineer, Member of 
Council; L. A. Wingfield, Esq., M.C., D.F.C., A.R.Ae.S., Solicitor to the Roya! Aeronautical 
Society; Air Cmdre. Sir Frank Whittle, C.B., C.B.E., M.A., D.Sc., F.R.S., Hon.M.1.Mech.E.. 
F.R.Ae.S., Technical Adviser on Engine Design and Production to the Controller of Supplies 
(Air), Ministry of Supply. 
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FLUTTER AND STABILITY 


by 


W. J. DUNCAN, D'Sc., F.R.S., F.R.Ae.S., M.I.Mech.E. 


The 768th Lecture was read before the 
Royal Aeronautical Society on 24th 
February 1949 at the Institution of Civil 
Engineers, London S.W.1. Dr. H. Roxbee 
Cox, D.LC., F.R.Ae.S., F.1.Ae.S., President 
of the Society, introduced the Lecturer, 
Professor W. J. Duncan, D.Sc., F.R.S., 
FR.Ae.S., M.I.Mech.E., who, in 1929, with 
his colleague Dr. Frazer had been responsible 
for R. & M., 1155, a great classical work 
which continued to influence flutter theory, 
and who was now Professor of Aerodynamics 
at the College of Aeronautics. Professor 
Duncan did certain experiments with some 
apparatus he had brought with him and 
explained particular flutter phenomena. 


INTRODUCTION 


The theme of flutter has not often been 
the subject of lectures before this Society, 
but the earliest of these“) was given as long 
ago as 1929 by my friend and former 
colleague Dr. R. A. Frazer. At that time 
a secure foundation for flutter theory had 
been laid and the principal methods for 
flutter prevention were already known. 
However, both theory and practice have 
developed much since 1929 and flutter has 
become of greater relative importance. The 
main reason for this is that flutter is a trouble 
which is most liable to occur in flight at 
high speeds, while speeds of flight have 


increased by a factor of three or more in the 
interval. 

Since the date of Frazer’s lecture other 
aeroelastic effects besides flutter have 
become important and it is no longer possible 
to regard an aircraft as rigid when develop- 
ing the theory of stability and control in 
relation to flight at high speeds. Further, 
the compressibility of the air can no longer 
be ignored in the theory of either flutter or 
general stability. Indeed the two subjects 
are becoming entangled to an ever increasing 
extent. 

My first aim in this lecture is to elucidate 
the subject of flutter by giving a strictly 
non-mathematical account of it, and with 
the emphasis on the physical causation. 
Needless to say, this approach is of limited 
utility and the serious student must study 
the mathematical theory. I feel that even the 
expert and the specialist gain from a non- 
mathematical exposition and I hope that my 
treatment will be found illuminating. My 
second aim is to discuss some aspects of 
general aircraft stability and of control. Here 
I have been frankly selective and I do not 
pretend to give anything resembling a general 
discussion. Lastly, I have hazarded a few 
opinions and stated fewer facts about some 
urgent problems of the day. Here, even 
more than elsewhere, I hope to learn from 
my audience. Throughout the lecture I shall 
try to keep the light touch. 


PART I 
A NON-MATHEMATICAL SKETCH OF FLUTTER 


WHAT IS FLUTTER? 


_ Flutter may be defined as an oscillatory 
instability occurring in an aircraft in flight, 
where the flexibility of the structure plays 
an essential part in the instability. As we 
now know well, the flexibility of the structure 
always influences the general “classical ” 
stability of an aircraft and sometimes 
profoundly, but there is usually a recognis- 


able distinction between the classical aircraft 
instabilities and flutter. As an instance, we 
do not class “snaking” as flutter, for the 
essential freedoms are rudder motion about 
the hinge and the general bodily freedoms 
of the aircraft, while distortion of the 
structure is secondary. But if the flexibility 
of the fin, for example, became great 
enough to influence the snaking to an 
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important extent, the distinction from flutter 
would cease. 

This example illustrates the fact that the 
boundary between flutter and _ general 
instability, at one time sharp and clear, is 
tending to become vague, as has been 
emphasised by Professor Collar in_ his 
paper) ““ The Expanding Domain of Aero- 
elasticity.” As a consequence, it is now 
necessary, in my opinion, for all scientists 
and technicians concerned with general 
stability to be well acquainted with flutter 
and the techniques of flutter theory; like- 
wise the flutter specialists must be well 
versed in general stability theory. 

The phenomena of typical flutter can be 
well illustrated by experiments with a 
cantilever wing mounted in a wind channel 
with rigid support at the root. We may 
suppose that the aileron is locked to the 
wing, whose essential freedoms are bending 
and twisting; alternatively we may suppose 
that the aileron can move on its hinge under 
the constraint of a spring, when the wing 
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may be torsionally rigid. In still air and 
when the wind speed is low the wing is quite 
stable, as shown by the fact that a 
“ damped ” oscillation, ultimately disappear. 
ing, follows any disturbance imposed on the 
wing. But, unless the wing has been 
designed to be flutter-free, a critical wind 
speed will be reached where a steady 
oscillation can just maintain itself, following 
a disturbance. At a still higher wind speed 
any slight disturbance will initiate an 
oscillation which grows in amplitude until, 
possibly, fracture of the wing occurs. Two 
facts of cardinal importance can easily be 
demonstrated by experiments of this kind: 


(a) The critical speed for flutter can be 
profoundly influenced by changes in the 
distribution of mass. For example, 
flutter of the flexural-aileron type can 
be prevented by arranging that the 
centre of mass of the aileron lies on or 
near the hinge axis. This is called mass- 
balancing the aileron. 


Fig. 1. 
The flutter engine. 
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Fig. 2. 
Phase-changing gear of flutter engine. 


() The critical flutter speed can be raised 
- by suitable stiffening of the system. For 
example, with the torsionally rigid 
model, the critical speed can be raised 
by a sufficient increase in the stiffness 
of the spring constraining the aileron 
movement. 


THE PHYSICAL CAUSATION 
FLUTTER 


Experiments with models, such as those 
described in the last section, serve to demon- 
strate that flutter is entirely distinct in 
nature from forced or resonance oscillation 
since the flutter occurs in the absence of any 
periodic external agent. Thus we have 
to recognise that a fluttering wing is a kind 
of air-engine which abstracts energy from 
the passing air stream. The problem is to 
understand how this happens. 


We shall reach the answer to this question 
by means of experiments with a piece of 
apparatus’*) which Mr. D. L. Ellis and I 
have called the “‘ flutter engine.” This con- 
sists of a rigid aerofoil mounted on a spindle 
running spanwise and able to move in pitch 
about the spindle so as to change its angle 
of incidence (see Fig. 1). The spindle itself 
can move in roll about a fixed axis at the 
wing root, so the aerofoil has two degrees of 
freedom, namely roll (corresponding to 
flexure of an elastic cantilever wing) and 
pitch (corresponding to torsion of an elastic 
wing). Oscillatory motions of these two 
kinds are imparted to the wing by a pair of 
cranks with connecting rods suitably linked 
to the wing and transverse spindle. The 
cranks are on a single shaft, but the angle 
between the cranks can be varied by means 
of a special gear which is an adaptation of 


OF 


the ordinary differential gear (see Fig. 2). 


In this way the phase difference between the 
aerofoil oscillations in roll and in pitch can 
be changed at will. 

Let the cranks be set so that the pitching 
oscillation lags 90 degrees behind the roll. 
This implies, for example, that the positive 
change of incidence is a maximum at the 
middle of the upstroke in roll and that the 
negative change of incidence is a maximum 
at the middle of the downward stroke in 
roll. It results from this that the lift on the 
aerofoil does positive work upon it during a 
complete cycle and if the wind speed is high 
enough, the engine will run under its own 
power. But when the pitch and roll are in 
phase the net work done by the lift is zero 
and the engine is unable to run under its 
own power. The diagram (Fig. 3) will make 
clear how the work done varies throughout 
the cycle in the two cases. 


We draw the conclusion from this 
experiment that flexural-torsional flutter of a 
wing can occur when the torsional oscillation 
lags sufficiently behind the flexural oscilla- 
tion and this inference is supported by 
experiment and by detailed mathematical 
analysis. Just at a critical flutter speed, 
where a flutter of steady amplitude can 
persist, the energy fed into the wing per 
cycle by the lift just balances the energy 
dissipated by the various damping actions 
during the cycle. 

We have still to explain the occurrence 
of critical flutter speeds. It is not possible 
to provide a complete explanation without 
the help of mathematics, but the broad 
explanation is given by two facts:— 


(a) The angle of lag of the torsion behind 
the flexure depends on the wind speed 
(see, for example, Fig. 4). 


(b) For a given angle of lag the energy input 
and the dissipation depend on wind speed 
in different ways (see Fig. 5). Other 
things being unchanged, the energy 
input is proportional to the square of 
the wind speed, whereas the dissipation 
is proportional to the wind speed. This 
goes far to explain why flutter occurs 
usually at high speeds. (See Appendix). 


Another factor which influences the energy 
balance and which varies with wind speed is 
the ratio of torsional amplitude to flexural 
amplitude (see, for example, Fig. 6. Figs. 
4, 5 and 6 refer to the same wing, of which 
the particulars are given in Reference 3). 
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TORSIONAL DISPLACEMENT LAGGING 90° BEHIND FLEXURAL DISPLACEMENT 


(TORSIONAL DISPLACEMENT IN PHASE WITH UPWARD FLEXURAL VELOCITY] 
FLEXURAL AND TORSIONAL WING DISPLACEMENTS DURING A COMPLETE CYCLE. 


The diagrams show the wing tip in end view at intervals of one eighth of a period, the relative 
wind blows from right to left. 


Fig. 3. 
Diagram showing lift force and work done during a cycle. 


The conclusions about the causation of 
flutter which we have drawn from our very 
incomplete, but essentially sound, arguments 


are important because they apply, in 
principle, to all varieties of flutter and to 
oscillatory aerodynamic instabilities in 
general. 


Experiments on flexural-torsional flutter of 
a cantilever wing and with the flutter engine 
have served to focus attention on the force 
tending to produce flexure which is brought 
into play when the wing twists. This is a 
particular instance of the coupling of two 
degrees of freedom and here the coupling is 
of aerodynamic origin. We say that two 
kinds of motion A and B are coupled when 
the existence of A tends to induce B or 
when the existence of B tends to induce A 
and, in general, both of these coupling 
actions will be present. Usually the term 
is restricted to mean direct coupling of A 
and B; two motions which have no direct 
coupling may yet influence one another 
because both are coupled to some third 
motion. Couplings are of three types, 
namely, inertial, elastic and aerodynamic 
and in order to isolate motion A from 
motion B it would be necessary to abolish 
all the couplings. 
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The experiments with the flutter engine 
suggest and mathematical analysis confirms, 
that flutter is essentially caused by aero- 
dynamic coupling of motions which could 
occur independently, i.e. be uncoupled, in the 
absence of aerodynamic actions. There is 
one exceptional case, namely stalling flutter, 
which only occurs for angles of incidence 
near the stall and which is an oscillatory 
instability essentially in a single freedom, 


involving change of sign of a damping | 


coefficient **), Coupled flutter is sometimes 
called classical flutter. 


THE NATURE OF THE FLUTTER 
PROBLEM 

The flutter problem presents itself alike 
to the scientist, to the designer and to the 
airworthiness authority under three main 
headings:— 

(i) Understanding and analysis 
(ii) Prediction of critical flutter speeds and 


of flutter characteristics i 


(iii) Flutter prevention. 


The outcome of understanding and analysis 
is a workable theory which becomes the 
basis alike for prediction and for the rational 
design of effective measures of prevention. 
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We must therefore consider this workable 
theory at some length but, in this place, 
without avoidable technical intricacy. 

The theory of flutter is the theory of a 
particular kind of mechanical oscillation 
complicated by aerodynamic actions. Hence 
the theory must have three main aspects 


(a) Kinematic 
(b) Mechanical 
(c) Aerodynamic 
These are considered separately in the 
following sections 


THE KINEMATICS OF FLUTTER 


We may legitimately simplify our problem 
somewhat by confining our attention to the 
flutter of constant amplitude which occurs 
just at a critical flutter speed. Thus we 
shall assume that in the flutter all parts of 
the affected structure have simple-harmonic 
motions of common frequency. 

The main kinematic problem, which is in 
many instances the most difficult portion of 
the whole flutter problem, is to decide 
a priori what kinds of movement are 
important in the flutter. In other words, the 
problem is to find a rational and reliable 
means for the selection of the degrees of 
freedom to be used in the calculations. 

In some of the simpler kinds of flutter this 
question can easily be answered as the result 
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Fig. 4. 
Variation of angle of lag with wind speed. 


of theoretical studies and numerous experi- 
ments. For instance, in the case of an 
unswept cantilever wing not carrying engines 
or other large concentrated masses and 
having rigid support at the root, it suffices 
to treat the wing as having just two degrees 
of freedom, one for flexure and one for 
torsion. This amounts to treating the wing 
as “semi-rigid,” i.e. as having fixed modes 
of distortion in flexure and in torsion. We 
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Variation of amplitude ratio with air speed. 
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can easily generalise the notion of semi- 
rigidity and suppose the wing to have two 
or more freedoms in both flexure and in 
torsion”). Completely successful applica- 
tions of this method have been made to the 
flexural-torsional flutter of cantilever wings 
carrying engines: * 7), 

Another treatment which has much to 
commend it consists of assuming that the 
motion in the flutter is compounded of a 
few of the natural free or resonance modes 
of the structure. As will be pointed out 
later, this results in eliminating a large 
number of the couplings between the distinct 
kinds of motion. 


The freedoms or kinds of motion already 
listed are not, in general, in themselves 
sufficient. An additional freedom must 
certainly be introduced. for each control 
surface, unless it can be shown that it takes 
no part in the flutter considered. This will 
be true when the control surface has really 
effective irreversible operation, or when it 
is not coupled to the other motions. An 
example of the last is that a single 
symmetrically placed rudder takes no part 
in symmetrical wing flutter. The other 
important freedoms not already mentioned 
are the general “rigid-body freedoms” of 
the aircraft as a whole, which require 
separate discussion. 

It is obvious that if the wings of an air- 
craft flutter antisymmetrically, so that the 
maximum upward bend of the starboard 
wing and the maximum downward bend of 
the port wing occur simultaneously, there 
will be a periodic rolling moment applied to 
the fuselage which will therefore have a 
periodic rolling motion. Theory and experi- 
ment have shown that the rolling freedom 
of the fuselage may have a considerable 
influence on the critical speed for anti- 
symmetric wing flutter, and in general this 
freedom must be retained in the theory of 
such flutter. Likewise the general freedoms 
of pitch and of translation normal to the 
plane of the wings must, in general, be 
retained in the theory of symmetric: wing 
flutter‘ 

We may summarise our conclusions by 
saying that the following freedoms must be 
used in calculations on flutter: 

(a) Appropriate resonance modes or other 
simple modes of structural deformation. 

(b) The relevant rigid body freedoms of the 
aircraft as a whole. 


(c) The relevant control movements. 
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This still leaves much uncertainty in the 
selection of freedoms and the choice will be 
guided by experience of similar problems 
When a problem of a new kind has to be 
faced we may begin by selecting a likely set 
of freedoms and calculate the lowest critical 
flutter speed and the characteristics of the 
flutter. Then we may bring in, one at a 
time, other freedoms which are judged to 
have a possible share in the flutter and find 


whether their introduction appreciably 
influences the critical speed and so on‘, 
This procedure, if intelligently and 
thoroughly carried through, will yield 
reliable results but it is exceedingly 
laborious. 


I have recently suggested elsewhere() 
that the additional freedom for trial should 
be selected on the following basis. A 
proposed freedom should be retained if the 
energy balance at an already calculated 
critical speed is sensitive to its introduction, 
unless it can be shown that the amplitude 
of oscillation in the new freedom is very 
small. The factors leading to smallness of 
amplitude are smallness of the couplings 
affecting the new freedom and high direct 
impedance to motion in that freedom at the 
flutter frequency. 

The selection of suitable freedoms for use 
in calculations about flutter is the choice of 
a conceptual scaffolding. Obviously the 
physical happenings are independent of this 
scaffolding, but the actual motion, in 
general, involves a very large (infinite) 
number of freedoms. Hence, when we 
choose a finite (small) number of freedoms, 


in general we are distorting the facts to some | 
extent and so the results of the calculation © 


are not independent of the choice. 


THE MECHANICAL ASPECT OF 
FLUTTER 


The mechanical effects in flutter may be 
classed as inertial and elastic, the latter 
including structural damping effects. Clearly 
the resistance of the structure to deformation, 
or to be precise its stiffnesses, plays a vital 
part in flutter. Subject to certain reasonable 
simplifications, it can be shown that the 
critical speed for flutter is raised by the 
factor nm when all the relevant stiffnesses 
are increased by the factor n. As a conse 
quence, the current airworthiness require 
ments for wings include a criterion for a 
minimum permissible stiffness in torsion. 
The emphasis is placed on the torsional 
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giffness since detailed calculations show 
that the flexural stiffness of an unswept 
cantilever wing is relatively unimportant in 
relation to flutter. 

The inertias of the aircraft are just as 
important as its stiffnesses, although perhaps 
less obviously so. In the first place the 
inertias, in conjunction with the stiffnesses, 
determine the natural frequencies of the 
gystem when oscillating in still air. A 
proportionate increase of all the material 
densities of the system reduces these 
frequencies, likewise the flutter frequency. 
Now it follows from the equation of energy 
that, other things being the same, an increase 
of frequency brings with it a proportional 
increase in the energy dissipated per cycle, 
but does not change the energy influx per 
cycle (see the Appendix). Thus an increase 
of frequency is a stabilising influence. 
Indeed, the relation between natural 
frequency and critical flutter speed is so close 
that Kiissner''') proposed a relation of the 
t 
V.=Kne, 
for flexural-torsional flutter, where K is a 
non-dimensional coefficient, n is the natural 
frequency of the “ dangerous” wing oscilla- 
tion in still air and c,, is the mean wing 
chord. Although K is by no means 
constant"*: 4) its variations, for wings which 
are not mass-balanced, are not very wide. 


The importance of natural frequency has 
also been emphasised in Pugsley’s simplified 
theory of flexural-torsional flutter’*), in 
which certain aerodynamic coupling coeffi- 
cients are neglected. It follows from what 
we have already said that a general increase 
of density, while other things remain con- 
stant, tends to reduce the critical flutter 
speed. However, detailed calculation shows 
that as the density rises the critical flutter 
speed ultimately approaches a_ constant 
asymptotic value’*) and this can be 
explained by the phase differences tending 
towards zero"). In the asymptotic condi- 
tion the “equivalent” critical flutter speed 
is independent of air density. 


Mass distribution has a profound influence 
on flutter and it is well known that the 
Participation of an aileron in wing flutter 
can be prevented by mass-balancing it. This 
Means, in essence, reducing to zero the 
Inertial coefficients which couple the motion 
of the aileron on its hinge to the bending and 
twisting of the wing. In practice this is 
achieved with sufficient precision by arrang- 


ing that the centre of mass of the aileron lies 
on or a little ahead of the hinge line, but 
the balance masses should be distributed 
along the span of the aileron and, ideally, 
should be large enough to counterbalance 
the virtual products of inertia. These are 
attributable to the motion impressed on the 
surrounding air by the oscillating aileron. 

Mass balance of hinged controls is a 
generally effective means of preventing them 
from taking an essential part in flutter. The 
spring tab is a special case and has peculiar 
pitfalls. The first of these is the “ limiting 
length ” of the balance arm, first discovered 
by Frazer and Jones"). I cannot enter 
fully into this question here and shall content 
myself with stating that, when the limiting 
length of the arm is exceeded, the balance 
mass increases, instead of decreasing the 
products of inertia of the tab and this 
follows from quite straightforward kine- 
matics'*), It is now recognised that effective 
mass-balancing of both tab and aileron is 
difficult and tricky and the favoured method 
of preventing flutter of the tab is to make it 
extremely Jight, in accordance with what has 
been said about the advantage of low 
density"). The aileron itself should be 
mass-balanced, but with a very light tab this 
is not difficult. 


It is of historical interest that the value 
of mass-balance of the aileron’ in 
suppressing flutter was pointed out by von 
Baumhauer and Koning’) as long ago as 
1923. 


THE AERODYNAMIC ASPECT OF 
FLUTTER 


In the theory of flutter we use the method 
of representing aerodynamic actions which 
was introduced by Professor G. H. Bryan‘??? 
in the discussion of general aircraft stability. 
Our main interest lies in the question 
whether flutter can be initiated by some 
small disturbance, such as a slight air bump 
or the movement of a control by the pilot, 
so the problem is to determine the stability 
of a small motion of deviation from a steady 
state of equilibrium. If the small motion is 
unstable it will grow into a large one to 
which our theory does not apply, but the 
matter of importance from our point of view 
is whether the small motion is stable or not. 
Since the deviations from the state of 
equilibrium are small, following the method 
introduced by Bryan, we may take the con- 
sequent changes of the aerodynamic forces 
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and moments to be linear functions of the 
displacements, velocities and accelerations of 
deviation and the coefficients of these 
quantities in the expressions for the forces 


and moments are called aerodynamic 
derivatives. 
In the treatment of general aircraft 


stability it is usual and often legitimate to 
calculate the aerodynamic derivatives as 
though each transient state of the system 
were a_ steady state. Derivatives so 
calculated are called quasi-static and their 
values can be estimated from the results of 
wind tunnel tests in which the model is 
stationary. However, quasi-static derivatives 
are not applicable when the changes of 
circulation round the wings, and the like, 
are relatively rapid for, in such circum- 
stances, vorticity is deposited in the wake of 


sufficient intensity to induce important 
velocities of downwash at the lifting 
surfaces. Analysis of the problem‘'**) shows 


that the critical quantity is 
Va 

where c is the chord, V the true air speed, 
a the amplitude of the incidence change and 
x the amplitude of the time rate of change 
of incidence (in unidirectional deviations we 
should substitute instantaneous values for 
amplitudes). It appears that quasi-static 
derivatives can be used with an error of less 
than five per cent. when 


x<0.02. 


Now for all ordinary flutters , is much 
greater than this bound (it is usually not far 
from unity), so quasi-static derivatives are 
definitely inapplicable to flutter. If we con- 
sider oscillatory deviations with frequency 
f we may easily show by dimensional 
analysis that the value of any given deriva- 
tive depends on the non-dimensional! 
frequency parameter 

2afc 


V 


which is equivalent to y, on Reynolds 
number 


R=Ve/v 
and on Mach number 

M=V/a 
where the symbols have their usual 
meanings. Measurements of the derivatives 


of flutter made in wind tunnels have amply 
demonstrated their dependence on_ the 
frequency parameter'?*: *4), 
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A theory of the aerodynamics of flutter 

has been worked out on the following 
basis: — 
(a) The flow is two-dimensional 
(b) The fluid is incompressible 
(c) The amplitude of motion is very small 
(d) The aerofoil is very thin 


(e) The velocity at the sharp trailing edge 
is always finite (Wagner’s generalisation 
of the hypothesis of Joukowski), 

(f) The Reynolds number is very great, 5 
the boundary layer and wake ar 
exceedingly thin. 


This theory has been worked out by a 
number of investigators and by a variety of 
methods, but the first success was scored 
by Glauert'’*) who based his work on 
Wagner’s general two-dimensional _ theory 
of transitory aerodynamic phenomena’*, 
Glauert obtained the aerodynamic deriva- 
tives for an aerofoil oscillating steadily in 
pitch, while Duncan and Collar'’) extended 
Glauert’s results to cover independent 
motions of normal translation and _ pitch, 
thus completing the set of two-dimensional 
derivatives for flexural-torsional flutter, and 
included divergent motions as well as simple 
harmonic oscillations. Theodorsen'**) first 
gave the derivatives or “air-load coefl- 
cients” for simple harmonic motion of a 
hinged flap and showed that the air-load 
coefficients for simple harmonic motion in 
two dimensions could be expressed by 
Hankel functions. A number of Continental 
investigators, notably  Kiissner'**) and 
Cicala®’), have contributed to this subject. 

Although the theory which has just been 
discussed has obvious limitations, there can 
be no doubt that it is soundly based and 
gives a rational explanation of the depend: 
ence of the derivatives on the frequency 
parameter. There is already a considerable 
volume of experimental evidence to show 


that the derivatives for the aerofoil proper | 


are fairly closely predicated by the theory, 
at least for thin aerofoils. The theoretical 
derivatives for the hinged flap exceed 
(numerically) the experimental results sub- 
stantially and this can be attributed to the 
thick boundary layer near the trailing edge. 


The first attempt to extend the foregoing 
theory to compressible flow was made by 
Possio") who gave a linearised theory for 
subsonic flow. His treatment has been 
improved in detail by Frazer’) and by 
Dietze**), Although experimental checks 
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are lacking, there is no apparent reason to 
doubt that this theory is substantially correct 
for thin aerofoils in two-dimensional flow at 
Mach numbers not exceeding say, 0.7. The 
corresponding linearised theory for super- 
sonic flow developed independently by von 
Borbely’**) and by Temple and Jahn®®’, 
has not been confirmed by the measurements 
made recently at the National Physical 
Laboratory by Bratt and his colleagues'*®). 
It is too early to pass a final judgment on 
this question but some calculations by 
Prichard Jones suggest that the thickness 
of the aerofoil, even when small, has a 
profound influence on the supersonic 
derivatives. 

Two other extensions of the theory remain 
to be mentioned. The first seeks to provide 
a reliable method of calculating derivatives 
for finite wings of arbitrary plan form and 
with any assigned modes of deformation. 
The simplest and frankly approximate 
method, which is in current use, is to 
suppose the wing to be divided into fore- 
and-aft strips and to calculate the changes 
of normal force and twisting moment on 
these by using two-dimensional aerodynamic 
coefficients. possibly with some empirical 
modifications. There are three modifications 
of this system: 


(a) Constant basic aerodynamic coefficients, 
independent of the frequency parameter 
but suitable for a mean practical value 
of it, are adopted. 

(b) The coefficients of two-dimensional 
vortex sheet theory are used, with a 
frequency parameter corresponding to a 
mean wing chord, or chord measured at 
a reference section. 

(c) Two-dimensional vortex sheet theory is 
used with the frequency parameter based 
on the local chord and therefore, in 
general, variable along the span. 


Attempts have been made to provide a 
more correctly based theory by several 
authors.  Lyon‘*’) gave a very complete 
analysis for rectangular wings but did not 
succeed in producing a workable scheme for 
calculating the derivatives, while Cicala(**) 
developed a more workable method based 
On a partial representation of the vortex 
system. More recently Prichard Jones* 4°) 
has given methods which appear to te both 
theoretically adequate and workable. The 
second extension of the theory aims at 
dealing with thick aerofoils and with angles 
of incidence which are not small, ie. it 
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aims at making due allowance for a thick 
boundary layer and wake. A beginning has 
been made by Schwarz’) and by Prichard 
Jones*), but much further research will be 
needed to provide a satisfactory theory and 
method of calculation. 


FLUTTER PREDICTION 


Critical speeds for flutter and the nature of 
the motion in the flutter are calculated on 
the basis of the dynamical equations, which 
are always in the form yielded by the 
method of Lagrange, with one equation 
corresponding to each degree of freedom. 
In order to construct these equations it is 
necessary to do three things:— 


(a) Select the degrees of freedom and assign 
the mode of motion for each freedom. 

(b) Calculate (or measure if a prototype or 
dynamically similar model exists) the 
mechanical constants, i.e. the inertial 
and stiffness coefficients. 

(c) Calculate the aerodynamic coefficients. 

We shall briefly consider these items in turn. 


For routine problems the choice of free- 
doms will have been reduced to a system. 
For example, in calculating the flexural- 
torsional flutter characteristics of an 
unswept cantilever wing not carrying engines, 
experience has shown that it is adequate to 
use just two degrees of freedom, one for 
fiexure and one for torsion. The modes of 
displacement in flexure and in torsion may 
be taken to agree with the corresponding 
fundamental free modes, if known; other- 
wise simple displacement functions of 
appropriate types may be assumed. It has 
been shown that, when the reference section 
where the flexure and twist are measured is 
at or near 0.7 of the overhang from the wing 
root, the calculated critical flutter speed is 
not sensitive to the choice of modes of 
displacement**} 4), On the other hand, 
when we are faced with novel designs the 
choice of freedoms may not be easy and 
there is always the need to minimise their 
number, since the labour of making the 
calculations rises rapidly with this number. 
This question has already been discussed 
under “ The Kinematics of Flutter.” 


The inertial coefficients of the system are 
calculated by the Lagrangian method from 
the now known displacement functions and 
the known mass distribution, with suitable 
additions for the virtual or aerodynamic 
inertias which are relatively important for 
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light structures and high air densities. If 
the resonance frequencies of the system are 
known it may be possible to calculate the 
stiffnesses from the now known inertias; 
otherwise the stiffnesses must be calculated or 
measured in static loading tests. The modes 
of deflexion in static tests should agree with 
the assumed modes and it should be noted 
that stiffnesses obtained from static tests 
with loads applied at a single locality are 
always too low. Normal modes have the 
advantage that the inertial and elastic 
couplings between such modes are absent. 
If the virtual inertias were not included when 
the normal modes were determined they will 
introduce, in general, couplings, but this 
complication does not arise when the normal 
modes are experimentally determined in still 
air. 

The aerodynamic coefficients must be 
calculated by one of the methods mentioned 
under “The Aerodynamic Aspect of 
Flutter.” 

Details of the methods of calculating the 
critical speeds are beyond the scope of this 
lecture, but may be found in the papers 
listed at the end‘?* 4% 1°), 

There is another method of estimating 
critical flutter speeds which, although not 
pretending to accurate prediction, is of high 
value. This is the statistical method, 
initiated by Roxbee Cox, and now-a-days 
made the basis of airworthiness rules on 
torsional wing stiffness. In the early days 
of the investigation of flutter, before mass- 
balance of ailerons had been adopted in 
practice, Dr. Frazer and I did not favour 
attempts at numerical predictions of the 
critical flutter speeds of wings. Our view 
was based on the fact that, with unbalanced 
ailerons, the wing flutter would be of the 
ternary type, so it would be necessary to 
know the values of upwards of 30 coeffi- 
cients in order to make a prediction. Since 
many of these coefficients cannot be 
estimated at all closely from drawings of the 
aircraft, we concluded that calculated critical 
flutter speeds would be very unreliable. This 
led Roxbee Cox**) to approach the problem 
of prediction, or at least of prevention within 
the range of flight speeds, on the basis of 
experience. He framed a suitable criterion 
for stiffness and assigned the limiting 
numerical value of the criterion for 
acceptance on the data for aircraft which 
had fluttered and for aircraft which had not 
fluttered. This approach to the question has 
been of great utility and has contributed 
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much to the achieved freedom from flutter 
of British aircraft. 

In the course of time the Wing stiffness 
criterion has been elaborated and modified 
The latest form, due to Collar, Broadbent 
and Puttick‘*’), covers the influence of the 
following factors:— 

Taper of the wing 

Position of inertia axis on the chord 

Compressibility of the air (conventional 
allowance). 


FLUTTER PREVENTION 

There are two basic methods of preventing 
flutter:— 

(a) Isolation of the motions in the several 
freedoms. (Not applicable to stalling 
flutter). 

(b) Raising the lowest critical flutter speed 
safely beyond possible speeds of flight 
by increase of the relevant natural 
frequencies. This is usually effected by 
designing for increased stiffnesses with 
a less than proportional increase of 
mass or, preferably, without increase of 
mass. 

Flutter speeds can also be raised, or flutter 
prevented, by a sufficient increase of the 
internal damping of the system, but this 
cannot be regarded as a generally applicable 
method. Further remarks on it are given 
later in relation to control surface flutter. 

As a matter of experience, most flutters 
occurring at comparatively low speeds of 
flight have essentially involved the oscilla- 
tion of at least one control surface. Hence 
the isolation (or uncoupling) of the control 
movements is the first and most important 
of the measures for flutter prevention and 
here the ailerons must receive first attention, 
since wing flutter is even more disastrous 
than tail flutter. By far the most important 
and effective measure of uncoupling is mass 
balance of the surface about the hinge, ie. 
elimination of the inertial coupling. For 4 
control surface with a reasonable degree of 
aerodynamic balance the aerodynamic 
couplings will be fairly small, although 
usually not strictly absent. 

When a control surface has been properly 
mass-balanced there will be no tendancy for 
it to move relatively to the wing (or tail) at 
the place where its operating lever 3 
connected, when the aircraft is subjected to 
resonance tests on the ground; this is the 
best proof of the efficiency of the mass 
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preferably be spread along the span of the 
control, for an isolated balance mass may 
only be fully effective for one of the possible 
modes of oscillation of the system and will 
not be fully effective if the torsional stiffness 
of the control surface is low. Remote or 
geared balance masses should only be used 
when the design has survived the most 
searching scrutiny and on no account must 
a nodal line lie between the masses and the 


control. 


Mass-balancing of controls was keenly 
disliked by designers when it was first intro- 
duced and that dislike persists on account 
of the weight of the balance masses. This 
has led to the investigation of dampers as 
devices for suppressing the flutter of control 
surfaces, but the results of the enquiry have 
been disappointing. It is not difficult to 
suppress flutter of a control by a damper 
anchored to the main structure and con- 
nected to the control, but then the added 
resistance to movement of the control by the 
pilot is intolerably great. A damper not 
anchored to the main structure but having 
a massive plunger controlled by a spring and 
tuned to the flutter frequency, can sometimes 
be effective and adds nothing to the pilot’s 
operating force. However, tuned dampers 
are unreliable because any one damper is 
not effective for all air densities and there 
would be great difficulties in the efficient 
tuning of the dampers and in making the 
damping independent of temperature. 
Moreover, it certainly is not probable that 
the damper would be lighter than the 
balance mass it is intended to replace, so 
there is really no case for it. 


The only effective alternative to mass- 
balancing of controls is to make them 
irreversible, as first suggested by South- 
well'*®) in 1927. Simple as this solution of 
the problem appears at first sight, it has not 
yet been applied in practice, except for 
power-operated controls, although a number 
of attempts to produce irreversible pilot- 
operated controls have been made. 


Once the controls have been effectively 
eliminated from the flutter, the problem of 
flutter prevention becomes one of providing 
structural stiffnesses which are adequate to 
keep the lowest critical flutter speed safely 
outside the range of possible flight speeds 
in all circumstances. It may be possible by 
clever design to reduce materially the 
magnitudes of the required stiffnesses. 


THE EFFECT OF FRICTION 
FLUTTER 


Friction at the hinges and in the operating 
circuit of a control may have an influence 
in suppressing flutter, provided that the 
disturbances are always very small. Any 
masking of flutter by friction, however, 
would be extremely dangerous as flutter 
could be initiated by a large disturbance at 
a flight speed exceeding the true critical 
speed. The flutter would therefore develop 
with an almost explosive violence. In some 
cases friction may tend to promote flutter, 
as has been demonstrated by experiment"®?. 
The theory of flutter involving “solid” 
frictional resistance is very complicated *®). 


ON 


RESONANCE TESTING THE 
GROUND AND IN FLIGHT 


Ground resonance tests can provide most 
valuable information on the modes and 
frequencies of oscillation of the aircraft 
structure which can be used in flutter 
prediction. In addition to the use of such 
data in calculations of critical speeds, the 
situation of the nodal lines may have 
diagnostic value. Kiissner long ago 
recognised that a torsional wing oscillation 
with the node at about the three-quarter- 
chord position near the wing tip was 
“dangerous” and he used the resonance 
frequency for such an oscillation in the 
flutter speed formula already quoted. This 
aspect of flutter prediction technique is now 
being pursued at the Royal Aircraft 
Establishment. 


It is important that the aircraft should not 
be unduly constrained by its supports when 
resonance tests are made on the ground. 
Ideally, the supports should have neither 
stiffness nor mass in their moving parts, but 
this is unattainable. I favour support by air 
bags inflated at a low pressure and the aim 
should be to make the pressures in the bags 
as nearly independent of the displacements 
as possible. 


Resonance tests made in flight would give 
particularly valuable information, but are 
clearly not easy to do and possibly are 
dangerous. Von Schlippe®°) proposed to 


make resonance tests at gradually rising 
speeds of flight with the object of detecting 
the approach to a critical flutter speed’. 
This technique is too dangerous for general 
use, 


but cautiously conducted flight 
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resonance tests made at speeds of flight up 
to some stated maximum would certainly 
be a most conclusive proof test. 


THE INFLUENCE OF AIR DENSITY 
ON FLUTTER 


For wings of low density, such as were 
common in the biplane era, the effect of air 
density on the critical flutter speed depends 
on the particulars of the design and the 
manner in which the critical speed depends 
on altitude requires detailed examination'*?. 
With the heavier wings which are now 
usual the lowest true critical air speed occurs 
with the highest air density and, apart from 
compressibility effects on the derivatives, it 
usually suffices to examine the flight condi- 
tion for which the ‘“ equivalent air speed ” 
is greatest. 


DUNCAN 


APPLICATIONS OF MECHANICAL 
ADMITTANCES IN FLUTTER 
THEORY 


The concept of a mechanical admittance 
or dynamic flexibility is usually restricted to 
simple-harmonic motions in undamped or 
conservative systems, but it can be extended 
to general systems whose equations of 
motion are linear and to any motions of the 
“exponential type,” for which the displace. 
ments are proportional to the real part of 
expAr with A real, pure imaginary or 
complex'’*). In general each admittance 
then becomes a complex number. I antici- 
pate that such generalised mechanical 
admittances will find useful applications in 
flutter theory, especially in dealing with the 
influence of general freedoms on the flutter. 
Valuable bibliographies of flutter and allied 
matters will be found in References 69, 70, 
71 and 72. 


PART II 
SOME TOPICS OF AIRCRAFT STABILITY AND CONTROL 


TOLERABLE AND 
INSTABILITIES 


The mathematical concept of stability is 
clear and definite but we must not, on that 
account, make a fetish of stability, for the 
things that really matter are the safety, 
comfort and handling qualities of the air- 
craft and these are not all reducible to 
mathematics, unless possibly statistical 
mathematics. Here we have a contrast with 
flutter where the attainment of stability 
(freedom from flutter) is the necessary and 
sufficient objective. Clearly an automatically 
controlled aircraft is unsafe unless the 
complete system aircraft cum auto-pilot is 
completely stable. Hence complete stability 
is here a necessary condition for air- 
worthiness, but it is not in itself a sufficient 
criterion for satisfactoriness as I shall point 
out in detail later. When the aircraft is 
controlled by a human pilot it need not be 
completely stable and it may easily be over- 
stable. 

It is obvious that the satisfactoriness or 
otherwise of an automatically controlled 
aircraft is a quality of the complete system 
consisting of auto-pilot and aircraft and this 
depends on the matching of the qualities of 
the two components of the system. Just so, 
the satisfactoriness of an aircraft with a 
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human pilot depends on the correct matching 
of the qualities of the components, namely 
aircraft and pilot. We must therefore 
enquire what are the most important relevant 
qualities of the human pilot. I suggest that 
they are:— 
(a) Muscular strength 
(b) Proneness to fatigue 
(c) Reaction time 
(d) Sensitiveness to acceleration and to the 
force reactions on the controls 

(e) Adaptability 

Items (c) and (d) are of particular 
importance in relation to the stability of the 
aircraft. Experience shows that an instability 
may be quite tolerable provided that it is 
sufficiently mild and we can measure mild- 
ness by the time required for a disturbance 
to be doubled in magnitude. When there 1s 
an unstable constituent of the motion of an 
aircraft, the disturbance itself or its ampli- 
tude when oscillatory will be proportional 
to an exponential function of st, where ¢ 1s 
the time and , a constant characteristic of 
the aircraft in its particular state of flight. 
The time 7 required for the disturbance to 
become doubled is given by 
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Now the most significant quantity as regards 
the tolerability of an instability is the non- 
dimensional quantity 
Pilot’s reaction time 
| suggest very tentatively that an instability 
will not be troublesome to an experienced 
pilot when this ratio is 50 or greater. “The 
ime T, for similar aircraft in similar circum- 


sances, will be proportional to the 
gerodynamic unit of time 

a W 

= geVS 


where W=weight of aircraft 

g=acceleration of gravity 

p=air density 

V=true speed of flight 

S=wing area. 
Since W depends on the cube of the linear 
dimensions the aerodynamic unit of time 
tends to increase with the size of the aircraft. 

Reaction time depends somewhat on the 
nature of the stimulus, on the individual and 
on his physical and mental condition. The 
reaction time of an alert and healthy person 
to a tactual or auditory stimulus is about 
0.12 sec. and about 0.17 sec. to a visual 
stimulus. Fatigue and lack of oxygen 
increase reaction time. 

There are two common kinds of aircraft 
instability which usually would be judged 
tolerable according to the above criterion. 
These are the phugoid oscillation (longi- 
tudinal-symmetric motion) and the spiral 
divergence (lateral-antisymmetric motion). 
these instabilities, unless 
unusually severe, may be ignored in aircraft 
to be operated by experienced pilots, but I 
would emphasise again that they would be 
altogether intolerable if they occurred with 
an automatic pilot. 


It is worth pointing out that the motion 


. which would follow the release of an 
untrimmed control by the pilot may, in its 


initial Stages, be violent even in a completely 
stable aircraft. 


SYNTHETIC STABILITY 


An aircraft with a suitably matched auto- 
matic pilot will be quite stable even when 
the aircraft, with fixed controls, is unstable. 
Thus the auto-pilot applies such correcting 
displacements to the control surfaces that the 
inherent instability of the uncontrolled 
aircraft Is Overcome. This same method of 
stabilisation can be used even when the air- 


craft is under the control of a human pilot. 
It is only necessary to gear the control 
surfaces so that their displacements are 
compounded of parts contributed by the 
human and automatic pilots Thus the 
human pilot controls the aircraft in the 
ordinary way, but when he holds his controls 
steady the automatic pilot takes charge and 
maintains stability. 

It is probable that synthetic stability will 
find useful applications for flight in the sonic 
and supersonic regions. 


STABILITY AND STABILISATION 


A civil aircraft with automatic control, in 
the interest of the passengers’ comfort, 
should respond as little as possible to atmos- 
pheric disturbances and the same is true of 
a military aircraft which is to serve as 4 
platform for bombing or gunnery. The 
achievement of the greatest steadiness has 
been called stabilisation by Gates and it goes 
beyond the mere attainment of stability. 
There is, in fact, a clear and important 
distinction between optimum stabilisation 
and maximum stability, as has _ been 
emphasised by Sudworth and others. 


A system is said to be stable when the 
motion of deviation or perturbation caused 
by an applied disturbance always becomes 
zero ultimately. Thus stability is concerned 
with the ultimate consequences of a disturb- 
ance and we should properly say that the 
stability is a maximum when the deviation 
ultimately dies out most rapidly*. It might 
well happen that a disturbance becomes 
largely amplified in the early stages of the 
motion when the stability is a maximum 
and that this amplification can be reduced 
or prevented when the condition of 
maximum stability is discarded. 

Accordingly the problem of optimum 
stabilisation consists in so designing and 
gearing the auto-pilot that the greatest 
excursions which follow disturbances shall 
be minimised. This problem is not easy to 
formulate in a precise way, since there are in 
practice many different kinds of disturbance 
and many different components in the 
response. In any given case it would be 
necessary to fix on those components of the 
response which are particularly objection- 


*This means that the numerically smallest real 
negative root or real negative part of a complex 
root of the characteristic equation is as numeric- 
ally great as possible. For stability no real root 
or real part of a complex root can be positive. 
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able and to minimise these for the most 
frequent or important kinds of disturbance. 


RESPONSE CALCULATIONS 


Response calculations are concerned with 
predicting the motion of the aircraft induced 
by operating its controls in a_ prescribed 
manner, or caused by gusts varying in time 
and space in some known way. The 
problem consists in solving the relevant 
dynamical equations with initial conditions 
corresponding to the undisturbed state of the 
aircraft and with the applied forces and 
moments corresponding to those produced 
by the controls or gusts. Two distinct cases 
arise accordingly as the deviations from the 
initial state always remain small or become 
large. Whenever the deviations become too 
large for their squares and products to be 
neglected, the dynamical equations become 
non-linear and, in general, are not soluble 
analytically. Then response calculations 
must be done by the step-by-step numerical 
integration of the equations of motion or 
with the help of a suitably arranged differen- 
tial analyser. But when the deviations are 
small enough for linearised dynamical 
equations to be used there is a number of 
applicable analytical methods, while the 
differential analyser is again available. 

A considerable variety of methods exist 
for the analytical solution of the linear 
differential equations of response theory, for 
example, Heaviside’s method or the allied 
method of the Laplace transform and 
methods making use of matrices**). I wish 
to draw attention to a method which appeals 
to me particularly, on account of the 
simplicity of its physical basis. This is the 
method of impulsive admittances*) and to 
explain it I shall take one definite simple 
case. Let it be required to find the 
manner in which the pitch angle varies with 


time when the elevator is moved in a 
prescribed manner, the aircraft being 
initially trimmed and flying steadily. Let 


us suppose that we know the pitching 
moment corresponding to each elevator 
setting. Then from the prescribed elevator 
movement we determine the way in which 
the pitching moment due to the elevator 
varies with time*, say M(t). Now suppose 


*Calculated on the assumption that the speed of 
flight is constant. The speed of flight really varies 
but the effect on the pitching moment 7M, due 
to the elevator is a small quantity of the second 
order. 
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that we find by solving the equations of 
motion that an instantaneous impulsive 
pitching moment of unit impulse applied 
when t=0 causes the angle of pitch to 
deviate from the undisturbed value by the 
amount A(t) at time t. Then the deviation 
of the angle of pitch at time ¢ under the 
_ of the pitching moment M(j) is given 
y 


6(t) A(t- T)M(T)dT. 


In order to see this we note that the 
pitching moment at time T is M(T) and that 
its impulse during a short interval dT js 
M(T)dT. If this were the only impulse 
acting on the aircraft the angle of pitch at 
a later time ¢ would be A(t-T)M(TMdT. 
Now the total effect of the applied pitching 
moment throughout the time interval 0 to 
t is obtained by adding the effects of all the 
small impulses which may be taken to 
replace the applied moment and this is just 
the integral already written. 

The function A(t), giving a_ particular 
component of the response to a particular 
kind of unit impulse, is called an impulsive 
admittance and we see that it is one com- 
ponent of the solution of the dynamical 
equations of free motion following on the 
administration of a particular kind of unit 
impulse to the undisturbed aircraft at the 
zero of time. If, for the type of motion con- 
sidered, the number of degrees of freedom 
is n then there will be n? independent basic 
impulsive admittances in all. We can extend 
this idea and say that if some quantity Q 
in which we happen to be interested 1s 
linearly expressible in terms of the deviations 
such as (tf) and the like, then it can be 


expressed by means of an integral of the | 
same type as that already given, or by a} 


sum of such integrals if there are several 
forces or moments acting simultaneously. 
For example, we could in this manner 
derive an expression for the normal accelera- 
tion at the C.G. caused by an arbitrary 
movement of the elevator. 


We have remarked that each impulsive 
admittance gives an expression for @ 
particular component of displacement, 
velocity or acceleration in the free motion 
which follows the administration of 4 
particular kind of unit impulse. Since we 
are dealing with a free motion we caf 
express the impulsive admittance as a linear 
combination of the displacements and the 
like, in the basic free modes of the system 
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and Mitchell'*”) has called the coefficients 
of these modes in the admittance modal 
response coefficients. 1 should add that 
Mitchell’s approach to the matter is very 
diferent from mine and that he does not 
ye the terminology of admittances; he 
arrives at his results by the method of 
“variation of parameters.” However, the 
connection between the methods is as stated. 


THE INFLUENCE OF DISTORTION 


It is usual to suppose that the influence 
of distortion upon aircraft stability is a very 
“modern” subject and it is therefore of 
interest to recall that the Wright Brothers 
were not only aware of this influence, but 
deliberately used it to improve the stability 
of their aeroplanes. The subject has only 
fairly recently become prominent again on 
account of the greatly increased speeds of 
modern aircraft. So far as | am aware, the 
discovery of the sensitiveness of the 
characteristics of ordinary flap controls, and 
especially of their hinge moments, to distor- 
tion of their surfaces is recent, although 
it has long been known that the hinge 
moments of Frise ailerons are very sensitive 
to changes of the nose shape. 

The thoroughgoing- and potentially com- 
plete method of investigating the effects of 
distortion is to introduce new dynamical 
co-ordinates to represent the distortional 
displacements'**: °*). This implies that there 
is an additional dynamical equation for each 
new co-ordinate and consequently the labour 
of forming and solving the dynamical 
equations is much increased. There is no 
substitute for this method when the rates of 
distortion in the motion considered are so 
large that the corresponding inertia forces 
of the structural elements are important. 

In many problems of aircraft dynamics 
the structural distortions occur relatively 
slowly and the corresponding inertia forces 
are negligible. For an oscillatory motion of 
frequency f this will certainly be true when 
the ratio f/f, is less than 0.1, where 
fs the fundamental structural frequency **). 
Then, for any one type of general motion of 
the aircraft, it is legitimate to assume that 
the distortions are linear functions of the 
deviations of the general motion from the 
equilibrium state. Consequently, the 
changes in the aerodynamic forces and 
Moments caused by the distortions are 
linearly related to the deviations and may 
be absorbed into the expressions for the 
aerodynamic derivatives. Thus, in this 


STABILITY 


method, which has been developed by Gates 
and Lyon), distortions are allowed for by 
modifying the values of the aerodynamic 
derivatives. These modifications depend 
importantly on the speed of flight. 


COMPRESSIBILITY 


The compressibility of the air for flight 
at higher Mach numbers gives rise to 
important changes both of trim and stability. 
Trim changes are associated with changes in 
the distribution of aerodynamic force over 
the surface of the aircraft and these are 
attributable to changes in the general regime 
of flow and to wake and downwash effects 
at the tail. The influence of compressibility 
on stability is attributable to changes in the 
values of the aerodynamic derivatives and 
there are several causes for such changes:— 
(a) The lift, moment and drag coefficients 

and their rates of change with incidence 
depend on Mach number. 


(b) The fact that these non-dimensional 
coefficients are not independent of the 
speed of flight for a fixed incidence 
introduces additional terms into the 
expressions for the derivatives. 

(c) Wake and downwash effects at the tail 
are not independent of Mach number. 


We should add that scale effects, i.e. depend- 
ence of the non-dimensional force and 
moment coefficients on Reynolds’ number, 
are mixed with the compressibility effects. 


Several of the troubles attributable to 
compressibility (and to distortion) can be 
overcome by using suitably designed power- 
operated controls. The power-operated all- 
moving tailplane without elevators is to be 
preferred for high speed flight as it gets rid 
at one stroke of the difficulties of aero- 
dynamic balance and of tailplane and 
elevator distortion. 


USE OF QUASI-STATIC DERIVATIVES 


It is usual to employ in stability calcula- 
tions values of the aerodynamic derivatives 
obtained from wind tunnel experiments on 
models made at a number of fixed angles 
of incidence or yaw. In any such test the 
circulations at all parts of the wings and tail 
are steady and there is no net vorticity with 
transverse axis of spin in the wake. This is 
no longer true when the aircraft has a 
motion of disturbance and the transverse 
vorticity in the wake induces velocities at 
the wing and tail which influence the aero- 
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dynamic forces. Consequently the quasi- 
Static derivatives obtained from model tests 
in the manner described are not applicable 
to flight. However, detailed analysis shows 
that the error in quasi-static derivatives 
should be less than five per cent. in an 
oscillation for which the non-dimensional 
frequency parameter 
(2zfcem)/ V 

is less than 0.02. This criterion is satisfied 
by the phugoid oscillation but not, usually, 
for the rapid longitudinal oscillation’?*). 
However, the resulting errors in the deriva- 
tives for such oscillations are probably not 
very serious. (See Part I under ‘“ The 
Aerodynamic Aspect of Flutter.’’) 


THE MEASUREMENT OF AERO- 
DYNAMIC DERIVATIVES 


It is possible to estimate the values of 
some of the aerodynamic derivatives of 
stability theory by pure theory, or on the 
basis of wind tunnel tests on aerofoils and 
other simple bodies. Such estimates require 
experimental confirmation and other deriva- 
tives cannot be predicted with any certainty. 
Consequently the measurement of derivatives 
is a matter of importance. 

There is one simple general maxim 
applicable to all measurements of derivatives, 
namely, “use the dynamical equation in 
which the derivative occurs,” and we may 
note that any selected derivative only occurs 
in one of these equations. The experiment 
must then be so devised that the unknown 
derivative can be found with the greatest 
accuracy and convenience. This implies that 
with the help of the dynamical equation the 
derivative can be expressed in terms of 
experimentally measured quantities, such as 
amplitude ratios, phase differences and 
frequencies, which be accurately 
determined, and other already determined 
coefficients, such as moments of inertia and 
stiffnesses. The value of the derivative 
should not depend on that of another 
derivative of doubtful magnitude. 


Whenever possible there should be only 
one degree of freedom in the experiment. 
All “ direct ” derivatives can be found from 
such experiments and some “cross” 
derivatives (aerodynamic couplings) can be 
found from the results of single freedom tests 
employing in succession two or more rather 
widely separated axes of rotation. In 
planning experiments it must not be for- 
gotten that derivatives depend on _ the 
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“frequency parameter” and on such aero. 
dynamic factors as transition point in the 
boundary layer. Unless the experiment jg 
made full-scale (or possibly in a compressed 
air tunnel) it will not be possible to achieve 
the true Reynolds’ number for flight. 


Experiments on models can be made in 
a wind channel or on a whirling arm. Wind 
channel tests on static models can be made 
the basis of estimates of quasi-static deriva- 
tives, as already discussed. Dynamic tests 
in wind channels may be made with con- 
tinuous rotation on a rolling or spinning 
balance or in oscillatory tests®*). The latter 
include free and forced oscillation tests, 
while the forcing may be mechanical or 
electrical. In the technique developed by 
Bratt at the National Physical Laboratory 
damping derivatives are measured by 
coupling the oscillating model to an electrical 
device which can provide a controlled and 
measured amount of negative or positive 
damping’*’). The electrical damping is 
adjusted until the system is neutrally stable, 
as shown by the persistence of a simple- 
harmonic motion and then the aerodynamic 
damping is equal and opposite to the 
electrical. The complete history of a force 
or moment during a movement can be 
obtained with the aid of a suitably arranged 
Strain gauge (magneto-striction or wire 
resistance type). Whirling arm tests permit 
the measurement of forces and moments in 
continuous rotation in pitch or yaw. These 
tests are made difficult by the swirl set up 
by the whirling arm and the model itself and 
by centrifugal effects. 


Full-scale measurements of derivatives, in 
general, are difficult but are much to be 
desired, since the uncertainties about scale 
effect in model tests are avoided. Since, 
however, it is relatively difficult to make a 
full-scale aircraft effectively rigid, it 1s 
important that all relevant distortions should 
be measured. It is also of the utmost 
importance to ensure that no control surface 
moves when it should be fixed, as even 
extremely small control movements may giv¢ 
rise to large moments. Further, the 
propulsive system must be strictly controlled 
in accordance with a predetermined scheme. 
Measurements of some quasi-static deriva- 
tives in flight can be made by applying 
moments, e.g. by shift of large masses or by 
means of drogues. 


A new and interesting technique has been 
developed in America and is described by 
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Milliken’. In this method a_ simple- 
iarmonic oscillation is impressed on a 
control surface by means of a_ suitably 
nodified automatic pilot and the response of 
the aircraft is fully recorded. 


AERODYNAMIC REVERSAL OF 


CONTROL 

In a paper’’’) read before the recent 
Congress on Applied Mechanics, R. A. Shaw 
described an aerodynamic reversal of the 
control force provided by a conventional 
hinged flap, occurring at a high subsonic 
Mach number (M in the region 0.85 to 0.87 
for the aerofoil tested. This was EC. 1250 
with a 25 per cent. control flap). The 
reversal of control was accompanied by 
reversal of hinge moment. Similar effects 
had previously been observed by Gothert on 
an aerofoil 9 per cent. thick and by 
Lindsey on an aerofoil 19 per cent. thick. 
This is unquestionably a potentially trouble- 
some and dangerous phenomenon, but it is 
evidently not an inevitable one since aircraft 
have flown right through the ciritical region 
without experiencing trouble. 

It is clear that the usual cure for reversal 
of control, namely increase in torsional stiff- 
ness of the main surface, is of no avail. 
There are indications that the trouble may 
be avoided by using aerofoils whose 


PART III 
SOME URGENT QUESTIONS 


IS THERE A “FLUTTER BARRIER”? 


Since nothing definite is known about the 
aerodynamic derivatives of flutter theory for 
fight in the neighbourhood of the speed of 
sound, it is reasonable to enquire whether 
there is an “iron curtain” of flutter to 
prevent passage through the sonic region. 
Theory can give little guidance here but the 
equation of energy would point to the con- 
clusion that flutter could be prevented, even 
in the sonic region, by the well-tried method 
of stiffening, with the proviso that no direct 


| aerodynamic damping derivative takes the 


unstable sign. Such experimental evidence 
as has been obtained for high subsonic and 
low supersonic speeds may be taken to 
indicate that reversal of sign in the dampings 
18 less likely to occur with higher values of 
the non-dimensional frequency parameter 
than with low, so here again the indication is 
in favour of high structural stiffness. 


-aerodynamic reversal of control already 


STABILITY 


maximum thicknesses are not too far back 
and by making the trailing edge angle of the 
control flap small. It appears that reversal 
is more likely to occur with flaps of large 
chord than of small and it is possible that 
there is a connection between this phenome- 
non and “aileron buzz.” 


AILERON BUZZ 


Aileron buzz is a rapid oscillation of the 
aileron occurring at high subsonic Mach 
numbers. It has been experienced on some 
American aircraft and the experiments of 
Smilg’°*) appear to demonstrate that it is 
essentially an oscillation in a single degree 
of freedom (flap angle). Smilg concludes 
that the hinge moment is out of phase with 
the displacement in such a sense that 
instability results and this must be true if 
no other degrees of freedom are concerned. 
It may be that aileron buzz is related to the 


mentioned. 

There is no evidence so far of the occur- 
rence of aileron buzz on British aircraft, but 
the cause of this immunity is not known. 
British requirements for the stiffness of 
aileron control circuits are more stringent 
than the American, but it seems improbable 
that this in itself would account for the 
immunity. 


Fortunately we can already say categoric- 
ally, as the result of experience, that flutter 
in the sonic region is not inevitable, for 
certainly both model and full-scale aircraft 
have successfully flown at speeds in excess 
of the sonic. Because of these successes, 
however, we must not assume lightheartedly 
that there are no particular dangers of flutter 
in the sonic region. Certainly no full-scale 
aircraft should be permitted to attempt flight 
near the speed of sound unless its stiffnesses 
are quite satisfactory. 


FLUTTER OF AIRCRAFT WITH 
HEAVILY SWEPT WINGS 


Knowledge of the flutter of flying wings 
and other aircraft having heavily swept wings 
is scanty. It is already clear that for such 
aircraft the flutter and general stability 
characteristics are so entangled that the 
distinction ceases to be useful. Thus, the 
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general freedoms of the aircraft, such as 
normal movement of the C.G. and pitching, 
play an essential part in the flutter. For 
flying wings the stabilising influence of free- 
dom in roll on flutter of the antisymmetric 
type is so great that it is very unlikely to 
occur. 

Experiments on the flutter of model flying 
wings made at the N.P.L. by Lambourne“®*) 
show that the general freedoms may reduce 
the critical speed for symmetric flutter. 
Flutter can occur in which the important 
freedoms are normal translation of the C.G., 
pitch and wing bending. It is remarkable 
that wing twist plays no essential part in this 
flutter and it is effectively replaced by the 
general pitching motion. It follows that the 
usual method of preventing flutter by raising 
the torsional wing stiffness is here altogether 
ineffective and that increase of the flexural 
stiffness is the curative measure. However, 
we must not assume that torsional wing 
stiffness is of no account, for it is required 
for the prevention of reversal of control and 
of wing divergence. Further, if the torsional 
wing stiffness were very low the flutter would 
no doubt have a torsional component and the 
critical flutter speed would be affected. 


BUFFETING IN HIGH SPEED 
FLIGHT 


Ordinary buffeting is a more or less heavy 
vibration of the tail caused by an eddying 
wake impinging upon it and transmitted 
from the tail to the whole aircraft'**). 
Usually the eddying wake springs from a 
wing root or wing-nacelle junction and is 
only serious for incidences near the stall. 
The buffeting is prevented in a fundamental 
way by good aerodynamic design of the wing 
junctions. It can be mitigated by raising the 
tail out of the wake and by ensuring that no 
natural frequency of the structure agrees with 
that of an important frequency in the wake 
“spectrum.” Leakage of “dead” air into 
a wing junction, or into the region between a 
fuselage and nacelle, may cause severe 
buffeting. 

Recently buffeting has become a serious 
nuisance in flight at high Mach numbers and 
has prevented the use of the higher lift 
coefficients which the wing could provide. 
Both the cause and the means for preventing 
this kind of buffeting are largely unknown; 
it is not even clear whether the buffets caus- 
ing the vibration are administered to the tail, 
or to the wing itself, or to both. Knowledge 
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of the influence of tail position and of the 
structural frequencies is lacking. My guess 
is that this species of buffeting will Only be 
avoided when there is clean flow all ove; 
the aircraft, for it is obvious that a very small 
fractional disturbance is important when the 
speed of flight is very high. Clean flow mus 
be sought in the usual ways but perhaps it 
may be unattainable at higher incidences 
without the aid of suction through slots or 
porous regions. 


POWER-OPERATED CONTROLS 


It is not my object to discuss this matter 
fully and such a discussion is rendered 
unnecessary by the recent able lecture by 
Lyons’*”), I should like to say a few 
words about flutter-prevention with powered 
controls. 

Power operation implies additional weight 
and it is therefore important to dispense with 
the balance masses on control surfaces if this 
can be done safely. There is here just one 
requirement:— 

When the pilot’s contro! is fixed the 

control surface shall either be truly locked 

for all possible loads on the surface, or so 

heavily damped that flutter cannot occur. 
The phrase “truly locked” implies that, 
when the pilot’s control is fixed, there is no 
measurable free movement of the control 
and that the elastic stiffness against such 
movement is so great that the natural 
frequency of the control would be much 
above that of any possible flutter. This 
requires great stiffness in all the following 
items:— 

The contro] surface itself and the supports 

of its hinges. 

The power unit and its supports. 

The connection between the power unit 

and the control. 
Thus it appears to me that success is largely 
a matter of good structural design and good 
“engineering.” In my opinion no aircraft 
should be allowed to fly without balance 
masses until the controls of the prototype 
have been thoroughly tested on the ground. 
The tests should include forced oscillation 
tests over a wide band of frequencies and 
with unusually large force amplitude, as well 
as measurements of backlash and stiffness 


GUST ALLEVIATION 

It is desirable to alleviate the effects of 
gusts in order to reduce the stresses in the 
structure and to improve the comfort of 
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passengers —perhaps we should call this 
“disgust alleviation.” Clearly two elements 


are required:— 


(a) A gust detector 


(b) An alleviator governed by the detector. 


both purposes. 


It is possible for one contrivance to serve 


Alleviation can be obtained by a change 
of incidence or of flap angle of magnitude 
and timing suitably related to those of the 
gust, and timing is the crux of the matter. 
However, the alleviation probably cannot 
be complete without spoiling the stability 


and control 


alleviators must 


of the aircraft. 
be carefully scrutinised 


All gust 


from the aspect of flutter prevention. 


APPENDIX 
THE EQUATION 


OF ENERGY FOR 


FLEXURAL-TORSIONAL WING 


FLUTTER 


For the flutter which occurs exactly at the 
critical flutter speed, and which is a simple- 
harmonic motion, the intake of energy per 
cycle is exactly balanced by the output of 
energy (dissipation by damping actions). 
Detailed analysis“) shows that the condition 
of energy balance can be expressed by the 


equation 


(K, - sin 
=o B,+(, +B)(®) cos Gy | 


where the ‘symbols have the following 


significance : — 


§=amplitude of torsional motion 
g=amplitude of flexural motion 
8=angle by which torsional oscillation lags 


behind flexure 


»=2z times frequency of the oscillation 


Aerodynamic stiffness deriva- 


—M,) tives, proportional to V?. 


K,=-Ls 
B,= - L, 
- 
- 

-M. 


Aerodynamic damping deriva- 
; | tives, proportional to V. 


L= flexural moment 
M=torsional moment 


Note—K, is large and positive while C, is 
zero or very small. 
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B, and J, are the direct flexural and 
torsional dampings, respectively, and 
both are normally positive. 

J, and B, are cross-damping coeffi- 
cients. 


The expression on the left hand side of 
the equation is proportional to the energy 
intake per complete cycle and the expression 
on the right is proportional to the energy 
dissipated per cycle. We note in particular 
that the energy input is proportional to 


(a) The sine of the angle of lag P. 
(b) The torsional amplitude. 


On the other hand the energy dissipated 
per cycle is proportional to », i.e. to the 
frequency. Hence, in general, an increase in 
frequency has a stabilising influence. It will 
be seen that, if w, 8 and @/¢ were all constant, 
the energy input would be proportional to 
the square of the wind speed and the 
dissipation proportional to the wind speed. 
Hence, as V rises, there is a tendency for 
input to overtake dissipation, but equality 
may not be reached within the range of flight 
speeds or indeed at any speed (see the curves 
for the mass-balanced wing in Fig. 5). 


A result of interest follows on substituting 
the appropriate non-dimensional aero- 
dynamic coefficients. We then have 


KE,=e 
€.=p exc. 
B,=p VS*cemb, 


mj, 
mb, 
J,=p Vse*mjs 
and the equation of energy can accordingly 
be reduced to 
(£28) cos B+j (st 


(k, sin 8 


where s=semi-span 
and c,,=mean chord. 


The expression on the right hand side of the 
last equation is non-dimensional and can be 
identified with K /2z7 (see the text under “The 
Mechanical Aspect of Flutter”). 


The stabilising influence of an increase in 
frequency follows from a simple general 
argument. Let the amplitudes of oscillation 
and the phase relations be supposed fixed. 
Then the velocities are proportional to the 
frequency and so also are the damping forces, 
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but the lift force in phase with the twist is 


independent of the frequency (if we neglect 
the minor dependence of the non-dimensional 
aerodynamic derivative coefficients on the 


frequency parameter). 


Thus the forces 


responsible for the influx of energy are 
independent of frequency, whereas the forces 
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DISCUSSION 


Dr. H. Roxbee Cox (President, Fellow): 
Flutter theory and flutter experience had 
been developing now for a long time. When 
Professor Duncan started, wings were more 
or less rectangular and they had more or 
less straight spars so that they could have 
a set of flexural centres which made a 
reasonably straight line in the wing. 
Structures had become more and more com- 
plicated so that the loci which they traced 
out on wings had got more and more curved 
and the application of theory had become 
more difficult. Not only that; they had 
gone faster and faster so that they had 
introduced into an already sufficiently 
difficult problem compressibility — effects, 
which meant that their knowledge of aero- 
dynamic derivatives was once more almost 
back to zero. They had introduced swept- 
back wings which made things still more 
complicated and in addition there was now 
a new field—turbines—in which flutter had 
to be studied. They were beginning to work 
on the flutter of blades of the gas turbine, 
which would bring its own complexities. 
Fortunately the degrees of freedom in that 
case seemed to be limited and it looked 
as if there would be not more than three and 
possibly only two, degrees of freedom. 
Furthermore, the possibility of flutter 
problems with the blades of axial com- 
pressors in which there was a complex aero- 
dynamic condition with a plurality of rows 
of blades had to be considered. 

There would be developments of the 
theory for a long time to come, so that it 
was most important that they should have 
a clear understanding if they were going to 
study the various new aspects of coupled 
oscillations. 

He had always had some difficulty in 
dealing with the “ Centre of Independence ” 
—which in a well-behaved wing could be 
translated into a nice straight Axis of 
Independence. He was not clear whether 
or not that axis bore a definite relation to 
the aerodynamic axis (“quarter chord 
axis”). He was never sure whether they 
were one and the same or just reasonably 
near together! He hoped that Professor 
Duncan would clear the matter up for him. 

Professor A. R. Collar (University of 
Bristol, Fellow): At the end of 1929 the 
Superintendent of the Aerodynamics Depart- 
ment of the National Physical Laboratory 
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instructed the most junior member of th 
staff to go and see Mr. Duncan (as he then 
was), because Mr. Duncan was in th 
middle of a flutter problem and needed a 
office boy and general helper. He was 
bound to say that he (Professor Collar) ha 
never regretted it. 

Since then a great deal of air had flowe 
through the wind tunnels; neither the Super. 
intendent, nor Professor Duncan, nor he 
himself were still members of the staff of 
the Nationai Physical Laboratory; and 
where they used to spend a great deal of 
time and labour in measuring the aero. 
dynamic derivatives corresponding to a 
binary flutter problem, and occasionally 
even a ternary case; or even looking further 
afield and wondering academically whether 
anyone would ever go on to a quarternary 
case, or higher. Nowadays flutter research 
workers were quite happily taking in their 
stride six, seven, eight or even 12 degrees 
of freedom in flutter problems. He thought 


they could still all learn from such exposi- | 


tions as had been given to them by the 
lecturer. 

Professor Duncan had commented on the 
attitude of designers to mass balancing and 
their dislike of it when first introduced 
because of the weight which it entailed. 
That dislike was quite understandable; and 
it was understandable, too, that sinc 
nowadays it involved considerable mass, it 
was still to some extent disliked. But he 
thought it should be looked at mor 
rationally; the balance mass fitted to prevent 
flutter was really an essential part of the 
structure, as much so as the material in the 
spar booms or the wing skin. One of the 
objections to it was that, in the early days 
at least, it was common for the balance 
weight to be made of lead; and no one 
associated any strength with lead. Never 
theless that piece of lead contributed to the 
strength of the aeroplane, and still did. But 
even nowadays there was a tendency to 
put on balance mass as an after-thought 
rather than to regard it as an essential part 
of the structure, when it would be possible 
to consider making the metal of the 
balance mass do some useful structural work 
in the construction of the control surfaces 
themselves. 

Professor Duncan favoured air bags for 
supporting aircraft during resonance tests, 
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FLUTTER AND STABILITY 


but he thought they had some way to go on 
the question of supporting aircraft during 
resonance tests, and he would like to hear 
other comments on that. On the last 
gcasion on which he had seen air bags 
used they had not proved very successful. 
Air bags introduced a considerable amount 
of damping and they also tended to act 
rather as roller bearings, and the aircraft 
displayed a tendency to roll off. It was 
necessary, therefore, to put some stiffening 
diaphragms into ihe air bags, which intro- 
duced undesirabie stiffnesses from the point 
of view of the experiment. 


The problem of sweepback was urgent. 
Could Professor Duncan tell them what, in 
his view, would be the most fruitful 
researches which could be made in the 
sweepback area? At the moment where 
did they most need information? 

In modern flutter investigations, research 
workers took into consideration very large 
numbers of degrees of freedom; but in order 
to do that they had to use theoretical 
expressions for the air forces. In the early 
days there were no theoretical expressions 
available and Dr. Frazer and Professor 
Duncan had measured aerodynamic forces 
and related them to the actual flutter 
experiments being made. They had seen 
some small experiments at the lecture, not 
perhaps particularly abstruse or scientific in 
their nature, but nevertheless more convinc- 
ing than a theoretical exposition would have 
been. Did Professor Duncan not feel that 
the time had come for a return to a more 
experimental approach to flutter? Did he not 
consider that the present approach was too 
theoretical, and that more time should be 
devoted to the experimental, rather than the 
theoretical, aspects? 

Professor A. C. Pugsley (University of 
Bristol, Fellow): One of the main features 
of the paper and lecture was an attempt to 
simplify the presentation of flutter for the 
average engineer. It was interesting to look 
back on previous work in that field and 
notice the different ways in which it had 
been done. Professor Duncan’s way struck 
him as particularly successful. It was a 
return in part to a way which was attempted 
by Sir Ben Lockspeiser some ten or fifteen 
years ago, when he outlined the nature of 
flutter in terms of the energy and phase 
charges involved. Professor Duncan had 
given a new and very realistic twist to those 
ideas. By his flutter engine and like 


experiments he had brought life to the 
energy conception, and by the neat form 
of his energy equation he had given the 
conception a wider power. 

He believed that this lecture was about 
the sixth of its type in the history of the 
Society and the second by Professor Duncan 
—about the history of flutter theory. 
Looking back over the 20 years since Frazer 
gave his first lecture on the subject to the 
Society, he wondered whether there had 
been as much fundamental advance as they 
would have liked. That was not to decry 
the work that had been done. In one sense 
it was, he thought, a great compliment to 
the work originally done by Frazer and 
Duncan, and to the flutter “ Bible.” It was 
hard to find many fundamental advances 
since that time. Most of the advances in 
which many of the audience had taken part 
had been primarily elaborations of the 
earlier work. He very much advocated, as 
did Professor Collar, a return to a more 
experimental approach. This should provide 
the fundamental knowledge now particularly 
needed in relation to high speed flight, and 
at the same time, might lead them through 
the forest of “degrees of freedom” now 
pressing upon them. 

J. D. North (Boulton Paul Aircraft Ltd., 
Fellow): The discussion of stability and 
control in Part II of the paper very rightly, 
in his view, drew attention to the fact that 
the complete aeroplane with the human 
operator was, in fact, a servo system which 
might be, and often was, stable even when 
the aeroplane was not stable with the 
controls fixed. He was somewhat surprised 
at certain of the suggestions which were 
given as to the criteria with regard to 
controllability. 

Reaction time was only one element and 
not a major element of the total lag of the 
circuit which in its effect on stability was 
exactly the same as any other lag, although 
because of the nature of the loop transfer 
function of an aeroplane, there was no fixed 
lag. The effect was roughly to add a con- 
siderably larger quantity on to the pilot’s 
reaction time and it was immaterial whether 
that lag were in front or behind the pilot. 
It was the time taken for the correction to 
appear at the same point as the disturbance 
which mattered. That meant that since the 
aerodynamic unit of time not only affected 
the rate of divergence, or convergence, but 
also affected the rate of response, it put a 
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DISCUSSION 


slightly different picture on what could be 
handled in the case of large aeroplanes. 

He was also a little surprised that 
Professor Duncan only referred to the rate 
of increase of a disturbance in the case of 
instability as being a criterion. There was 
good evidence that frequency played an 
important part. If the frequency were high 
—certainly if it were of the order of ten 
radians a second or higher—unless it were 
of a purely rythmic character it was almost 
impossible to handle by manual control, 
even if the logarithmic increment were so 
slow that the disturbance doubled itself 
within a half hour. 

The reaction time as Professor Bartlett 
had pointed out, was not the same in 
individual and serial tasks. In the latter 
case it was considerably longer and_ its 
magnitude depended on the difficulty of the 
task as well as the stimulus. In quite 
elementary tasks the average lag which could 
be extracted by analysis from circuits 
indicated that something in the neighbour- 
hood of .3 seconds for visual stimulus was 
quite usual, even where only one degree of 
freedom was being controlled. 

Professor Duncan had selected the force 
response and excluded the displacement 
response of a pilot. The displacement 
increment probably played at least an 
important role on the force increment, 
indeed the force itself might be of more 
value as trim reference. 

He was not quite clear what Professor 
Duncan meant by “ match” in connection 
with the automatic pilot. There were certain 
restrictions which had to be satisfied, but 
he would not have thought that the 
characteristics of the automatic pilot had to 
be matched to the aircraft within narrow 
limits. If it could be treated as a linear 
problem then it was just one of the elements 
in the loop transfer function and, provided 
that the final loop transfer function satisfied 
certain stability criteria, the result was 
satisfactory. As in flutter it was the limits 
in one direction that mattered. 

Professor Duncan had given an interesting 
suggestion for the evaluation or expression 
of loop transfer functions: did he know of 
Professor Tustin’s method of expressing 
them by serial operators? This method was 
fully described in the JOURNAL OF THE I. OF 
E.E., Vol. 94 Ila 1947. The possibility of 
dealing with response problems by operators 
was discussed in the JOURNAL OF AERO- 
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NAUTICAL SCIENCES—September  194g— 
Application of Performance Operator to 
Aircraft Automatic Control. Perhaps 
Professor Duncan was not familiar with the 
original Cornell University work which 
formed the basis of Milliken’s paper, 
reference 58, as the simpie harmonic 
oscillation disturbances were abandoned in 
favour of step function disturbances because 
of the great time saving. 


D. J. Lyons (Royal Aircraft Establish. 
ment, Assoc. Fellow): Professor Duncan had 
dealt briefly with the use of flight resonance 
testing. It seemed quite true that there were 
certain dangers in this practical method of 
testing for the onset of flutter, but he 
submitted that it was far safer to use this 
than just flying the aircraft to its limiting 
speed, albeit gradually, to see if flutter 
occurred. The need for flight resonance 
testing was most critical when a new or 
strange system was proposed. Certainly full 
flutter calculations must be done to show 
that the system was theoretically free from 
flutter before flight tests could be envisaged, 
but those responsible for the flight tests 
naturally felt a little dubious about testing 
the new system in flight, because of the 
unknown accuracy of the theory as applied 
to the new system, bearing in mind the 
scarcity of information on the aerodynamic 
derivatives important in flutter. They had 
one such system in view at Farnborough, 
where it had been shown theoretically that 
it was possible to remove all mass balance 
weights from both control and tab of certain 
servo tab systems. This was a radically new 
method and doubts as to its effectiveness 
must arise. Flight resonance testing, if 
carefully applied, gave them a good chance 
of noticing the onset of an unpredicted 
flutter by the growing amplitude of oscilla- 
tion provoked by the vibrating machinery, 
which otherwise might be encountered 
catastrophically. He urged that experience 
with the use of flight resonance testing 
should be gained in Great Britain. 

He felt rather strongly on the insufficient 
knowledge of the aerodynamic derivatives 
affecting flutter. He would like to suggest 
that some flight measurements, at reasonable 
Reynolds numbers, of the dominant 
derivatives should be undertaken. It should 
be possible to make a sort of flying wind 
tunnel, by mounting a model over the top 
of an aircraft, and to measure the required 
derivatives on this model more or less free 
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FLUTTER AND STABILITY 


of interference from both the aerodynamic 
and vibration eflects of the parent aircraft. 

During the section dealing with stability 
Professor Duncan had said in one place that 
powered controls gave them a means of 
overcoming some stability troubles arising 
from distortion, and cited the case of the 
powered all-moving tailplane as overcoming 
ihe distortion troubles associated with tail 
twist. He differed from this opinion. All 
that powered controls could do was to over- 
come control balancing problems—they 
could do no more than that. They could not 
affect the basic stick-fixed changes in stability 
—but could only affect the stick-free changes. 
Also, the all-moving tailplane did not 
obliterate tail twist on an aircraft which was 
to travel in the sonic region. For the aero- 
dynamic centre would move, say, from the 
1'to the 4 chord point as the Mach number 
was increased, and with the fixed flexural 
axis of the tailplane, twist could only be 
avoided in either subsonic or supersonic 
conditions, not both. 

C. G. A. Woodford (Consulting Engineer, 
Assoc. Fellow): He suggested that the mass 
of the motor which might be put in for the 
electrical power operation of controls should 
be used as part of the mass balance of the 
wing, and that the system should therefore 
be designed at the time that the wing design 
was begun and not at the end. Did 
Professor Duncan think that feasible? 


Not only had the aeroplane to be flown 
manually, but it had also to be flown by 
automatic pilotage, so that they must con- 
sider all the reactions from that angle and 
that, he imagined, would present a more 
exacting limiting case than possibly manual 
flight control. In such a case the problem 
was intensified by the fact that the tendency 
now with automatic pilotage was to intro- 
duce what he would call ‘“ anticipatory 
sensitivity.” That was to say, they were 
tending to try to create “rate of change” 
tesponse rather than change of direction 
response. They were working to the second 
degree of the variant and by virtue of 
introducing anticipatory sensitivity it seemed 
to him to emphasise the nature of the 
control system that they had to put in. That 
should be borne in mind as to the different 
effect to what would normally be the case 
with the manual pilot, and it gave weight, 
in his view, to the need to consider the 
Power-Operated flight control with a view 
'0 using the mass of the power-operated 


mechanism, particularly the motors, and to 
deploy that mass in a way favourable to give 
the stability which Professor Duncan had 
called for in the overall design of the wing. 


There was also the rather vexed point in 
the case of power-operated flight controls, 
as to the relative merits on the one hand of 
self-sustaining drive or non-reversible screws 
or of introducing damping, which was only 
another way of adding inefficiency and 
therefore weight and size to the operating 
mechanism or, on the other hand, of infinite 
stiffness,” or the approach to infinite 
“stiffness ratio” in the servo chain which 
carried out the power-operated control. 
Taking such a case as the one-tenth of a 
degree error between the torque shaft of the 
aileron or the rudder or elevator control 
surface, and the controlling signal setting, 
whether from manual or auto-pilot or gust 
detector, which had been visualised, a degree 
of “ stiffness” would be proposed such that 
this smal! error signal would produce the 
full restoring torque. That did represent a 
fairly high degree of “stiffness ratio” and 
would tend as he saw it, to set up a rather 
rapid correction, or set up a rather fine 
correction hunt about the ideal line as 
opposed to a slower coarser hunt which 
might be the case with damped controls or 
friction controls, or even manual controls 
using aerodynamic assistance. If they 
looked at the problem from that point of 
view they had to consider that rather ‘more 
rapid case as typical of the flight control 
problem which they were facing at the 
present time, and the chance it offered of 
reducing the power and weight by eliminat- 
ing inefficient screw damping and friction. 
What did Professor Duncan recommend 
from the point of using the efficient 
un-damped control as against the damped, 
or highly frictional control? 

Dr. P. Jordan: There were one or two 
points of minor importance which he 
thought should be mentioned in view of the 
importance which Professor Collar attached 
to the swept-back wing. Professor Duncan 
mentioned two effects which were essentially 
related. One was the fact that a wing, if 
its density were increased, gradually 
decreased its flutter speed, but only to an 
asymptotic limit. The other was that, con- 
sidering an aircraft flying at different 
altitudes and hence in air of different density, 
the smallest true critical speed of such an 
aircraft was its critical speed in dense air, or 
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Both facts were true in 


at a low altitude. 
general, but not in those cases where flutter 
with only one degree of freedom was 


possible. They knew of two such instances 
at high speeds. One of those, the case of 
low supersonic speeds, was founded on a 
doubtful theoretical basis and had not been 
confirmed by experiment. 


Professor Duncan mentioned that if the 
body of an aircraft with swept wings were 
allowed to oscillate, flutter might occur with 
a frequency which was related to the bending 
frequency of the wing but independent of its 
torsional frequency. But under certain 
circumstances the swept-back wing could 
develop pure bending flutter at high sub- 
sonic speeds, which meant that the wing 
which was infinitely stiff against torsion 
could flutter even if it were rigidly fixed at its 
root. This pure bending flutter depended on 
the wing being heavy compared with the 
surrounding air. In fact in all those 
questions the important quantity was the 
ratio between density of wing and density of 
air. As the swept wing at high speeds was 
held very stiff by the surrounding air, it 
could flutter only if its own density, or mass, 
were large enough to overcome this stiffness 
of the air. The damping forces of the air 
however, were negative in the present case, 
and thus endeavoured to make the wing 
flutter. Hence a swept wing, the density of 
which was increased, or which was taken to 
higher altitudes at the same true speed, 
might be in greater danger of flutter than 
it was before. 


Professor G. Temple (Department of 
Mathematics, King’s College, Fellow): He 
had a sound reason for being grateful to 
Professor Duncan—because he had been so 
kind to him. In his paper, in a masterly 
piece of understatement in relation to the 
theory of supersonic flutter derivatives, he 
had said that “the corresponding linearised 
theory for supersonic flow developed inde- 
pendently by von Borbely and by Temple 
and Jahn has not been confirmed by the 
measurements made recently at the National 
Physical Laboratory by Bratt and _ his 
colleagues.” 

The truth of the matter was that those 
calculations not only gave derivatives of the 
wrong magnitude, but of the wrong sign! He 
could only say that Professor Jahn and he 
had done their best, using the linearised 
theory, which completely prevented them 
from studying the effect of wing thickness 
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or incidence. As far as the linearised 
theories went they treated all wings as 
though they were flat plates oscillating at 
linear incidence. 


Looking through the paper he noticed that 
the British tradition, he might almost say the 
Teddington and Farnborough tradition, in 
flutter had been to establish the complete 
system of aeroelastic and dynamic theory 
and support it by initiating an extensive 
series of experimental measurements of 
derivatives. There did not seem to be much 
left for other people to do. But they knew 
that flutter had been studied in America and 
Germany and Italy and he thought it might 
be interesting if Professor Duncan added a 
little note at the end of the paper to indicate 
the kind of contribution which had come 
from those countries. 


R. H. Chaplin (Hawker Aircraft Ltd., 
Fellow): Much that they had heard at the 
lecture and which was in the printed paper 
referred to aircraft which existed. The 
aircraft designer’s problem was with respect 
to aircraft which did not exist. How was 
he to tackle the problems on complete 
new aircraft? What was the lecturer's 
advice to the aircraft designer as to the way 
in which he should go about it? It rather 
looked as if with the advent of the swept- 
back wing the empirical requirements for 
stiffness were even now less significant than 
before and it was therefore difficult to know 
what degree of stiffness should be built inte 
a design initially. Professor Collar talked 
about experiments as a means of examining 
the subject. He wondered whether small- 
scale models could be used to establish the 
stiffness and mass distribution which would 
be satisfactory in a new design? 


Professor Collar had told them that the 
aircraft designer should accept mass balanc- 
ing of control surfaces with good grace. 
Nevertheless there had been a number of 
successful aircraft in which there was no 
mass balance on some of the control surfaces 
and those aircraft had no record of any 
trace of flutter. There seemed to be some- 
thing wrong with present requirements 
which called for mass balance on such 
surfaces and which experience had shown to 
be quite unnecessary. 


Dr. R. A. Frazer (National Physical 
Laboratory, Fellow), contributed: A papet 
of this nature was badly needed, and would 
be welcomed equally by the Aircraft 
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FLUTTER AND STABILITY 


Industry and by all who were engaged in 
aerodynamic research. 


By way of introduction Professor Duncan 
alluded to a lecture on flutter which he had 
delivered to the Society just twenty years 
ago. That lecture was based, to a general 
axtent, on two years of flutter research which 
he had done with Professor Duncan’s close 
collaboration. One speaker had expressed 
disappointment during the Discussion at the 
poor progress that had been made with the 
subject. As the criticism had been worrying 
him for the past twenty years, it was 
reassuring to hear now from the present 
lecturer that even in those early days the 
foundation for fiutter theory was securely 
laid and that the principles of flutter 
prevention were already understood. 


He appreciated that it was much more 
important to know how they stood now in 
relation to flutter than how they stood in 
1929. The answer must be judged from the 
facts which Professor Duncan had so ably 
summarised. He had reviewed the essentials 
of the flutter and stability problems, had 
outlined the present state of development of 
the subjects, and had drawn attention to their 
growing inter-relationship. | He was also 
careful to point out that flutter was a trouble 
most liable to occur at high speeds, and that 
speeds of flight—already trebled since 1929 
—would continue to increase. Although his 
account made it clear that many useful 
advances had been made, he did not think 
that anyone would view the present position 
with complacency. On the one hand, aero- 
planes must be designed and built to fly at 
high speeds, and they must fly safely. On 
the other hand, the present trend in the 
analysis of aeroelastic oscillations was in the 
direction of increasing complication, while 
experimental work on flutter was becoming 
more and more difficult and costly. As this 
position would certainly not be remedied 
during the next few critical years, experi- 
mental flights at high speeds would need 
exceptional care. 


Professor Duncan’s lecture would do much 
to encourage the interest which many aircraft 
firms had shown in flutter investigation in 
recent years. He thought it was now widely 
recognised that many ad hoc flutter 
Investigations could be made advantageously 
by aircraft firms. The full use of such 
resources would relieve establishments more 
Properly engaged in fundamental or long- 


term researches and so would improve the 
flutter position generally. 


PROFESSOR DUNCAN’S REPLY 


One rather general theme had run through 
a large part of the Discussion, namely the 
request tor more experimental work on 
flutter, in particular on aerodynamic deriva- 
tives and especially, on full-scale. He was 
in complete sympathy with this. If he had 
any authority at ali he would put it one 
hundred per cent. behind that request. He 
did not know whether there was anyone 
present connected with the Ministry of 
Supply, but possibly there might be and 
possibly some echo of his remarks might 
arrive there, but he made a strong plea to 
that Ministry to increase the experimental 
facilities for work on flutter. He had been 
a member of the Oscillation Sub-Committee 
of the A.R.C. since it was formed and had 
been the chairman for quite a long time. He 
knew that the Committee had repeatedly— 
he would not like to say how many times— 
made requests for more wind tunnels for 
work on flutter and they had emphasised 
times without number, the urgent need for 
more work on flutter derivatives at very 
high speeds, subsonic and supersonic. 


Dr. Roxbee Cox: For a straight, uniform 
and unswept wing the “axis of independ- 
ence” would coincide with the spanwise 
axis through the aerodynamic centre. In 
other cases the agreement would not be 
exact but as a rule the difference would be 
slight. 


Professor Collar: Professor Collar had 
made an excellently clear and fair case for 
mass-balance of controls and he commended 
this viewpoint to designers. | He thought 
there might be some points about the 
technique of the use of air bags as supports 
for aircraft in resonance tests, which had not 
been completely worked out, but he knew 
that they had been successfully used at the 
R.A.E. and adhered to his view about their 
value. 

The most fruitful researches which could 
be made on the flutter of swept-back wings 
were, he thought, mainly experiments on the 
nature of the flutter of aircraft having swept- 
back wings and on the measurement of the 
derivatives for such wings. These researches 
were most urgently needed. He agreed 
whole-heartedly with Professor Collar’s plea 
for more experimental work on flutter. 
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DISCUSSION 


Professor Pugsley: He agreed that the 
approach to the flutter problem adopted by 
Sir Ben Lockspeiser in his paper “ A Simple 
Approach to the Wing Flutter Problem” 
(JOURNAL of the Royal Aeronautical Society 
July 1933) was in some degree akin to that 
of the lecture. The energy equation used so 
extensively in the lecture was developed in 
R. and M. 1155 (1928) where the importance 
of the phase relations of the motions had 
also been stressed. 


Mr. North: He agreed that the reaction 
time of the pilot was only one item in the 
total lag, but he could not agree that it was 
of little importance as it provided the 
measuring rod for all other lags, being itself 
irreducible. He also agreed on_ the 
importance of the frequency of any motion 
which the pilot sought to control. It was 
certainly correct that the pilot could not 
follow and control a periodic movement 
unless the frequency were low. The simple 
criterion which he had put forward was only 
intended to apply to aperiodic motions and 
to those of very low frequency. He remained 
of the opinion that the force response of the 
pilot was more important than his displace- 
ment response. 

In speaking of the “matching” of an 
auto-pilot and an aircraft he did not mean to 
imply that one type of auto-pilot might not 
be suitable for many different aircraft. He 
only meant that the auto-pilot should be well 
suited to the aircraft. He thanked Mr. 
North for the references he had supplied. 

Mr. Lyons: He agreed that resonance 
tests made in flight were a valuable method 
for investigating the flutter characteristics of 
a prototype aircraft and he hoped that this 
method would be used extensively. He 
welcomed the suggestion that the aero- 
dynamic derivatives of flutter might be 
measured on models carried by aircraft in 
flight and hoped Mr. Lyons would be able 
to follow this up. He agreed that he had 
somewhat over-stated the advantages of the 
power-operated all-moving tailplane, but 
considered that they were very substantial. 

Mr. Woodford: The interesting suggestion 
that an electric motor used to operate a 
control surface might be used as a balance 
mass appeared to him to be worthy of 
exploration. He regretted that he could not 
give a general answer to Mr. Woodford’s 
final question, as the success of any scheme 
depended so much on the ability and 
ingenuity of the designer. 
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Dr. Jordan: He agreed that his remarks op 
the influence of wing and air densities on 
flutter did not apply to that species of flutter 
caused essentially by the change of sign of 
a principal aerodynamic damping coefficient. 
He was interested to hear of a case of pure 
bending flutter of a swept-back wing of high 
density having rigid support at the root and 
would welcome further details. 


Professor Temple: Professor Temple had 
made humorous reference to the failure of 
the linearised theory of derivatives for super- 
sonic flow. He would like to emphasise that 
this was in no way the fault of the 
investigators; it was a fundamental fault in 
the basic assumptions of the theory and the 
failure was all the more surprising as the 
corresponding theory for subsonic flow was 
successful. 


It would be legitimate to say that British 
work on flutter had been carried on with 
vigour since 1924 or thereabouts. _ British 
investigators had always been very conscious 
of the fact that flutter was a problem with 
important kinematic and mechanical as well 
as aerodynamic aspects, whereas there had 
been a tendency abroad to place almost the 
whole emphasis on the aerodynamic aspect. 
Further, the British investigators had placed 
greater emphasis than most others on the 
experimental study of flutter. He would say 
that the most important contributions from 
abroad had been: Wagner’s theory of 
unsteady motion, Possio’s theory of the two- 
dimensional aerodynamic derivatives for 
subsonic flow, Kiissner’s emphasis on the 
importance of frequency in flutter, Stiider’s 
investigation of stalling flutter and von 
Schlippe’s work on forced oscillations in 
flight. 


Mr. Chaplin: He could not offer Mr. 
Chaplin any panacea of advice to the 
designer of new types of aircraft about the 
flutter problem, and could only say that, if 
the aircraft embodied any features which 
were not already well-tried in relation to 
flutter, it would be essential to have them 
thoroughly investigated. The investigation 
might begin with calculation and should end 
with resonance tests in flight. Any proposal 
to dispense with balance masses on controls 
should be dealt with in this manner. 


Dr. Frazer: He agreed that there was no 
cause for complacency with regard to the 
problem of flutter. No sooner was one 
problem solved than two or three new ones 
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He hoped that Dr. Frazer was right 
hat the lecture would do some- 
age firms to take an active 
n flutter investigations. In the past 


the Industry had depended to an undesirable 


STABILITY 


extent on the official research establishments 
to solve their flutter problems. It certainly 
should be the main business of the research 
establishments to evolve the methods which 
could be applied by the firms themselves. 
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THE GENERAL THEORY OF 
CYLINDRICAL AND CONICAL TUBES 
UNDER TORSION AND BENDING 

LOADS 


Single and Many Cell Tubes of Arbitrary Cross-Section with 
Rigid Diaphragms 


by 


J. HADJI-ARGYRIS, D.E., A.F.R.Ae.S., and P. C. DUNNE, 
B.Sc., A.F.R.Ae.S. 


PART VI—(concluded) 


Parts I to IV of this paper were published in the February 1947 JOURNAL, Part V 
in the JOURNAL for September and November 1947 and the first half of Part VI 
in May 1949. 


6.5. THE OPEN TUBE 


6.5.6. The open tube with no St. Venant torsional stiffness. The torsion-bending 
stress system gr hr 


As stated in the introduction to Section 6.5 the practically important case of 
an open section tube with no St. Venant torsional stiffness (see Fig. 38) will be 
treated as the limiting case of the open tube with J as J—>0. It will be shown that 
the stress system g,h, over the open tube degenerates into a stress system gr hr 
which is statically equivalent to the torque Ts about the shear centre. This stress 
system is for a uniform cylindrical tube with no J identical to the Wagner-Kappus 
torsion-bending stress system in open tubes. 

An important relation between A, and J of an n-boom open tube as J—>-0 can 
immediately be derived from the secular equation (562). 

Thus, 
as J—0, \1°7A,—>0 . : (591) 


so that the smallest root A, tends to zero and, 
lim == ByB.Ba | be, ed |*/4 B,ByB.By | ab, be, ed |?=1/1 
J30 (2) (592) 


where A, and A,, are given by equations (563) and (564). The justification of the 
notation I" will appear below. The result that A,” tends to zero as J is true, however 
large the number of booms and hence applies also when the walls have direct stress 
thickness f.. 

From the normalisation condition (78) and the equations (555) it follows that, 


Also from the normalisation condition (557) and equation (593) 
limJO,?7=1 . ‘ (594) 
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Walls carrying direct stress and shear stress. 
Walls carrying only shear stress. 
Fig. 38. 
Open tube with no St. Venant torsional stiffness. 


Equations (592) to (594) show that 4,9, tends to a finite limit as J tends to 
zero. The properties of this function will now be investigated. Multiplying equation 
(535) by ©, and letting J—>0 one obtains since ,°h,O, vanishes, 


ds ) =-§ lim 0, ) 
or with (594) 
djidy\ § 
where 
hr =-S*lim(h,9,) ‘ . (596) 
and 
p=limd, . . (597) 
J5>0 
The boundary conditions (555) and (556) become for the hr function 
[he sin = [ie cosvds=0 . (598) 
Cc Cc 
and 
| he pds=S* : ‘ . (599) 
Also at the free ends of the section, ° 
hy (0)=hr (S)=0 . (600) 
Integrating the differential equation (595) one obtains with conditions (600), 
dhy 


and, 


(602) 
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and the constant of integration in (603) is determined by the condition, 


In symmetrical cross-sections the origin of integration in (603) is at the centre of 
symmetry on the periphery. Re-writing equation (599) in the form, 


Cc Cc Cc 
dhy 
and substituting for — from (601), 


Cc 
which is the well-known formula of Wagner'*") and Kappus"*) for the torsion- 
bending constant I’ about the origin of p: In the notation of the above papers 2A, 
is denoted by w, which is proportional to the warping due to unit rate of twist about 
the origin of p. 


To find the origin €, 7 of p one uses the zero moment conditions in the form, 


dhr 
Cc 
With equation (601) and the iin. 


one obtains, 


8 


(p ( | )ds 
0 


(| ids Jas | 


Cc 
Inspection of equations (605) shows immediately that €, are the co-ordinates €¢s, 
nes Of the shear centre E,, which is in agreement with the Wagner-Kappus theory. 
Thus, 


Equations (597a) and (604) define now completely the fr distribution of 
equation (602). 

Denoting the limiting form of the stress-function F, = ¢,h, by Fr = grhr the next 
step is to consider the function g;, which by virtue of equation (596) is given by, 


The differential equation (299) in the first mode may be written, 


dg, = d Tr, 
or dividing by — S°0, noting that as J—>0, u,—>0 and one obtains, 


for the conical tube, 
d p dgr ] T3s 607) 
dr \t* dr )=- - 


and 
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and for the cylindrical tube, 

(Tr 

dz’ \t® dz’ ) ~ §2 ) (607a) 
To integrate these equations it is necessary to have the boundary conditions. For 
a tube built-in at r=r,, free at r=r., the boundary conditions (310) and (311) become, 


[gr] r=r, =0 
(608) 
dr |r=r, 
and the solutions of equation (607) are, 
__ Tos 
‘ ‘ . (609) 


If the boundary conditions are different from those stipulated above it is necessary 
to consider the general solution of equation (607) which may be written in the form, 


r 

1 Tx: 

where Dy is a constant determined from the conditions at the two ends of the tube. 
Using equations (602) and (609), which determine completely the stress-function 

Fr=go hy in a tube built-in at r=r, and free at r=r., in the expressions for the 

stresses given by equations (28) one obtains, for the conical tube, 


Tr r 


(00) 
_ _ _ Tas 0 | 


and for the cylindrical tube, 


4 


iy 


(61 1a) 
A striking analogy may be demonstrated between the above formule and the 
familiar expressions for the stresses in a built-in tube under shear loads. Thus to 
obtain the formule for the torsion-bending stresses due to T,s from the bending 
Stresses due to S,, read in formule (272) and (276a): — 
2A, for x 
for /, 


25 1,4,ds 


and in the case of a cylindrical tube 
T xs for S, 


| Tzsdz’ for M, 


and in the case of a conical tube 
Txs/p for QO. 


| 
| 
dr for M, | 


0 
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Only the first and third of these substitutions are basic, since the second and fourth 
derive from the first and third respectively. 
Thus, 


['=4 | 1.4.2ds I, =| 


Cc Cc 
r 


and M,=p [2 dr 


The first substitution shows also very clearly the interesting similarity between the 
following conditions, 


Cc 


t,Ads=0 and | 
Cc 


| 1.A.ds=0 and lr. ry ds=0 
Cc Cc 
In the derivation of the deformations of the open tube, which will be given in 
Section 6.5.11, some special limits occur which will be discussed here. 
Thus from the orthogonality and normalisation conditions (557) in an open tube 
with J one obtains, 


| ds =S°0, 
1-JO,? 
| ) 
Cc 
Similarly from the orthogonality and normalisation conditions (77) and (78), 
dts de | 
ds 


1 (dhy\?, 


Finally application of equations (284) to an open tube with J yields, 


hr hs 


Cc 
J t J>0 


6.5.7. The hr -function for an n-boom open tube 


The h - distribution and the shear centre E, of an open n-boom tube of arbitrary 
cross-section may be obtained either from the analysis of Section 6.5.3 by letting J 
tend to zero or from the general method developed in Section 6.5.6. 

By the first method the results are, 

[=4 & [B,B,B.B, | ab, be, cd |?]/ = [B,B.B,| be, cd|*7]. 
(3 (3) 
(hp)e. [B,B,B.B, | ab, be, cd | | pq, rs |]/2 [B,B,B.B, | ab, be, 


n 
4 


where x= ) ,l<a<b<c <d<n (see Fig. 36), and pq, rs, 
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are the two walls of the four-boom open tube (a, b, c, d) which do not include 
R, R+1 and which are taken in a clockwise order starting immediately after the 
wall including R, R+1. Thus if R, R+1 is included in be, | pq, rs| becomes 
icd, ab. 
The co-ordinates of the shear centre are given by, 

= —2 [B,B.Ba | be, cd | | be, cd |x]/= [B,B.Ba | be, ed |*] | 
(3) 3 


4 615 
Nes = (B,B.By | be, cd | | be, cd [B,B.B. | be, cd | 
(3) 3 


where L<b<c<d<n. 
If the shear centre has already been obtained with equations (605) or (615) it 
is preferable to use equations (602) to (604) for the determination of Ar and I. 


Thus, 


T=45 B, (Aar)’ . (604a) 
and 
282 R 
(hr )r, = — Rhy . . (602a) 


where the origin of the radius vector of the swept areas A., is at the shear centre 
E, and the point « on the periphery is determined by the equation, 


B.A,,/ = « . (6035) 
For a four-boom tube (see Fig. 39) formule (592a), (614), and (615) take 


particularly simple forms. 
Thus, 


4 2 


1 


Fig. 39. 
Open four-boom tube. 
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(617) 
A oss Asa A 412 
where =,, is given by equation (548). 
Also, 
Ags 
— | 23, 34 | 13, 34 | 12, 24 | 
4 2 
Aus ] s 
(618) 
23) 49 
| 12, 23 |,+ B) 


4 
Note that €:s and 7s in these formule are referred to any orthogonal co-ordinates 
x, 

For a singly-symmetrical section of the type shown in Fig. 37 one obtains by 
letting J and A—~>0 in Section 6.5.4, 


Nes = 2 =4 = ‘ . (619) 


The formula for I’ and Ar are most easily obtained from equations (604a) and 
(602a). Thus, 


=82 . : : ‘ . (620) 
1 
and, 
292 


In these equations the areas A,» derive from a an vector centred on Es. 

It is interesting to consider open tubes with J=0 which have no resistance to 
torque although built-in at one end. Obviously the three-boom open tube of 
arbitrary cross-section has this property. In the case of tubes with four or more 
booms or direct stress carrying walls the exceptional case arises when all walls 
lie on straight lines through a point which is clearly the shear centre. 


6.5.8. The eigenfunctions for an open tube 


The eigenfunctions for an open tube with J=0 are obtained from the general 
method set out in Section 6.5.2 with certain simplifications. Thus equations (556) 
and (557) are replaced by: 


Cc 
as=1 (557a) 


6.5.9. The eigenfunctions for an n-boom open tube of arbitrary cross section 


The eigenfunctions for an n-boom open tube of arbitrary cross section may be 
obtained immediately from the analysis for the n-boom open tube with J, given in 
Setion 6.5.3, by putting J=0. Thus the secular equation becomes, 


and has (n—4) roots A,. 
564 
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Equation (566) still holds and the co-factors Qp, Qn, Qn.,, and Qn,, are exactly 
as in equations (567) to (570). 
Hence the eigenfunctions and © are, 


1 
O=- ‘ . (623) 


n-1l 


The centres of twist €, and the values b,, by are as given in equations (573) and 
(573a). 


6.5.10. The eigenfunctions for a symmetrical multi-boom open tube. 


The results may be obtained at once from Section 6.5.4 (see Fig. 37). The 
eigenfunctions for flexure in the plane of symmetry YZ are the same as given in 
case (a) of Section 6.5.4. 

The eigenfunctions for twist and flexure in the XZ-plane are immediately 
obtained from case (b) of 6.5.4. Thus : — 

The secular equation is, 


The co-factors and » are the same and the eigenfunctions and © become, 
Q 
[2 = dp, p41 
(625) 


1 


6.5.11. The stresses and deformations for an open tube 
The expressions for the stresses are easily obtained from equations (587) by” 
letting J—>0 and noting that the first mode stresses degenerate into those of the 
l'-system. Alternatively they may be derived directly by the methods of Section 
5.4.4 using the stress function, 
F=gxhx+eyhy+erhr + . (626) 


and the first relations in equations (612) and (612a). By either method, 
(627) 


T=TetTr+ Gabe 


where +, and 7, are the generalised engineers’ theory stresses for loads through 
the shear centre. The G-functions are given for many loadings on a built-in tube 


| hrhy QO, 1 Vv 2 pe 
Gs | ; ds. | 2 dr+ | dr 5 * - (628) 
c r 


| 
) 
| 
| 
| 
| 
in Parts IV and V and are defined by equation (315a) of Section 6.5.5. ae 
The deformations of a built-in open tube may be obtained from equations (588) i 
and (589) by letting J—>0 and using the second equation of (612) and the relations i. 
(613). It is simpler, however, to use the method of Section 5.4.4 with the stress- i: fe 
function (626) and the first equation of (612) and relations (612a). By either method, Se 
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G 


Cc 


The deflection v, is obtained from equation (629) by interchanging X and Y and 
substituting —(&—€xs) for xs). The deformations for a cylindrical tube may 


r 


be written down at once. In the exceptional case when 0;=0 the substitutions given 
after equation (589) in Section 6.5.5 must be used. 

The meaning of the various terms in equations (628) and (629) will now be 
discussed. 


The first and second terms in (628) represent the twist due to the direct and 
shear stresses respectively of the torsion-bending theory. The third and fourth 
terms are twist resulting from the shear stress resultants Q, andQ, parallel to the f 
principal axes and acting through the shear centre Es. The remaining terms are 
the contribution to the twist of the eigenloads, which, as in a closed tube, is in 
general small compared with the other terms. Particular attention is drawn to the 
fact that the second, third and fourth terms can collectively be described as the 
twist due to the shear strains of the complete statically equivalent system. These > 
terms are neglected in the Wagner-Kappus theory. They are of considerable | 
importance in the relatively short open tubes encountered in aircraft structures. 
Note that in point symmetrical cross sections there is no twist due to Q, andQ, | 
and the third and fourth terms vanish, while in sections symmetrical about the 
X-axis, say, there is no twist due to Q,. 

Comparing the expression (629) with the corresponding expression (359) for 
the deflections of a closed tube one notes that, if E in (359) is taken as the shear 
centre Es, the forms of the two expressions differ only in the additional term 

¥ 
1 hxhr Tas 
GS | ds.p| dr 

Cc 
in equation (629). This term assumes importance in short open tubes. 

The warping of an open tube cannot be expressed solely in terms of the warping 
functions W,*. In fact there is necessary in addition a warping function which may 
be obtained as the limiting form, 


Wr*=-S lim(O,W,*) . . (630) 


where W,* is the warping function corresponding to the first eigenfunction in the 
open tube with J. Using equations (590) and (299) in (630) and proceeding to the 
limit, one obtains 


4 
s 


where gr is given by equation (609a) and W*;. is the warping at r=r,. For the 
particular case of a tube built-in at r=r, and free at r=r, equation (631) reduces 
with equation (609) to, 


jet? 
fo T> 


Note that these equations can also be derived directly from, 


Wrt= | hr ds= wds 
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where = is obtained from equation (26) with F as defined in (626). 
cs 


When considering the general boundary conditions in Section 6.6 it is necessary 
to have expressions for the actual warping w* of a tube. Equation (362) shows 
that, apart from rigid body movements, the warping w* of a closed tube or an 
open tube with J may be expressed as follows : — 

W* dh; 

i A? ts ds (63 ) 
In the case of an open tube the expression in the first mode becomes, using 
equations (601) and (630), 


w*= — 


and hence, w* /1 dh 


It should be noted that open tubes have a distinct difference from closed tubes 
in that they possess a certain kinematic freedom. Thus an open tube unconstrained 
against twisting may be subject to warping at one end proportional to the |" system 
direct stresses without inducing any stresses (see equation (631)). An imposed 
warping W,* at the end of an open or closed cylindrical tube dies out quickly 
together with the induced rate of twist, but an imposed warping Wr * at the end 
of an open cylindrical tube remains unchanged and induces a constant rate of twist. 
It is for this reason that the end conditions of an open tube have a much greater 
influence upon the deformations than is the case in a closed tube; indeed, the open 
tube can not resist torque and develop a statically equivalent stress system in the 
absence of suitable warping constraints. 


6.6. THE BOUNDARY CONDITIONS AT JOINTS BETWEEN TUBES WITH DIFFERENT 
CROSS SECTIONS ; 


In Parts IV and V the boundary conditions considered relate only to discon- 
tinuities in the lengthwise distribution of the material (i.e. discontinuities in ¢* 
and/or t,*) or loading. To represent wing structures with some degree of accuracy 
it is also necessary to consider tubes in which at one or more cross-sections there 
is a sudden change in the crosswise distributions f and/or ¢, of the material (e.g. 
cut-outs). The purpose of this chapter is to develop the boundary conditions for 
such cases. Each section of the tube in which ¢ and ¢,’ vary as in equations (38) 
and (40) is called an element and the cross-section at which two adjoining elements 
meet a joint. The analysis will be confined to the consideration of joints in which 
the direct stress-carrying material on one side is placed opposite some or all of 
the direct stress-carrying material on the other side (see Fig. 40). A necessary 
preliminary investigation is to determine the /j-functions, 0; and €, », in each element 
of the tube by the methods of Parts III to V. The calculation of the stress distri- 
bution is completed with the determination of the g,-functions from the given loading 
and boundary conditions. 

Section 6.6.1 develops the boundary conditions to give the necessary simul- 
taneous equations in the unknown D,-coefficients of the g,-functions in each element 
of the tube. The work may be shortened considerably by making use of the 
information available in the 05.08 series of the Stressed Skin Structures Data Sheets 
of the Royal Aeronautical Society. In these data sheets diagrams are given for: — 


(a) the G-functions of tubes, under certain important transverse loadings, 
built-in at the root and free at the tip, and 


(b) the G-functions of tubes under unit self-equilibrating end load systems. 
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(q-1)th element. 


Walls carrying direct stress and shear stress. 
————— Walls carrying only shear stress. 


Fig. 40. 


Part of wing structure composed of three joined conical or cylindrical tube elements 
q—1,q, and qg+1. 


Using the diagrams (a) the stresses and warping of each element are found, 
assuming each element to be built-in at its inner end and free to warp at the other 
end. Applying now the unknown multiples e; of the unit self-equilibrating end load 
systems, given in diagrams (b), necessary to equalise the direct stress flow and 
warping on each side of the joints, one obtains simultaneous equations for the 
coefficients e;. These equations have been developed in Section 6.6.2. 

In Section 6.6.3 the important particular case of joined four-boom elements is 
developed. 


6.6.1. The boundary conditions for the g,-functions 


Consider a tube built up by joining a number of conical or cylindrical elements 
each of which is governed by equations (38) and (40). The cross-sections of each 
element are arbitrary except that at each joint the direct stress-carrying material 
on one side is placed opposite some or all of the direct stress-carrying material on 
the other side. The taper of each element may also differ. Before developing the 
boundary conditions some clarification of the notation used will be given. Thus:— 


(a) All quantities belonging to the gth element have the suffix q. 
For example, 


hu, Pasir teas Mx, etc. 


(b) The values of quantities belonging to the gth element at the joints (y— 1. q) 
and (g, g+1) are denoted by the additional suffices c and o respectively. 
This notation is consistent with the previous use of o for the root and c 
for the tip, and the reference cross section is at end o. 

For example, 

81,941, Ci T sigo> M,, a+1, o€te. 


Section A-A : Section B-B Section C-C 
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Note that poo =Pa+1,o=1 and that s, is measured over the periphery of the gth 
element at o. The asterisk in the warping function W; has been omitted in this 
section for ease of printing. 

The boundary conditions at the joint (g, g+1) will now be investigated, 
assuming for the moment that each element has at least one closed cell or a-torsion 
constant J and that the (q+ 1)th element has direct stress-carrying material over at 
least as much of the periphery as the qth element (see Fig. 40). Thus the number 
Ny. of eigenmodes j of the (q+ 1)th element is greater than, or equal to, the number 
N, of eigenmodes i of the gth element, i.e. 


The first set of boundary conditions expresses the fact that at every point of 
the periphery of the joint the direct stress flow must be equal on each side. Hence. 


foo=fa 
or 
dhyx, dhyaq v dhiq | 
1 dhx,q +1 


Note that »,,=1, but that it is convenient to retain the expressions p,, in the 
equations as the resulting symmetry is of use later. Now multiplying both sides 
of this equation by, 


and integrating w.r.t. s,., over all walls of the (q¢+1)th element one obtains N,., 
equations of the form, 


M yqo M,. oO 
Pao qo i=l 
where, 
ts dh 
tag anj,q+1 H 
Yo (Ai, a5 ds,., 
and 


1 Cdhg | 


If the gth and (¢+1)th elements are multi-boom tubes with walls effective only in 
shear equations (636) and (637) become, 


1 ma B 2 
(Aj, q- 1) =1 Bas l,m (636a) 
1 
= m A j,q+1,m 
Yu (A;, m=1 Bas 1. m h } 
and, 
(A;, q- i¥ in=1 B, +1,m 


Where n, is the number of booms of the gth element. Note that the summation in 
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relations (636) extends only over those booms that are common at the joint 
(q, q+ 1). 

If the shear walls of the gth tube are covered by the shear walls of the (¢+1)th 
tube, equation (637) may also be written, 


qa+1 
The first two terms on the left hand side of equation (635) may be expressed 
in the alternative form, 
1 f 1 q+l1 


where fi, is the direct stress flow at the end (g, 0) due to the engineers’ theory. 
In the case of multi-boom tubes this equation becomes, 


1 a+1 


where (Peqo)m 1s the engineers’ theory load in boom m at (q, 0). Note that AAj, q+1, m 
may also be written P «1, m, Which is the mth boom load for unit g; in the (q+ l)th 
tube. 

In the special case when the direct stress-carrying material on both sides of 
the joint (g, g+1) is equal or proportional, so that f.,/f., q.1=const., 
then 

and equation (635) simplifies to, 


A 


q 


M 


i 
If in addition the cross-section geometry is the same and, 
=Const. 
equation (635a) reduces to the particularly simple form, 


which is a generalisation of equation (312). 

The second set of boundary conditions expresses the fact that the warping 
displacements w must be equal at points of the joint where direct stress-carrying 
material exists on each side. Hence, at these points, 


Using relation (632) in (638) one obtains, 
( 1 o«) 1 o+1 
t 


ll 


) +rigid body movements 


i=l (Ain)? sq ds, ts,q41 
(638a) 
Multiplying this equation by Be and integrating over all walls on? obtains 
the N, equations, 
Nq4+l 
j=1 


where ci; is given by equation (637). Note that although equation (638) applies only 
over the direct stress-carrying material of the gth element the integration on the 
right hand side of (639) has been performed also over the direct stress-carrying 


material of the (q¢ + 1)th element; this is correct since d ng is zero over these regions. 
q 
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Re-writing equation (639) with equation (364), and multiplying by G, 


| (Pao) dr le | (pu dr 


If there are no further discontinuities and the tube is built-in at (¢+1, 0) and 
free at (g, c) there will be (V,,+N,.,) equations of the form (639) and (635), Ny 
conditions (310a) and N,,, conditions (311a) to determine the 2 (N,+N,,,) unknown 
D’sin the gi, and g;,4,1 functions. Similarly, if there are more than two different 
elements it will be found that equations of the type (635) and (639) provide the 
necessary number of conditions to determine the unknown D’s. 

If the opposite conditions apply to those described by equation (634) and the 
preceding paragraph, that is if, 


equations (635), (636), (637) and (639) have to be rewritten by interchanging, 
i and j | 
gandqg+l;} . ‘ . (640) 
o and c 


The important practical case when tiere 1s a completely open element on one 
side of a joint will next be considered. It will be assumed that the gth element is 
an open tube of the type discussed in Section 6.5.6, the (q¢+1)th element has at 
least one closed cell, and that condition (634) holds (see Fig. 41). In this case the 
function g,, with two unknown D’s of the case just considered degenerates, as J—>0, 
into the function g;, with only one unknown Dr« (see equation (610)). A second 
unknown in the |’ system now appears as the W;.,,. of equation (631). This difference 
arises from the kinematic freedom of an open tube discussed in Section 6.5.11. 
The boundary conditions may be obtained by a limiting process from equations 
(635) and (639), which hold for an open tube with J, or ab initio, using in the gth 


(g-1\th element 


qe] 


Section A-A Section B-B Section C-C 


—_——. Walls carrying direct stress and shear stress. 
———— Walls carrying only shear stress. 


Fig. 41. 


Part of wing structure composed of three joined conical or cylindrical tube elements of which 
the central one q is open. 
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element the stress function (626) and the warping expression (633). By either method 
the equality of direct stress flow condition yields the N,., equations of the type, 


qo 
M, My Naq 

qo qo q i= 

Tac 
where 


If the gth and (q¢ + 1)th elements are multi-boom tubes with walls effective only 
in shear equation (642) becomes, 


l 
ipl, Ba +1.m 
The third term on the left hand side of (641) may also be written, 


2A 1, m . . (642a) 


or for a multiple-boom tube, 


i PraoOhi.a +) 


where frao and (P,..),, are the direct stress flow and the load at the boom m due to 
T ys in a tube built-in at r=r, and free at r=r.. 

Similarly the equality of warping yields the (N,.,) conditions of which the N, 
are exactly as in equation (639) and the last condition is, 


Na+ 


j=l 


If there are no further discontinuities and the tube is built-in at (¢+1, 0) and 
free at (g, c), Di. =0 and equation (643) is not required; in fact in such a case the l° 
stress system is statically determinate. There will be in all N,+N,., equations of 
the form (639) and (641), N, conditions (310a), and N,., conditions (31la) to 
cetermine the 2(N,+N,.,) unknown D’s in the g;, and gj,,+1 functions. 

When the (q+ 1)th element is open and the gth element is closed and condition 
(634a) holds, the conditions at the joint (g, g+1) may be obtained from equations 
(641), (642), (639) and (643) with the substitutions given in relations (640). If there 
are no further discontinuities and the tube is built-in at (¢+1, 0) and free at (q. c), 
Wr.gs1. =O and there will be 2(N,+N,.,)+1 equations to determine the 
2 (N,+N,.,) and the coefficient Dr, q +1. 

Similarly, if there are more than two elements, it will always be found that 
equations of the type (641), (643) and (639) provide the necessary number of 
conditions to determine the unknown coefficients. 


6.6.2. The boundary conditions expressed in terms of the undetermined 
e-coefficients and the G-functions 


As stated in Section 6.6 the stresses in each element may be expressed as the 
stresses due to the applied external loads with the end o built-in and the end c 
free, plus the stresses due to self-equilibrating end load systems. 

Thus the stresses in the gth element, assuming that there is at least one closed 
cell, may be written, 


N N 
q qa 
+ = T + HiqLaig CiqoG cial | 
i= i=1 
(644) 
Ta=Tpa + = + 


a c o 
HMiqLriq CiqoGria] 
1 
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The G-functions are those corresponding to the applied load with the end o 
built-in and the end c free. These functions are given for various loadings in 
Parts IV and V and also in the 05.08 series of the previously mentioned data sheets. 


The G and G functions are those corresponding to unit self-equilibrating end 
load systems applied at the ends c and o respectively with the other end free. These 
functions are given in Section 5.5 and in the data sheets. The L-functions are 
defined in equations (340). The coefficients e;,. and ej,, are the unknown multiples 
of the ith unit end load systems applied at the ends c and o respectively. 


The warping function W;, due to the applied load with the end o built-in and 
the end c free is given by equation (590), 


T 


The functions W;, and W,, are the warping functions W,* of equation (372a) 
for unit end load systems at the ends c and o respectively. Thus, 


Hig 


Pa 
G ‘| 
oO 
| 
Note the reciprocal relation, 
The boundary conditions (635) and (639) can now be re-written in the form, 
Yu — sao T iqG siqo + Cigi@iqo = @ i. . (647) 
Pao Pgo i=1 i=1 ; 
CigeW igo Cig (C5, 41, Wi, o41,¢ +3, q~+l, ) 
i=l j=l 


(648) 


The only unknowns in these linear simultaneous equations are the e-coefficients. 
If the end (qy, c) is free the coefficients e,,. are zero and only the coefficients e,,, 
remain unknown in the gth element. If the end (q+ 1, o) is built-in there is no stress 
condition (647) and the warping condition (648) reduces to, 


If the gth element has more eigenmodes than the (q+ 1)th element the appropriate 
equations may be obtained from (647) and (648) with the substitutions given in 
relations (640). 


The boundary conditions will now be given for the case when the gth element 
is an open tube. The stress expressions (644) become for an open tube, 


2A, 
Pat sq 


[€iqG siq + (650) 


N 
q c 
45 2 High Tiq [@igeGriq CiqoGriq] 
=1 


Where o, and rp, are given by equations (611). 
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The condition (641) for equality of stress flow becomes, 


Tao 


N 
M, M, Cini 
xqo BSq Nv iqij 
Pao ao Pa i=1 Mig 
Tac 
x Cigi@iqo = €j, a4 . . (651) 


The N, + 1 warping conditions at the joint (g, g + 1) are given by the N,, equations 
(648) and equation (643) which in the present notation becomes, 
(Nas Na+1 


Nq+1 c 


In the simultaneous equations (651), (648) and (652) the unknowns are the 
e-coefficients, Drg and Wrao. When the (¢+1)th element is open and the gth 
element is closed and condition (634a) holds, the conditions at the joint (¢, q+1) 
may be obtained from equations (651), (648) and (652) with the substitutions given 


in relations (640). The left hand side of equation (652) becomes Wy... -1,., which 
however, can be expressed in terms of Dp, q.1 and Wy, «-1.. with relations (631) and 
(609a). Thus, 
To is r 
fo To 


where the suffix (¢+1) has been omitted for clarity in printing. Note that the 
unknown W;,. appears only in the two equations of the type (652) for the joints 
(g, q+1) and (q—1, q) and can hence be eliminated. The resulting expression for 
Dy. in terms of the e-coefficients can be substituted in equations (651). Thus the 
only remaining unknowns are the e-coefficients as in the case of the closed tube. 


If the (¢+ 1)th element is built-in at (¢+1, 0) conditions (649) hold and 


The practically important boundary conditions in a tube built-in at o over part 
only of the direct stress-carrying material requires special consideration. The two 
boundary conditions are that there is zero direct stress over the part C, of the 
boundary not built-in and zero warping w over the built-in part C. of the boundary. 


Assuming that the tube is open the zero stress condition may be written in 
the form, 


i 


over the part C, of the boundary. 
The zero warping condition is, 


c 9 
Wr 6 +e,,W,,]+rigid body movements=0 (655) 


over the part C, of the boundary. 


Conditions (654) and (655) are always sufficient to determine the unknown 
coefficients. Thus, for example, in an n-boom tube with m booms fixed, where 
m > 4, there will be (n— m) conditions (654) and m conditions (655) which, together 
with the (n—3) conditions at the end c, are sufficient to determine the 2(n—4) 
coefficients, e, D;, Wy., and the three arbitrary constants in the rigid body dis- 
placements at o. To eliminate initially the rigid body displacements from equation 
(655), the following procedure may be used. Consider the closed tube obtained 
by joining the two ends of the open tube C, by an arbitrary wall effective only in 
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shear. For this tube calculate the functions h;. Multiplying now equation (655) by 


a and integrating over C,, the integrals over the rigid body movements vanish and 
s 


one obtains for each j equations of the form, 


c 


where 
Af 
ds 
(657) 
ds ds 


The case when the partially built-in tube has at least one closed cell is obtained 
from the above equations by excluding the |’ terms. Note that if only three booms 
are fixed in a closed tube only equations (654) are required to determine the 
unknowns. 


6.6.3. The boundary conditions for the four-boom tube 


The boundary conditions for joined four-boom tubes may be obtained from 
Sections 6.6.1 and 6.6.2 by omitting the summations signs in the various boundary 
relations derived there. In the case of the four-boom tube, however, it is preferable 
to use in the boundary conditions the boom load function P, defined in equation (233), 
Then the appropriate warping function is w such that }Pw is the work done by an 
end load system P. Note that w is more suitable than W, of Section 6.6.2 and W* 
of Section 4.3.4, since for given boom displacements it depends only on the positions 
of the four booms and is the proper generalised co-ordinate to associate with the 
generalised force P. The various warping functions are connected by the relations; 

Pw = PW* /4A*p?= — 2,W,* 
but from (229), 


4 
1 
Hence for closed tubes using (199), 
4 = 
w= W*/44%p? . . (658) 
1 
and for open tubes using (617), 
Wr =W*/4A*p?= —Wy *=,,/A*S (6 58a) 


The case when all joined elements have at least one closed cell will first be 
considered 

If P,, and P,. are the unknown boom load functions at the root and tip 
respectively of the gth element the stress equations (644) may be written, 


(659) 
“1 Ba + TEq T [Pu.Grq 
where, 
E 1/2 4 (A, 
NG = 2A, “= ) 2 Is 
*/u A, [Is 


and W is given by equation (549) for a single cell or a multi-cell tube without 
openings and by equation (549a) for a multi-cell tube with an opening. 


) 
) - 
a 
) 
} 
i 
“i 
oo 
575 


J. HADJI-ARGYRIS AND P. C. DUNNE 


Note that the last expression in the first of equations (659) may be written for 
boom 1, for example, 


1 e 


B, 
in which form the relation is obvious. 
When ©, =0 one replaces in equations (658) and (660) by 
(Qxqob QyqoP a) 
The warping function w as defined in (658) will now be given for various 
loadings. Using equations (365) and (658) the warping function w, for a tube 
built-in at (g, 0) and free at (q, c) gi 


T 7 
T poof aG'rq | (661) 
where ©», is defined by, 
Thus for a single cell or mullti- cell element without openings, 
and for a multi-cell tube with an opening (41), 
. ; . (663a) 


The warping functions w,* for unit g are given by equation (372a). From this 
equation, using the conversion ance between P and g given in (233), and equation 


(658), the warping functions » w and w corresponding to P,=1 and P.=1 are, 


Note the reciprocal relationship, : 
Wac= — Wao . : : : . (664a) 


Using suffices E and T to denote engineers’ theory boom loads and axial 
constraint boom loads of an element, built-in at o and free at c, the corresponding 
boom load functions become, 


Py= Pr, — Pr. + Prs — Pr, =T3.9AG;. . (666) 


The condition of equality of boom loads on each side of the joint (g, g+1) 
may now be written in the form, 


=P Pro + — Pe, o41,-)= 1, (667) 
This is a special case of equation (647). Ata free end (qc), 


If the boom areas in the elements g and (q+1) are proportional, or if the 
elements are under pure torque, 


The equation for the equality of warping at the a (q, q+) is, 
and at a built-in end (go), 
Py Woot Pa (668) 
These equations are special cases of (648) and (649). 
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The solution of equations (667) and (668) yields all the unknown P’s and hence 
the stresses by equations (659). To solve these equations it is best to express the 
P,,’s in terms of the P,.’s by means of (667) and substitute for the P,.’s in (668). 
Thus one obtains “three moment” equations in the P,.’s (see equations (677)). Note 
that if the G-functions are available for both ends free equations (667) and (668) 
take a more symmetrical form, 


and 
Waot Was — Py+1, Was =Waq41,¢ — Woo (668b) 
Note that if Pra. —Pr,q+1.. =0 


and equation (6685) becomes a “three moment” equation in the Po,q+1’s. For 
rectangular and trapezoidal cross sections this simple form has been given by 
previous writers'®: 

The next step is to give the formule for the stresses and warping of a four-boom 
open element and to investigate the conditions at the joints of such an element with 
a closed element. 

The boom loads and shear stresses in an open four-boom element may be 
written in the form, 

+ Prm + Pom 


Tm, ms+1 =—TEm, m+1 + m1 ) 


where the suffix » refers as usually to the mth boom. Pr and rr are the torsion- 
bending boom loads and shear stresses in the element, when built-in at o and free 
atc. P.., is the constant load in the boom m derived from the unknown boom load 
function P. The relations between P.,, and P, are given in equations (232) and (233). 
The connection between P. and the constant of integration Dr of Section 6.5.6 is 
given by, 


(669) 


r+l 
if the opening is between booms randr+1. Using equations (601), (602), (611), (616) 
and (617) P,.. and zp for an element with an opening (41) are given by, 


Pr: Prs = Pri Pr dr 
where =,, is given by equation (549), and 
These equations may be compared with (232), (233) and (617). 
The warping function wy as defined by equation (658a) is, 
Tr T T 
P. | 1 {7 
Wr =Wr, 4E(A pte EAT) pte) p drdr . (672) 
To T> Ty 


Where the suffix |” in this equation has been used to denote the total warping of an 
Open element. wy, is the unknown warping at the end o and I" is given by 
equation (616). Equation (672) may be derived using equations (653) and (670) 
in (658). If the element is built-in at o and free at c the warping function wr is 
given simply by tie last expression on the right hand side of (672). 

If the missing wall is (12), (23) or (34) instead of (41), equations (671) are still 
applicable if all suffixes are cyclicly advanced by 1, 2 or 3 respectively, but 
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<----l-- —---st-- 
4 z va 
q+3 q+2 q+l q q-l q-2 
Fig. 42. 
Diagrammatic representation of joined tube elements. 
equations (672) must be transformed by replacing =,, by —=,., =., or —=,, 


respectively. 

Consider now a series of four-boom elements as shown in Fig. 42 in which the 
qth is an open element while the (q—1)th and (q+ 1)th have at least one closed cell. 
The following conditions at the joints (g, g+1) and (q¢—1, qg) may be found either 
directly or from the general analysis in Section 6.6.2. 


Joint (q, q+1) 
Condition of equality of boom loads, 
where P;,,, is the value of Py, at the end (qo). 
Condition of equality of warping, 


c ) 


Joint (q-1. q) 
The two relations are, 
i,c l,o +P, i,o =WTae . . (674a) 
Note that the only unknowns introduced through the gth open element are P,. and 
Wrap SINCE IS related to wr,,, by equation (672). 
Adding (674) to (674a) one obtains with (672), 


P, 1,cWa 1,0 +P, 1,0Wq 1,6 =P Wig oWig = Wa+1,¢ + Woe 
. (675) 
where, 


Ty 

To 


(676) 
To T 
1 Sa | 1 {2s 
wy To 


are the warping functions at the end (gc) due to unit P,. and the loading Txs 
respectively with the end (go) built-in. 

To illustrate the application of the above boundary conditions to the formation 
of the “three moment” equations in the P.’s, the four equations for the equality of 
warping at the joints (q—3, g—2), (¢q—2, q- 1) and (q+1, q+2) and the difference 
of warping between the joints (q, ¢+1) and (q—1, q) will be written down. 

Note that no separate equations are required for the joints (q—1. q) and 
(q, g+1) since P,. can be eliminated by (673). 
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Joint (q-3, q-2) 
c 
+(Wa-3,0 —Wa-2,¢ ) Pa-2,¢ — Wa—2,¢ Pa-1,¢ =Wa-2,¢ 


Juint (q-2, g-1) 


From difference of warping between joints (q—1, q) and (q, q+1) 


Joint (q+ 1,q+2) 


c c 
Wa+l, +(Wa41,0 — Wa+2,«) — Wa+2,¢ Pus =Wa42, 


Note that (677a) and (677d) do not involve the open element q and are identical in 
form; they are the standard “three moment” equation when all elements have at least 
one closed cell. The two equations (677b) and (677c) both involve the open element. 
The II functions are defined by equations (667) for the closed elements (q—3, q—2, 
q-1.q+1, q+2) and by equation (673) for the open element g. The coefficients w 
are given by equations (661) and (665) for the closed elements and by equations (676) 


for the open element; these latter coefficients are Wye and W,- in the above equations.. 
Near built-in or free ends equations (677) require modification. 


Built-in at (q +1, 0) 
Equation (677d) immediately reduces to, 


Note that equation (678) is identical with (668a) since from (667) 
= Po+2,¢ —II,., 


It is preferable, however, to keep P,2,- as the unknown since (678) then appears 
as a particular case of (677d). In this case II,,, may be fixed arbitrarily but is 
conveniently taken as 


Pr, a+l,o 
The remaining “three moment” equations are unaltered. 


Built-in at (qo) 
Equation (677c) becomes, 


Ww, 1,0Pq +(Wa-1.0 —Wec) Pati, ¢ = Wye + Wa-1, 0 + (Wa_1, 0 — Wace) Tq (679) 
II, may be taken as Pi). The remaining equations are unaltered. 


Free at (q—1, c) 
Instead of (677b) one has, 


which replaces (677c). The remaining equations are unaltered. 
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Free at (qc) 
In this case P,.=0 but from (667) 


which replaces (677c). The remaining equations are unaltered. 

In practice when stressing a wing the root may always be taken at the aircraft 
(L. If the loading is symmetrical the root is built-in. If the loading is anti- 
symmetrical the root will be free to warp and under no direct stresses. Denoting 
the root element by n the condition at the joint (m, n+ 1), whether element n is open 
or closed, is 

This is immediately obvious but may be obtained as a particular case from (667) or 
(673) noting that 
Pr, n+1l,c =P,..=0 
and 


If q=n—1 equation (677d) is replaced by (682), and if q=-n equation (677c) is 
replaced by (682). 


6.7. SOME FURTHER RESULTS ON THE THEORY 


This section gives some further developments of the theory which are useful 
for practical applications. 


6.7.1. The hy and hy distributions and the shear centre in multi-boom tubes 


A general method of finding hx, hy and the shear centre Ex is given in Section 
5.4.3. For single-cell n-boom tubes (see Fig. 25) and multi-cell n-boom tubes of 
the type shown in Fig. 16 it is possible to give immediately symmetrical formule 
for hx and hy when referred to Es. 

The hx and hy functions for any orthogonal system of axes through the centroid 
are defined by equations (275) and (275a) which may be written, 

where 


s 


For a single cell n-boom tube there will be n equilibrium conditions at the 
booms, of which the mth is, 


Only (n—1) of these equations are independent. Since hy is here referred to the 
shear centre there is no twist and hence, 


n 


= Om, m- iD xm, m+1 > (685) 


m=1 


The symmetrical solution of equations (684) and (685) is, 


n n 
m+1 >= > (6m Bese) /2 > by, r+1 (686) 
r=1 r=1 
and similarly, 


n 


Dia: m+1 = (Om Onn) B.y,/2 > Or, (686a) 
r=1 r=] 


from which hy and hy are given by (683). 
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It is often preferable to calculate D,,, and D,,, from the non-symmetrical 
formule, 


r=2 r=1 

and, n n 

r=2 r=1 


and to obtain the D,., and so on from the recurrence formule (684) and the 
responding one, 

When the tube is open between booms n and 1 formule (686) reduce to the obvious 

forms, 


m 


n+1 = = B. Yr (688) 
r=1 


m 


1 


The co-ordinates ss. zs of the shear centre may be determined from the 
equations, 


Ty pD,ds= 2 Am, m+1Dxm, m+1 
m=1 
(689) 
| = = pD,ds = —2 Am, m+1Dym, m+1 

For the derivation of these equations the reader is referred to the corresponding 
case of an open tube given in equation (605). Note that p and Am,m+1 may be 
referred to any origin and that ;s, xs will then be given w.r.t. axes through this 
origin parallel to the x, y axes. a 

For multi-cell tubes of the type shown in Fig. 16 a procedure similar to that 
given in Section 4.3.3 for the /g-distribution can be applied. Thus putting, 


D° xm, mer D' xm, n+1 = Dxm, m-+1 . (690) 


the equilibrium condition at the mth boom again appears as in equation (684). The 
zero twist condition of the cell between booms m and (m+ 1) may be written, 


Solving (690) and (691) for D°® and D' one obtains, 
D° xm, m+) m+1 m-+J Din, m+1 . . (692) 


The zero twist condition for the whole tube is, using (692) 


n 


m+ 1 D°xm, mil >= Su, m+1 Dxm, m+1 = (685a} 
m=1 m=1 
where © is defined by equation (206a). Thus the D,’s are again given by equation 
(686) using the 4’s of the equivalent single-cell tubes of Sections 4.3.2 and 4.3.3. The 
shear centre is readily found to be defined by, 


m=1 


n 


(689a) 

— = A‘ m, m+ iDsm m+1 

m=1 ‘ 

where the A‘’s are given by equations (207). For the multi-cell tube with an opening 

between booms n and 1 the D,’s of the equivalent tube are again determined by 
equation (688). 
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6.7.2. An alternative method for finding the centres of twist S, 


When the shear centre is known it is often preferable to use equations (284) for 
the determination of the co-ordinates &, 7; of S; 


— — Ss) +1, (n,- Nes)= ds 


hy 
— Ens) — Ixy OF — Nes) = ds 


In the case of an n-boom multi-cell tube of the type shown in Fig. 16, the integrals 
may be taken over the walls of the equivalent single-cell tube. For a four-boom | 
single or multi-cell tube the equations for &;, 1; become, | 


—Eps)+1, Nes) = 2 r+ iA Pp, p- r+1 


(693) 
i; (E, — Ens) — Ixy 1A +1Dyr, r+l [9p 


where p=r+2 and ©), is given by equations (663) or (663a). Note that for a multi- 
cell tube with an opening the term in the summation corresponding to the missing 
wall must be omitted. If 0,)=0, while the numerators on the right hand side of (693) 
are finite, the terms b, and bh, are obtained from, 


4 


r=1 


(693a) 


4 
= +6 r,r+l A Pp. ri | 
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6.8. NUMERICAL EXAMPLE 


This section presents the calculation of the stresses in a wing the half-span of 
which consists of seven joined four-boom conical or cylindrical elements as illustrated 
in Fig. 43. The overall dimensions and the loading (Fig. 45 and Table II) are based 
on an actual wing under a symmetrical gust but some changes in the structure have 
been made to make the example more illustrative of all types of four-boom cross 
sections. Each element has been adapted to satisfy the condition r* =1,* = 1 and for 
this reason there are discontinuities in the boom areas at the joints. The adaptation 
of actual boom areas to satisfy a t,* law is discussed in the Appendix. The given 
boom areas also include an allowance for the direct stress-carrying capacity of the 
covers. Methods of making such allowances appear in the Appendix. 

The geometry of each element is given in Table I and the cross-sections of 
elements 5 and 6 are illustrated in Fig. 44. Table II gives details of the loading and 
the position of the axis of reference for the torque T, . The calculations are arranged 
in tabular form and with the notes given are self-explanatory, with few exceptions. 
Detailed calculations are not given for every element but sufficient are included to 
illustrate the numerical procedure for each type of four-boom cross-section. Thus, 
for example, the detailed calculations of the cross-wise variation of the stresses are 
given for element 6 in Tables IV, [Xa and IX), but since element 4 has the same type 
of cross-section no such detailed calculations appear for this element. 

The shear stresses in the walls of the tube are sometimes in excess of the buckling 
stress and consequently, when calculating the 6’s, it is necessary to assume an 
effective wall thickness by reducing the actual wall thickness in the ratio G,/G 
where G, is the secant shear modulus. The ratio G,/G will vary both around the 
cross-section and along the tube and it is clear that only by successive approximation 
can the correct result be obtained. Such a procedure would be extremely lengthy. 
In the present example the actual thicknesses are used throughout and over the more 
heavily loaded elements 2 to 6 the shear modulus is taken as G,=0.8G; in element 
one no reduction in G is necessary. Since all the shear and torque is removed at the 
fuselage pick-up plane there is no shear in element 7 and hence G, need not be 


considered. 


The calculations have been made in terms of the boom loads and the shear 
stresses. The use of boom loads rather than stresses is preferable since they are in 
any case needed to form the “three moment” equations and afford a ready check 
at various stages of the calculations. All the derived cross sectional properties are 
computed to five significant figures. Such accuracy is not necessarily desirable in 
the stress office. The computations are not in general sensitive to small errors and 
slide rule accuracy is sufficient. However, if using a calculating machine, it is worth 
while retaining an extra figure since the checking processes are then more positive. 

All quantities, unless otherwise stated, are in ton inch units. Young’s modulus 
E is taken as 10’ ib./in.* and the shear modulus G as 3.85 x 10° Ib. /in.*. 


A list of the tables and figures follows with further explanatory notes where 


required. 


Fig. 43. 
Fig. 44. Cross-sections of elements 5 and 6. 


Fi 


g. 45. 


Plan view and cross-sections of wing. 


Lift loads and torque on wing due to symmetrical gust. 


Tables Ia and Ib. Basic geometrical properties of tube elements. 


Table 
Table 


Table 


Table 


II. 
IL. 


IV. 


External loading on wing due to symmetrical gust. 

Elements 1 to 7: Cross-section functions for the engineers’ 
theory boom loads. 

Element 6: Calculation of cross-section functions for the 
engineers’ and Bredt-Batho shear stresses. 

Element 5: Calculation of cross-section functions for the 
engineers’ shear stresses and the torsion-bending stress system. 
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Table VI. Element 1: Calculation of cross-section functions for the 
engineers’ and Bredt-Batho shear stresses. 

Table VII. Elements 2, 3 and 4: Cross-section functions for the engineers 
and Bredt-Batho shear stresses. 

The term “cross-section function” in Tables III to VII refers to those factors 
in the expressions for the boom loads or shear stresses which depend only on the 
root cross-section geometry of each element. 

Table VIII. Co-ordinates of shear centre. 

Table VIIla. Relative co-ordinates of shear centre with respect to reference 

point E and the zero-warping centre. 

Note considerable difference in positions of shear centres and zero warping 
centres; in element 3 the zero-warping centre is in front of the front spar. 

Tables [Xa and IXb. Element 6: Calculation of cross-section functions for the 

axial constraint boom loads and shear Stresses. 

Table Xa and Xb. Element 1: Calculation of cross-section functions for the 
axial constraint boom loads and shear stresses. 

Table XI. Elements 2, 3 and 4: Cross-section functions for the axial 
constraint boom loads and shear stresses. 

Table XII. Element 2: Calculation of lengthwise functions for boom loads 
and shear stresses of the statically equivalent stress system, and 
torque Tx,. 

Table XIII. Elements 1 to 7: Lengthwise functions for boom loads and 
shear stresses of the statically equivalent stress systems, and 
torques 

The term “lengthwise function” in Tables XII and XIII refers to those factors 
in the expressions for the boom loads or shear stresses which depend only on the 
loading and position from the root. It is convenient to approximate 7», by 
combining certain simple analytical forms for which the G-functions are available 
(see, for example, the 05.08 series of the Royal Aeronautical Society’s Data Sheets 
on Stressed Skin Structures). 

Table XIV. Elements 4 and 5: Calculation of boom loads and shear stresses 

of the statically equivalent stress systems. 

Table XV. Elements 1 to 7: Boom loads and shear stresses of the statically 
equivalent stress systems. 

It is advisable to check the calculation of boom loads in Tables XIV and XV 


4 4 
with the criteria and Prr=0. 
r=1 r=1 
Table XVI. Elements 1, 2, 3, 4 and 6: The G, G and G functions. 
Table XVII. Elements 1 to 6: Calculation of the functions I1,. 
Table XVIII. Elements 1 to 6: Calculation of the functions w,.. 
Note that w,,.=0 since there is no torque or shear on this element. 


Table XIX. Elements 1 to 7: Calculation of the warping functions due to 
unit boom load functions. 


Table XX. Elements 1 to 7: The “three moment” equations and calculation 
of the self-equilibrating end load systems. 


Since the tip of element 1 is free P,.=0. Also, since element 7 has no shear or 
torque w,.=0 and II,.=0 (see 6.6.3). Hence the “ three moment ” equations (677) 
to (681) are, 


At joint (1, 2) 
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At joint (2, 3) 
At joint (3, 4) 


(Wo w re), w Wace + ll, = Wee + II;) (E3) 
From difference of warping between joints (4, 5) and (5, 6) 


W oP yc + aw Woe) Pu” W re = We + Wee +w, ll, (W so We) II, — W ecll, 
(E4) 
At joint (6, 7) 
W cot (Woo— Wee) Pro=Wooll, (BS) 
Having P.., Pye, P,. and P,, the P.’s and P;. are obtained from equations (667) 
and (673) which in the present case are, 
P,,.=P;,,— Il, P,.=P,,- 
P,.=P,,— Il, ; Ul, (E6) 
P,.=F,.~ 
Table XXI. Element 3: Calculation of total boom loads and shear stresses. 
Table XXII. Elements 1 to 7: Total boom loads and shear stresses. 


4 
Check calculations of boom loads with = P,=0. 


Fig. 46. Boom loads P, and P;. 

Fig. 47. Shear stress in rear spar web, wall (12). 

Fig. 48. Shear stress in frort spar web, wall (34) or (34). 

As a matter of interest Figs. 46 to 48 also show the effect of considering the 
wing built-in at the fuselage pick-up plane, joint (6, 7). 

When calculating the boom stresses from the boom loads the actual boom 
areas corrected for the effect of the walls should be used. 

The maximum divergence from the statically equivalent stress system occurs 
each side of the open element 5. The boom loads are changed by up to 10 per cent. 
in the front spar and 40 per cent. in the rear spar. The corresponding changes in 
the shear stresses are 50 per cent. and 150 per cent. Considerable axial constraint 
stresses occur when the wing is built-in at the joint (6, 7). The constraint stresses 
are small in the outer wing (elements 1, 2 and 3) . This is due to the low torsional 
stiffness of elements 1 and 2, the small torque about the zero warping axis, and to 


the fact that the transverse variation of the boom areas is practically identical in 
all three elements. 


Shear flow in ribs 


An approximate estimate of the shear flows qx in the ribs at the joints between 
elements, may be made from the average difference between the shear flows in the 
covers at each side of the joints. 

Thus 


GR ia, 441) =2,240 x 0.048 (Fa41,¢ — Ib./in. 
With the -’s from Table XXII the following table is derived. 


Rib (g,¢+1) 0,1 3,4 3.6 6, 7 
— -1.18 -2.30 340 -204 -495 5.13 -—3.12 
— -118 -2.30 3.39 -204 -497 5.19 -—3.39 
GR Ib./in. -247 365 -219 -533 555  -345 


Note that at joints (0, 1); (1, 2); (3, 4) the differences are equal, since S, is continuous 
Over these joints. 
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TABLE 1a. 
BASIC GEOMETRICAL PROPERTIES OF TUBE ELEMENTS 
ELEMENTS 
WITH wil, Ixy in’ B. int | Xin in 
NOSE CELL 
0-88 | 2-764 | 22-494 
| 2 0:80 |-6- 
1182-7 230-8 | -81-30 
3 1-81 |-6-779 |-9- 746 
4 2-00 | 7-44 |-9- 745 
1-40 | 3-625 | 29-321 
| | 2-86 |-9-035 -42-949 
| 4 3-20 | 9-785 |-42-949 
2-44 | 4-522 | 34-809 
4 | 8485-2 | 1506-4 | -556-23 | 
3 4-3 [-4-078 271 
| 4 4-92 | 42422 
4-994 37 - 746 
| 2 3-0 |-4-039 | 37-746 
6 13,500 | 2462-4 |-946-86 L— es 
3 5-9 |-42-019  |-49 - 024 
4 6-6 | 13-271 |-49- 024 
ELEMENTS | 
WiTHOUT JI, Tyy int B. x in| y, in 
NOSE CELL 
1-84 4-038 | 32-856 
| -9-352 | 32-856 
3 |5297-2 4029-5 -358-49 | 
| 3 3-62 |-10-472 |-14- 564. 
4 4-05 | 40-958 564 
| 5-484 | 36-615 
| 2 2:9 |-40-836 | 36-645 
12,472 | 2484-7 |- 874-75 | 
| 3 5-2 |-20- 155 
4 5-7 | 43-474 |-20- 155 
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TABLE le. 
BASIC GEOMETRICAL PROPERTIES OF TUBE ELEMENTS 
| | | — | — | 0-048 | — |, 
a4 | 138-3 | 426-2 | 138-3 | 57-45 | 44-49 | 0-048 | 0-048 
| 67-7 | 408-2 | — | 32-25 | — | | — 
12 97-8 | 97-8 | — | 4493 _ 0-064 — 
| eee] — | 
243-2 | | 243-2 | 76-18 | 18-82 | 0-048 | 0-056 
H 154-2 | 490-3 42-76 — | 0-048 
42 448-5 | 4148-5 44°7 0 - 080 
«234-3 (260-5 ee | — = 
430-0 360-8 | 99-9 | 23-2 | 0-48 | 0-000 | 
“| 234-3 289-2 — | 
| 176-3 | 476-3 | — 16-02 | — | 0-080 | — 
EE | 278 [309-5 56-90 |oos | — 
34 439-4 | 4353-4 | 439-4 402-36 | 25-29 | 0-048 | 0-080 
278-4 | — | | — | | — 
| | 2 | Fa r 
123-+4 | 123-4 | 13°39 | 0-080 4 |294-4|450 -4| 444-0 |0-54037 
| 494-3 | 216-0 | 47-53 | 0-048 2 |390-0/294-4| 95-9 [0-75440 
306-7 | 306-7 | 24-43 | 0-064 || 3 50-2 [069097 
4+ 239-8 47-99 | 0-048 4 |586-2/533-8| 52-4 
476-3 476-3 | 46-02 | 0-404 5 |639-0|/586-2) 52-8 |0-91737 
| — || cveworicar | 48-7| 4 
4 6439-4 | 439-4 | 25-29 | 0-080 7 |CYLINORICAL | 63-0} 4 
44 278-4 | 343-6 | 57-45 | 0-048 NOTE: ELEMENT N°-7 1S OPEN 
AND HAS THE BOOM AREAS OF 
ELEMENT N°-6 
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TABLE Il. 

EXTERNAL LOADING ON WING DUE TO SYMMETRICAL GUST 
ELEMENT N°] STATION 
,{AIRCRAFT , on in. on. in. ton. in. 

515 150 -4 0-9 0 40 0 75 
=-10:776| 467 498-4 | -3-25 0-4 - 140 2-4 -20-8 
£,=-0050} aug 46-4 | 0-2 - 323 9-6 | -34-2 
m, =-0-445 

374 + 294-4 | -8-0 0-3 -654 24-6 - 47-5 
Be 323 | 342-4 | | | - 38-4 -15-7 
=-0-050 — 
m.=-0-445| 275-4+ | 390-0 | -44°5 0-5 - 4729 59-9 16-4 

3 275:-4- | 475-6 0-44 - 4729 59-9 265-0 
4337 

236-3+ | - 44-3 0-5 - 

46-039 4 4 2237 78-9 163°5 

2363- | 544-4 | 0-48 | -2237 78-9 378-0 

7 216-9+ | 533-8 | -42-3 0-54 - 2474 88-5 338-0 
1468) 490.7 560-0 | -43-05 0:56 | -2846 102-5 292-0 
7, =-17- 644 
£.=m.=0 | 4164-5+ | 586-2 | -13-8 0-64 - 3155 447-8 246-0 

5 464-5- | 586-2 | -40-3 0-54 | -3155 481-0 
€.- 460 
7, 2-19-20 438-4 642-6 - 44-95 0-605 - 3449 432-9 443-0 
| 444-74 | 639-0 | -13-6 0-67 | -3786 149-7 | 345-0 
6 44-7 - -40-4 0-59 | -3786 149-7 582:0 

= 4-60 
7, =-19-20 87-35 CYLINDRICAL! -42-2 0-665 | -4058 465-0 542-5 
| 63+ -14-3 0:74 | - 4380 182-4 443-0 

7 63 - CYLINDRICAL; 0 0 - 4380 182-4 0 

NOTES 


1. The signs + and - in the second column refer to the outboard and inboard sides 
of a station. 
The shears and torque pass into the fuselage at station 63. Thus element 7 


merely carries the bending moment across the fuselage. 


3. €, and 7], are co-ordinates Ww.r.t QX,QY) of the points where the ‘torque 
reference axis passes through the root cross-section of each element. £, and m, 
are the direction cosines of the torque reference axis wrt the right - handed 
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QY, QZ. 
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TABLE IV. 
ELEMENT 6: CALCULATION OF CROSS-SECTION FUNCTIONS FOR THE 
ENGINEERS’ AND BREDT-BATHO SHEAR _ STRESSES 
4 
12 23 34 4A = DERIVATION 
| rr 
| 200-3 | 4485-4 | 2132-5 | 1496-9 FORMULA (343) — 
Values of s andt | | 
$44 | 200-3 | 4185-4 | 275-3 | 496-9 | 2857-9 | Equation (206a) | 
creer | 30-642 | -2-475 | -73-387/ 14-202 Equations (687),G87q) 
and Equations (684), |), / 
5-402 | 118-340 | 6-098 |-119-460 (684) using Table 
De ney | 30°642 | -2-475 |-9-474 | 14-202 
Equations 692) |: 
A 
5-402 | 48-340 | 0-890 |-449-460 Equation (692) 
With suffix x 
= 6-002 replaced by y 7 
Ih, 0-45557 |-0-02094 |-0-08047 | 0-42017 
Equations (683) [ 
0-32448 Values of t, I, , | 
fh, | 0-48385 | 0-00137 |-0-48313 and Iy appear | 
— | — | 
176-3. | 309-5 | 557-4 | 343-6 | 4386-8 | Equations (2070 a 
| U'red 
5402 |-766  |~40,906 | 4,880 |-34390 From present 
899 | 36,626 | 3524 (-44,046 0 table 
° 
given in 
00 resent table | 
|0- 0029808 0-0049846 | 0-0044779\0 0049846 and tablet with 
i calculated from 
At 0-0003040; Equation (225) 
NOTES ers 
4 In Equation (225), viz., A =A /s , Ais obtained from Table Ka 
Thus AY =(1702-4)*/1386-8 = 2089-8 in?. 
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TABLE V. 


ELEMENT 5: CALCULATION OF CROSS-SECTION FUNCTIONS FOR THE 
ENGINEERS’ SHEAR STRESSES AND THE TORSION-BENDING STRESS 
SYSTEM 

| | 
34 | 44 42 DERIVATION 
r=3 | 
448) 15-354 | 34-424 Equations (688) and 688 a) 
“404-807 |-249-694 106-484 | using Table 
035187 | 0-14644 0-13830 | Equations 683) 
— Values of t, I, and I, 
40504 | ~0°36698 -0-08186 appear in Table I 
439° ined 343° 6 | 176-3 From Table I 
-27,000 5276 | 5,540 484 
From present Table 
052, 486 ‘|48,720 | -140,258 
NOT From Equation 548) 
REQUIRED | 87879 From Table I 
Ste 008708, 0 | 0- 010575 ‘Equation (647) and above row 
From Table | 
¥ 717-8 | 454-7 454-7 2345 
iA | 
Petree 166 207 477,670 | 39,759 | 36,274 449,907 |From above two rows , 
|0- 0010998 -0-0040998 | 0-0006967 |-0-0006967 | Equations 617)and 648) 
NOTES 


1, =, is obtained from equation (548) by advancing cyclically by 2 


2. ee the first row of equation 617) ‘s cyclic changes :- 


3. From equation 646) by cyclic changes :~ 


4, From the am TOW - equation 617) by cyclic changes :- 


the suffices and the elements ¥ the determinant. 
5 |34,44,42| =A +A. 326,343 


Thus, 42 °° 344 


hy 


326,343%/419,907 = 1,044,500 in’. 


Sens 
| 
) | 
I | 
| 
) | 
687) 
4), | 
| 
| | 
| 
| 
| 
| 
| | 
| 
ind b) | 
| 
| 
| 
le 
8 
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TABLE VI. 
ELEMENT 1: CALCULATION OF CROSS-SECTION FUNCTIONS FOR THE 
ENGINEERS’ AND BREDT-BATHO SHEAR STRESSES 
— 
23 34 44 DERIVATION 
r=2 
N° 
[665-2 | | Formula 343) 
7 Values of and t appear 
- |2956 | - | ‘in Table I 
[865-2 [2374 67-9 Equation (206g) 
| 747983 [-17-253 | - 2-434 | | (688 and 688a) 
| 177683 | 0-043 | -19-449 | Using Table Il 
747983 3-447 | - 2-431 | 
17-683 | 0-009 | 743-449 | Equations (692) 
Dye. a - | 0-034 | a With suffix x replaced by y 
hh, [St], 4] - 0-30845 | - 0-24944. 
| 
| -4-24896 | Equations (683) 
- Values of t, I, and I 
FS 0-348) 0-00046 | - 0-34259 | appear as in Table I 
| 
| 0-00060) - 
975 195-3 408-2 Equation (2074) 
e | From present Table. 
,0v Jo | | 
Notes 
1. The values of &,, 8,,, §,, and 6,, are not required in the present table; 
see however Table x 
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3 AND 4: 


TABLE VII. 


CROSS-SECTION 


FUNCTIONS FOR THE 
ENGINEERS’ AND BREDT-BATHO SHEAR STRESSES 


| 42 23 34 DERIVATION 
in, | 0-25548 — _|-0-68792 | 0-17653 
—  |-0-42487 — 
fh, |-0-47582 | —  |-0-49194 |-0-49262 | FOLLOW SCHEME 
2 OF TABLE UI 
by (ELEMENT 4) 
— | 0-02057%) 0 
[0-49087 |0-04720 |-0-54597 |0-47043 | SIMILAR, BUT SIMPLER, 
oe SCHEME TO TABLE IV 
3 [h,/Sd,.., }0-00748 |0-24792 | 000793 |-0-22440 | SEE NOTES 
0007057 |0-044762 | 0008822 | 0-044762 | NOTE = 4 
[O°46757 |-0-04947  |-0-08767 | 042859 
| — |-03557 | — 
= FOLLOW SCHEME 
00700 | 0-19949 | 0-00429 |-0-20186 | 
| 
by |0-00524 | @LeMENT 6) 
0 0035534) 0-0059223 | 0 -0053188|0-0059223 
NOTES 


1 Element 3 is a closed single cell tube. Thus. as in Table IV, find D 


from-(687) and (87a) and the remaining D,’s and D,’s from (684) and (684q ) 
The h/St functions follow then from (683). 
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TABLE VIII. 
CO-ORDINATES OF SHEAR CENTRE E 
ELEMENT 2 3 4 5 6 i 
reed 
3-459 20-530 4-054 1-638 24-203 | 4-838 | 
| 
-36:800 |-27-473 |-43-767 |-24-439 -24-473 | -26-206 
NOTES 
are the co-ordinates of the shear centre E, of | 
the root cross-section of each element wrt. the axes QX, QY lt, 
through the intersection of the diagonals (13) , (24) [See eg. Fig ‘ 
These co-ordinates are calculated from equations 689) or 689¢) 
Thus, for example, for the double-cell section 6, equation 689a) ee 
is appropriate, and with 1,,I, and I, from Table I; and aaa 
the summations on the RHS from Table IZ, one obtains ka 
after dividing by I, and I, respectively, me 
0-38457€.. + %, * 
+ 0070149. = 0 
TABLE VIIa. 
RELATIVE CO-ORDINATES OF THE SHEAR CENTRE WITH RESPECT TO mer 
: REFERENCE POINT E AND THE ZERO-WARPING CENTRE ie 
NCTES 
ELEMENT { 2 3 4 5 6 
2-26 19 34 -0 29 0-17 22 60 | 
| -2602 48 2-27 -6§ 53 -5:27 -7-04 
3.85 0-44 0-04 0:22 
84 - 16-49 -{4-77 86 
4 and are obtained from Tables VII. 
and Veg Pe obtained from Tables VI, K, X and XI 


2 
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ELEMENT 6: 


TABLE [Xa. 


CALCULATION OF CROSS-SECTION 
AXIAL CONSTRAINT BOOM LOADS AND SHEAR STRESSES 


FUNCTIONS FOR THE 


| 4 
| need 12 23 34 44 2 DERIVATION 
|309°5 557-4 | 343-6 | 4386-8 | FROM TABLE IW 
[22-42 [430-90 30-40 | 432-47 FROM TABLE 
[19842 [440-40 | 587-80 | 475-77 | 4702-39 
[278-4 439-4 | 278-4 | 1472-5 | FROM TABLE 
[22-12 [162-00 448-40 | 197-37 Ne 
95560 |-285,340 | 58,910 |-288,10 |-418,980 | 
i pat [45594000 | 68,684,000 | 12,606,000) 69350,000 531 | FROM ABOVE TWO ROWS 
2928,200 | 706,200 |-4,323,300|~4,091,700 |-4780.600 | AND Dye ev 
478600 |-33767,000 | 359,200 | 34,417,000 | 4,487,800 | FROM TABLE IY 
0.034058 -0-017184 | 0-015276 |-0 017484 EQUATIONS (96) AND (497) 
| | 0 017246 |-0-015244 EQUATION OF 574 
0°035998|-0-045244 | 0-008780 |-0-015244 EQUATIONS (204 5.74 
9 008436 ALSO THIRD LINE OF 
Lees 0°44998}-0-31717 018292 _|-0-31717 EQUATION(340) 
040545 WITH t FROM TABLE I 
aie To obtain K calculate first (A™)* from (240), noting that there is only one term under 


the summation. From Table I AS‘ = 944 Din? and from Table W 83,+85, 
hence (A™)* = (914)*/2448-6 = ‘ATin* From Table IV § 


thus with Ae = 4386: Bin2, Ke 1368 8- -[1386 82 + 2857- 


2. From Equation 492), N 


3. From Equation 663), @,= 


Also from Equation 662),@=@ 05 //196,234,000 = - 0-0987843 
rei t2K @where 2K = 0-0019404 


4. Use Equation (241) in the form, 


r,r+4 


roared 


= 2448-6, 
= 2857-9, 


= -0-44043 in? 


Use 204.) for the walls 34° and 34+ m the form, 


° t 


34 


= 39][278- 0 + 197-37)- (476-3x148-4 + 
- 418-980 
2x 1702-39 


| 
| 

; 
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TABLE 


ELEMENT 6: CALCULATION OF CROSS-SECTION FUNCTIONS FOR THE 
AXIAL CONSTRAINT BOOM LOADS AND SHEAR _ STRESSES 


Me 


r { | 2 3 4 DERIVATION 


ry 
> 


B 


r 


3:0 6°6 


FROM TABLE I 


717-8 | 747-8 «454-7 454-7 2345 


2 
Berean]! 456 440 | 474,750 | 35,043 | 31,326 | 394,260 |FROM ABOVE TWO ROWS 


| 20600 |-2-0600 | 14-3049 |-4-3049 EQUATION G40) 


— 4. From equation 493), with the summations given above and in Table Ika, 


=> / = 394,260/196,234,000 = 0-0020092 
Hence A= 0044824 in.” 
From equation @4),using the effective shear modulus G is 0:8G, 
= 0-044824 [0-8 x 3-85 x 105/s0x40°| 0-024876 in 
2. From equation (660) 
= 2A*/p 2345/0-024876 x 14,008 = 6-7296 in. 
3. Calculation of the co-ordinates €,,77, of the zero warping centre :- 
From equations 693), 
with the values given in TableI and Table Ka, 
0. 38457(€,-€,,)+ = 745-780 
Hence, and 75715: 863 in 
Thus using the values of €.,and 7, given in Table Vill , 
€, = 2055in. and n,= ~42°069in.. 
The important quantities calculated under Table [Xa and this table are tabulated below. 


K ®, | @ A | A | WM |] & | 


“0-043 | 44,008 | 37-69 |-423-05 |-0-0087843}0-044824 |0:024876/6:7296 |-0-059115) 2°055 |-42:069 


| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
1 
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TABLE Xa. 


ELEMENT 1: CALCULATION OF CROSS-SECTION FUNCTIONS FOR THE 
AXIAL CONSTRAINT BOOM LOADS AND SHEAR STRESSES 


| | ‘ 
2 23 34 DERIVATION 
975 |1953 | 1062 FROM TABLE YI AND EQU 639q) 
| FROM TABLE I 
feeeo | | EQUATIONS (9) 
"33,976,000 1626, 000 134348, ,000 |69,920,000| FROM ABOVE TWO ROWS 
te 658,400  -800 2,953,200 294,100 | FROM TABLE XL 
| 019059 006983 -0:019059 | EQUATIONS (55i) 
[0019059 -0-046974 019059 EQUATIONS 644) AND G42) 
|ALSO THIRD LINE OF FOU 
NOTES 
4. From equation (538), K-0, so that for rA1, 
2. The relation N=-Ay, follows from equation (546 with the suffices increased 
by 4 and the sign “changed. 
3. In equation 642) (must be found from equations 652) and (G49q) which in the 
present case take the form (See also equation (662)), 
From equation(240) and Tables I and W, (A")* = 287-9°/1492-5 = 55-535 in.*. 
Equation 548) becomes, 
Aaa (Ares 404K 226 - 195:3x 143-2 = 62,659 in*. 
169 Hence, @, = 62,659/2x55:535 = 564-44 
Also, = 69,920,000 +(62,6592/55,535) = 140, 617,000 in.*, - 44850 in® 
Thus, @ = 564-44/44,850 = 0047575 


4 
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TABLE Xs. 


ELEMENT 1: CALCULATION OF CROSS-SECTION FUNCTIONS FOR THE 


AXIAL CONSTRAINT BOOM LOADS AND SHEAR STRESSES 


4 
r 2 3 4 DERIVATION 
0-88 | 0-80 | 4-81 2-00 
FROM TABLE I 
res | 22670 226-0 | 143-2 143-2 738 
58,044 | 63,845 | 44,330 40,253 | 143,469 FROM TWO ROWS 
0-96474 |-0-96174 | 0-60938 |-0-60938 EQUATION (654) 
NOTES 


4. From equation (550), 
/ = 443,469 / 140,617,000 = 0-0040203 
Hence = 0°034942 
From equation (44), r= 0-034942 [3-85 x 10°/40x10°] = 0-049848 in” 
2. From equation (660) 
A= = 738-4 /0-019818 x 1,850 - 3-4424 in. 


3. Calculation of co-ordinates €,, 77, of the zero warping centre:- 
From equations (693), 
with the values given in Table I and Table Xa, 
0-35225 € -€..)+ (1,- = 7 14-490 
€,- 0-444 
Hence 1-240in.and ~M1-627in. 
Thus using the values of €,, and 1. given in Table YI, 
4:399in., and -48-427in.. 
The important quantities calculated under Table Xa and this table are 
tabulated below. 


| 


| 
| 
| 
1,850 |-138-3 564-44 (0-047575 0-034942 0-019818 3-4424 044949 | 4-399 
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TABLE XI. 
ELEMENTS 2, 3 AND 4: CROSS-SECTION FUNCTIONS FOR THE AXIAL 
CONSTRAINT BOOM LOADS AND SHEAR STRESSES 
12 23 34 AA DERIVATION 
49459 0-87048 0 FOLLOW SCHEME 
OF TABLE X 
— |-036986 | — (ELEMENT 4) 
9.9 SIMILAR, BUT SIMPLER SCHEME 
| 0-22479 |-0-04460 |-0:33154 FOLLOW SCHEME 
— 0-49470 | — (ELEMENT 6) 
a r 2 3 4 DERIVATION 
2 4-35520 |-4-35520 |0-85935 |-0-85935 
3 | 452758 |-452758 |0-96776 |-0-96776 | AS IN ABOVE TABLE 
4 174755 |-4-74755 | 4-40740 |-4-40740 
LEWENT| 
2 | 42,631 |-459-35 |-0-036367|0-041832 |0-023216 |4-4291 |24-384 | 32-307 
AS 
3 [9366-5 |-145-66 |-0-015554 0-063094 | 0-035045|4°9907 1-458 |-30:259 | 
4 [1,293 |-12-86 |-0-009994|0-05516 |0-03064 |5-7099 | 1-679 |-38-943 | 
NOTES 


1. The element (2) has a nose cell and no wall between booms 2 and 3. 
Hence the procedure follows that used for element 1. Note that all the appropriate 
formulae are obtained from those in the text of the paper by cyclic changes of 
two in the ‘wee? but no change of sign, for example in the present case 


2. The element (3) is a 
the rows from A 


toh 


closed single cell section. Thus in Table IX@) one needs only 


Gee 43 4)and finally L from equation 340) 
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GENERAL THEORY OF CYLINDRICAL AND CONICAL TUBES 
TABLE XV. 
ELEMENTS 1 TO 7: BOOM LOADS AND SHEAR STRESSES OF THE 
STATICALLY EQUIVALENT STRESS SYSTEMS 
STATION 
| 
55 | 0-33 40470 | -4-47 -0-567 |-2-350 |-0-254 |-0-568 
| [-2-67 | 278 | 9-70 | | NO /_|-4-426 [-5-898 | -0-439 |-1-338 
re REQUIRED 
[-5:26 | 6-59 22-86 | -23-49 /__|-4:868 |-7-287 |-0-982 |-1-737 
| | 49 38-56 |-39-24 | 56-40 - 2-082 |-7-637 4-566 |-4-93 
| 4-20 | 38-55 |-39-24 | 56-08 | 1-562 /| 4-683 |-9-090 | 0-149 
| 323 |-15-43 | 16.47 | 56-78 |-57-72 | 1°857 0-822 |-9-454 | 0-468 
| 2754+/-24-04 | 22-44 | 77-48 |-78-64 | 42-34 | 2-046 0-569 |-40-031 | 0-479 
STATION | Pee Peg Pe Tene | Teast Teas” Teas 
| 
| - 21-08 22°49 77-42 «| - 78°53 442-08 4-046 4°488 |-2-400 4-649 
2363+ |-25-20 | 26-91 | 92-25 |-93-95 | Nor | 3-410 | 2-684 |-4-378 / _| 2-862 
236-3- |- 25-20 | 26-91 | 92-25 |-93-95 [REQUIRED | 4.273 | 5-236 |-0-869 | 5-387 
| ae-g+ [-26-80 | 28-66 | 98-47 |-100-03 |142-74 | 3-778 | 4-428 | -4-433 
wry STATION | Pe, Pea Pes Pea Pe Tove | | Teast 
| [ae-g- |-29-60 | 34-41 | 96-47 |-98-28 [493-74 | 3-4 | 3-436 | 1-710 | 3-297 
4 [490-7 |-32-13 | 34-44 | 404-76 | -106-77 |pequineo| 2°786 | 2-605 |-2-897 | 4-242 | 2-760 
164-5+ |-34-37 | 36-56 | 442-10 |-114-30 | 55-47 | 2-564 | 2-229 |-2-926 | 0-903 | 2-382 
STATION 1 Peo Peat Pra Peat Pry Pe + Ceggt Urgy 
| | -37-64 | 39-87 | 109-99 |-112-26 | 444-77 | 6-927 /| 0+436 
5 | 138-4 | -54-84 | 57-30 | 105-42 |-107-58 5-808 | | 4-554 | | 0-448 
|-59°32 | 74-95 | 103-26 |-405-89 | 67-88 | 4-985 0-257 0-457 
STATION | £2 Peg Pea Pe Tet Tare C523] “esa | Case teat 
| - |-38-29 | 40-93 [422-91 |-125-55 | 3-542 | 3-658 |-3-435 | 2-400 | 3-875 
| | 87-395 |-40-99 | 43-90 | 131-68 |-134-58 | 3°740 | 3-428 |-3-848 | 1-682 | 3-673 
| 63+ | 47-42 [-445-31 [495-79 | 3-878 | 3-295 |-4-558 | 4-264 | 3-574 
| No STATION Peo Pes Pea Pe 
0T063- |-44-2) |47°42 | 442-44 |-145-34 | 195-79 0 0 0 


NOTES. The decaviec calculations for elements 4 and 5 appear in Table XIV. 

The calcuiations for elements 4,2,3 and 6 follow the scheme used for Element 4, noting that the 
variation functions are obtained from Table XIII and the cross-wise variation functions 
from Tables VI, Vil and IV. No special calculation is required for Element 7, it being as in Element 6. station 63+. 
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For cylindrical tubes, 


GENERAL THEORY OF CYLINDRICAL AND CONICAL TUBES 


TABLE XVII. 


ELEMENTS 1 TO 6: CALCULATION OF THE FUNCTIONS I, 


TIP. TORQUE UNIFORMLY DISTRIBUTED TORQUE] CONCENTRATED TORQUE Piao 
14100 | Tago x Tago go x Tago x Prgo | ton. | ton. | ton. 
4 2 3 Gogo 2 3 Gogo 3 Gogo ton. , 5 5 6 
[4-493 | 44-93 [64-2 ]0-794] 50-98] — — |  |62-94] 9-40 |56-40 |56-08] 9-42 


3 7.7640 370 0-625 |-64:38 | 180 |0-544 | 92-52 405-47 |-34-45 [442-74 [433-73 |-22-44 


4 +5*7065 487 |0-970 |o-559|-29-07] — — 1443-32 |-25-30 }455-47 144-77 |-44-60 

5 -94-50 [459-38 1169-24 |-104-36 

6 [5915/8 43 |0-822 1668-29] -43 10-463]-19-91] — — |648-38]-38-33 1195-79 |495-79 |-38-33 


NOTES. 4.The values of A and @are obtained from Tables X , XI and IX. 
2 The Tggo'S are the root values of the assumed Tg," given in Tables XII and XII. 
3.The values of Gy, appear in Table XVI . 
4 The value of ts obtained from equation G66) ,ViZ., = 
For the open element 5 the value Pr» is obtained from Table XIV. 
5 The values Pegg and are obtained from Table XIV and XY. 


6 The values of Tl, are obtained from equations @67)and 673), viz.:- 


* Prao® Pe grec 667), Ty Pe.get,c 673) 


TABLE XVIII. 


ELEMENTS 1 TO 6: CALCULATION OF THE FUNCTIONS Wee 


TIP TORQUE |UD. TORQUE | CONC TORQUE} 
agent 7, T T | “3/3 
‘{o-5t037/0 26048 | 369-2 [8.7093] 10 |-0-678|-6-78 [0-509 | 32-68] — | — [25-90] 6-75 | 16-75 |-4-459 


0-75410 |0°56867 | 649-4 |-3-4408| -64 |-0-400) 25-60 | — — |25-60 | 44-56 |-49-44 |-1-686 


4 0640-82929) 986 -9 487 |-0-487|-237-47 |0-227 |-41-80 — — 1-281 


2 
3 |0-89097/0 79383 848-4 |-0:72640| 370 [0-265 |-27-30 | 0-064 | 44-52 |256-73} 1-865 
4 
6 


{ 1472-5 |-0-38462} 843 | 0-555 |-451-22)0-444 |-6-49 | — — 1-357 


NOTES. 4. The values of gge are given in and the values of AG are obtained from Table I. 
2. The values ‘of ©, are given in Table X , XI and x. 
3. The Tage’s are the tip values of the Tg,’S given in Tables XII and XII 
The Tag's may be read directly from Table XVI 
4 The values of appear in Table 
5 The values of Wge are obtained from equation (664), viz.:- 


Gy 

5 For the open element 5 the voles Of Wage is obtained from equation 676), which in the 


Present case may be written, 
Wse EA* ro 
For ond [see TableV, for (Tas see Table XI and for A*and 9. see TableI 
Thus : wse= —6.1454 x 10-* x 6392{(1850/pe)[1 — po-+Inpe] +735 (Inp-)?} = 0.05681 in. 
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J. HADJI-ARGYRIS AND P. C. DUNNE 


TABLE XIX. 
ELEMENTS 1 TO 7: CALCULATION OF THE WARPING FUNCTIONS 
DUE TO UNIT BOOM LOAD FUNCTIONS 

N° x10? *9° in /ton. in-/ton. 
4 2:9737 - 4-004 0-040 2-976 -0-448 
2 {-5969 - 4-042 0-456 4-616 - 0-249 
3 083391 027 0-234 0-856 -0 495 
4 0:72875 | -1-074 0-383 0-780 = 0-279 | 
6 0-64559 | -4-195 0-654 0:774 - 0-422 


NOTES. 


2.The values G., and 


1 The values of pi and Aq are given in Tables X , XI and KX. 
zo are given in Table XVI 


3.The unit warping functions Wee and Wigo are obtained from equation (64) Thus, 
4 


° 


From Table XVI , 


tc 


s 


c 


° 
G 


c 


5 Ww 6 = 
tgo 4° Ag Ag Gs 


Wac 


Note that this relation is not generally true but happens to be 


true for the present case because t*=t®=4. 


and hence from equation(664) ~Wac 


4 For the open elements 5 and 7 only the W,.s are required. Thus from equation 676), 


Cc 
Wec 


Referring to Table I and YZ for element 5 andto Tables I and IX for element 7 (noting that 
the booms and areas are the same as element 6) one obtains, 


— 1° 7405 x 1075 x 55-109 = 6060x410" °x 63 4-042 x /ton 


ELEMENTS 1 


TO: 7: 


THE 


TABLE XxX. 


“THREE MOMENT” 
CALCULATION OF THE SELF-EQUILIBRATING END LOAD SYSTEMS 


EQUATIONS AND 


COEFFICIENTS OF UNKNOWNS IN-THE “THREE MOMENT” EQUATIONS x 40° RHS 
OF “THREE MOMENT! 
2c 3c Pac 6c Pre EQUATIONS x 10° 
4 
4-592 0-249 0 0 0 ~18-60 
- 0-249 2-472 0-495 0 0 25-88 
0 0-195 1-636 -0-279 0 40-79 
0 0 - 0-279 2-494 -0-422 ~ 99-48 
0 0 0 0-422 1-783 29-55 
ton. 0 -42-25 -26:57 
ton. - 12-84 59-44 1-76 


NOTES 4. The coefficients are calculated from equations €{) to 5) using Table XIX for the W's and W's. 


2. The R.H.S.’s are calculated from equations €4 to(ES) using Table XVII for the w's end Table XVIL 
for the TT’s. 


3. The framed values are those obtained immediately from the solution of the equations in the top table. 
4. The Pgq’s and are calculated from equations 6). 
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2.The detailed calculations for elements 4,2,4,and 6 are made as for element 3. 
3. The shear stresses are the same for elements Sand 7, as those given in Table XY 
The boom loads of Table XY are added to the boom loads obtained from P.. and P,. 
by means of relations (232). The resulting multiplying factors are the same as L, B/A 

in Table XXI. Thus,for example, the end load correction to boom 4 of element 5 is - 
0-3064 P,. = 0-3061x 59-44 = 48-09 ton. 
4.The maximum discrepancy in the boom loads at the joints between elements is 0:05 ton. 


or about 005%. This is a satisfatory check on the accuracy of the calculations. 


J. HADJI-ARGYRIS AND P. C. DUNNE 
TABLE XXII. 
ELEMENTS 1 TO 7: TOTAL BOOM LOADS AND SHEAR STRESSES 
IN. FROM 3 4 2 
No ¢ AIRCRAFT ton. ton. ton. ton. ton /in® | ton./in® {ton /in? | ton. /in’ 
55 |-0-33 | 0-33 | 4:47 [-4-47 “4-461 |-3-422 | 0-296 | 
467 [-3:24 [3-35 [9-33 [-9-44 [-6-003 |-0 357 | -1-422_ 
f-6-98 | 7-34 [22-40 |-2e-73 | 7  |-1-964 [-7407 |-0 697 
2] 374+ [412-24 | 37-89 |-38-57 -2:304 |-7:916 |-1-368 | -2-453 
|-4-56 | 42-25 | 37-89 |-38:58 | 0-590 2-453 |-10-°840 | 0-449 
2 323 |- 14-47 | 45°54 | 57-398 |-58-32 | 1-303 299 | 0-468 
2| 275:-4+ |-17-53 | 18-96 | 79-39 |-80-82 | 0-824 |~ 4-540 |-42-100 | 0-479 
ELEMENT 
275:4- |-17-52 | 18°93 | 79-38 | -80-79 | 5-24 | 3-404 |-1-776 
236:3+ |-24-54 | 23-22 | 9457 |-96-27 | 2-446 | 3-337 |-4-753 3-515 
236-3- |-24-54 | 23-22 | 94-57°|-96-27 | 4-874 | 4-679 |-0-548 4-830 
4] 246-9+ |-23-74 | 25-57 100-43 |-404-99 | 3-325 | 4-849 |-1-676 4-998 
ELEMENT a? ~? 
246-9- |-23-75 | 25-56 [100-18 |-101-99 | 3-576 | 2-808 |-2-603 [1-616 | 2-959 
4 190-7 |-24-84 | 26-85 |409-38 |-4414:39 | 1-364 [3-644 |-3-754 [1-528 | 3-766 
1645+ |-49-56 | 24-75 | 424-48 |-4123-68 |-1-285 | 4-945 |[-5-230 | 1-674 | 5-104 
ELEMENT 
yor | STATION P, 34 
164°5-|-49-52 | 24-78 | 124-45 |-423-72 | 6-927 3-084 0-436 
5 1438-4 |-36-75 | 39-24 | 446°58 |-449-04 | 5-808 1-554 “| 0-48 
gf 444-7+ |-54-23 | 53-86 | 44-72 |-447-35 | 4-885 0-257 0-457 
ELEMENT ~0 ~o 
17-7- | 53-66 | 4472 |-417-35 | 1-460 | 5-126 |-3-623 | 1-254 | 5-343 | 
6 87-35 |-49-90 | 52-84 | 126-03 |-420-93 | 2-824 | 4-055 |-4-056 | 4-324 | 4-300 | 
63+ |-52-34 | 55-55 | 136-96 |-140-46 | 4-135 |3 445 |-4-498 | 1-368 | 3 394 
ELEMENT 
7 ,{0 T0 63-|-52:34 | 55-55 [436-96 [-140-46 0 0 0 
NOTES. 4.The detailed calculations for element 3 are given in Table XXL. 


| 


| 
| 
3 
| 
| 
| 
| 
| 
Ft 
|. 
ie 
Lie 
610 
dhe 


GENERAL THEORY OF CYLINDRICAL .AND CONICAL TUBES 


‘fq pur woog 


‘Op 
00¢ 002 00! 0 
402 
d 40% 

409 

waysks ssaujs quajpainba 

Wood 400) 
402i 
40%! 


611 


4 
ki 

at 


DUNNE 


C. 


HADJI-ARGYRIS AND P. 


‘qom szeds UT sso.14s 


‘Ly 


00S 00b 002 001 
T T T 
wayshs 
$$2J}S 


/u09) 


$S31}S IDIYS 


612 


|: 
| 
| 
| 
|= 
it 
> 


GENERAL THEORY OF CYFINDRICAL AND CONICAL TUBES 


00S 00v 


10 (pe) ‘qom Jeds UI 


‘Sp 


002 


0 


T 


v 


L 


$1 qn} UDays 


waysks yuajDainba 
JO IDIYS 


$$24}S 


(ZI) ‘qom szeds 


ur sso13s aeous 


eg 
| 
| 
\ 
a 
4 


J. HADJI-ARGYRIS AND P. C. DUNNE 


APPENDIX 
APPROXIMATE METHODS 


This Appendix contains, 


(a) A method to obtain the statically equivalent stress system in conical or 
cylindrical closed or open tubes in which the similarity laws of the thick- 
nesses are not followed. 


(b) An approximate method for the calculation of the axial constraint stresses 
in conical tubes not satisfying the similarity laws of the thicknesses. 


(c) An approximate method of obtaining the diffusion stresses in direct stress 
carrying walls at cut-outs. 


No detailed derivations are given and the presentation has been kept as 
elementary as possible. The method (a) is the same in principle to other methods 
in general usage but advantage is taken of the various “ equivalent ” tube procedures 
developed in the main paper. The main problem of (5) is to find the torque dis- 
tribution giving rise to the axial constraint stresses. In (c) the diffusion problem in 
direct stress carrying walls near cut-outs is solved with the assumption that the 
areas of the adjoining booms are increased so that their stresses remain as calculated 
under (a) and (5). — 

In the whole of the Appendix the symbols x, y, B, 8 etc. refer to quantities at 
the or cross-section and not at the root only as in Parts I to VI, unless otherwise 
stated. 


A.l. THE STATICALLY EQUIVALENT AND CONSISTENT STRESS SYSTEM FOR CONICAL 
TUBES 
This section deals with the derivation of a statically equivalent and consistent 
stress system in conical tubes not satisfying the 1,.* and f* laws of equations (40) 
and (38). 


A.l.1. The n-boom multi-cell tube. 
The load in the rth boom is given by the ordinary engineers’ theory as, 


where M, and M, are obtained from equations (361). 
The shear flows q are calculated in two parts g; and gy where 


‘ : .  (A2) 


The shear flow qx is in equilibrium with the applied shear stress resultants Q,, Q, 
with the condition of zero twist throughout the tube. The additional shear flow qs 
balances the applied torque. gy, and qgz may be regarded as quasi-engineers’ and 
Bredt-Batho shear stresses. Note that while the boom loads are consistent with 
equation (272) for the engineers’ theory direct stresses, the shear flow gy is not 
consistent with the engineers’ theory shear stress equations (276a) and (2765). 

To determine gq, at a given cross-section of a multi-cell tube of the type shown 
in Fig. 16 the following procedure may be adopted. 

(a) Find at the cross-section the gradient of the boom loads which are 
denoted by, 


(b) Find the length thickness ratios 5°;,,,,and 8.4: and the equivalent 
single wall ratio 6 ..1from (206a). Also find 6= 2 rid 


| 
= | 
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(c) Find between any pair of booms m, m+1 the value gem. m+1 from, 


m—1 


m+1 = — (6,mDPx,)/6 a (A4) 
r=m+t 
Find the remaining q,.’s from the recurrence formule, 

(d) Find q;° and q;' from the relations, 
er, r410°r, r+1 = rer+l Er, r+1 Or. (A6) 

(e) Check by applying the condition, 

r= 


To determine qzg at a given cross-section, 


(f) Find the areas (A") from equations (207a), (207), (210a) and 
the area A” from equation (225) in the Corrigenda to Parts I to V. 

(g) Find the torque AT , of the shear flows gz about any convenient axis OE 
through the apex O from, 


n n 
= 2 (A°r, r+1 er, r+1 + Er, r+1 A’: r+1Q Er, r+1l. (A8) 
1 1 


Find the Bredt-Batho twisting torque T, from, 
(Note that in tubes satisfying the ¢*, t,* laws AT; would be zero if OE were chosen 
as the flexural axis.) 
(h) Find 
D' sr, r+] ar, r+1 (A®/ A?) 
from equations (226) in the Corrigenda to Parts I to V. 
(i) Find the shear flows gg from 


Ter 
r,r+1 > PD)° r, r+2 | 
(All) 
er, r41= D' sr, | 


If the tube contains an opening (n, 1) find gy immediately from 
der, r+1 = DP: e e e (Al2) 
1 


and hence qj, and q,' from (A6). To obtain gg use the equations 


If one or more of the single walls in a multi-cell tube are direct stress carrying, 
then over these walls equations (Al) and (A3) are modified as: — 


fe= tx + M, ty . (Ala) 
I, 

df: fr 

Df. (a (A3a) 


To illustrate the necessary modifications in equations (A4), A5) and (A8) assume 
that walls (1, 2) and (3, 4) in the four-boom tube of Fig. 15 in 4.3.2 are direct stress 
carrying. Then equation (A4) for the wall (2, 3) of the equivalent tube becomes, 
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4 2 
3 1 
where, 
2 
i= 
w= (4 
(A14) 
ds ds, 


Note that the integrals are taken over the equivalent single cell tube. Equation (A5) 
for the direct stress carrying wall becomes 
s 


des = + Df.ds . (AS5a) 


Finally equation (A8) takes the form, 


2 4 
( Gupds 2A’ + ( - . (A8a) 


w 


A.1.2.. The n-boom open tube 


The tube considered is as in Fig. 36 without J. In this case the complete 
statically equivalent stress system is given by the engineers’ theory and torsion 
bending stress systems. 

The engineers’ theory boom loads and shear stresses are found from equations 
(Al), (A3) and (A12). The torque 7, is now equilibrated by the torsion-bending 
shear flows and is given by equation (A9) with 


n 


1 


The boom loads Py are obtained from equation (611) which in the present case 
becomes, 


r 


where | is given by equation (604a). The origin of the radius vector of the swept 
areas A,, is the shear centre and the point « on the periphery is determined by 
equation (6035). Note that all cross-sectional properties appertain to the cross- 
section at r. The shear flows gv are determined from: — 


rT 
qrr, r+1 >= DPrm e (A16) 
1 


Note the correspondence between (A1), (A15) and (A12), (A16). For the special case 
of a four-boom tube Pr and qr are given by, 


ane gm | 


é 
it 
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where =,, is defined in (548) and all areas and =,, are now the values at the root. 
Note that only in the case of the four-boom tube are the shear flows independent 
of the variation of the cross-section when the ¢,* law is not satisfied. 


A2. THE AXIAL CONSTRAINT STRESSES FOR CLOSED FOUR-BOOM CONICAL TUBES 


This section gives an approximate method for the calculation of the axial 
constraint stresses in conical four-boom tubes (with or without direct stress carrying 
walls) which do not satisfy the ¢,* and ¢* laws of equations (40) and (38). 


A.2.1. Four-boom tubes in which the walls carry no direct stress 


The axial constraint stresses in a four-boom tube satisfying the ¢* and 7,* laws 
are calculated from the torque about the zero warping axis (see 5.6.1). In the present 
case the zero-warping axis does not exist. However, it is still possible to calculate 
the torque which must be applied, in addition to the given external loading, in order 
to reduce the warping of the tube to zero anywhere. Defining this torque as —T, 
and remembering that the tube is conical or cylindrical, it can be seen that under 
the external load and — 7, the stresses in the tube will be correctly given by the 
quasi-engineers’ and Bredt-Batho theories of Al. It follows that the axial constraint 
stresses in the tube are those due solely to a torque T, applied separately, the Bredt- 
Batho shear stresses from 7, being ignored. The torque 7, is given by the formula, 


4 


T, is given by equation (AQ), gx,, etc. are obtained from equations (A4) and (A5), 
and ©), is given by (663) for single cell or multi-cell tubes without openings or by 
(663a) for multi-cell tubes with an opening (4, 1). The areas A,, etc. are given by 
(191) and (191a) for a tube without openings and by (547) for tubes with an opening 
(4,1). Note that the terms in the summation in (A18) are taken alternatively positive 
and negative, the first term in 6,, being positive; for a tube with no walls between 
two booms the corresponding term in the summation is zero. 

A rough estimate of the axial constraint boom loads and shear flows may be 
made by applying 7, to an idealised tube of which the boom areas and skin thick- 
nesses are of the same order as in the actual tube and follow suitable t,* and /* 
laws. It is generally best to make the two tubes most similar at the cross-section 
where the axial constraint stresses are most marked, e.g. at the root if built-in. 

A more accurate estimate of the axial constraint stresses may be made by 
applying the “three moment” equations method of 6.6.3. For this purpose the 
tube is divided into a number of conical or cylindrical four-boom elements, of 
different cross-sections or tapers, each of which separately satisfies a t.* and t* law. 
The cross-section of each idealised element is made identical with that of the actual 
element at the middle section. The r* and 7,* laws are then chosen to make the 
cross-sections at the ends most nearly equal to those of the actual element. This 
procedure may also be used to find the statically equivalent stress systems; for 
example the wing structure analysed in Section 6.8 was treated by this method. 

By whatever method the boom loads and shear flows are found the final stresses 
should be calculated with the actual boom areas and wall thicknesses. 


A.2.2._ Four-boom tube with direct stress carrying walls 


In the case of the multi-cell tube of Fig. 15 with direct stress carrying walls (1, 2) 
and (3, 4) the torque T, is given by, 
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A fair approximation to the axial constraint stresses may be found by adding 1/12 
of the total direct stress carrying skin +stringer area to each boom and then pro- 
ceeding as under A.2.1. However, this 1/12 rule is rather inaccurate in wings in 
which the moments of inertia of front and rear spar differ widely. In such cases 
the following method should be applied. 

Denoting the effective area of the boom B, by B, +B, one has, 


AB, : (2+ +; Siotsia (2+ . (A20) 
where 


B,+OB, | 
(A21) 


7, Ay, B,+OB, 
and similar equations for \B, etc. From these equations it is easy to calculate 
AB, etc. by iteration which is rapidly convergent. If one assumes initially that 
AB, = AB,=AB,=AB,=0 the final values of AB, etc. are obtained with sufficient 
accuracy after two or three operations. 

For the singly symmetrical trapezoidal tube of Fig. 29, 


4— — | | 

(A22) 


and \B,=AB, and \B,=AB, may be obtained immediately from (A20). With 
the effective boom areas thus calculated the tube may be analysed as a four-boom 
tube by the approximate method given in Section A.2.1, or directly by Part IV if 
t,* and ¢* laws are satisfied. The resulting direct stresses are then assumed to vary 
linearly in the walls between two adjoining booms. The shear stresses are calculated 
from the direct stress gradient with the condition that the average shear flow is as 
calculated for the substitute four-boom tube. This method of finding the effective 
boom areas may also be applied when calculating the torsion bending stresses in 
four-boom open tubes with direct stress carrying walls. 


A.3. THE DIFFUSION STRESSES AT CUT-OUTS 


The method described in Section A.2 considers only the axial constraint stresses 
arising from the restrained warping of the four booms. The stresses arising from 
shear lag effects in direct stress carrying walls were not considered. These stresses 
are generally of major importance only at cut-outs in the direct stress carrying walls 
and in this section an approximate method for dealing with the problem is given. 

Consider a flat panel between booms | and 2. By the methods of Sections Al 
and A2 one obtains a stress varying linearly between booms 1 and 2. Ata cut-out 
the direct stress at the unsupported edge must be zero so that it is necessary to 
superimpose on the combined stress systems of Al and A2 a stress distribution 
varying from —c,, at boom | to —c., at boom 2. This stress distribution is 
balanced by boom loads. 


(7), +O a9) — (Wt,/ 12) at boom 1 
and 
(Fp +729) (019 — (Wt, / 12) at boom 2 


where w is the width of the panel. 

The diffusion problem is to determine the additional boom loads at any cross- 
section and the shear flows between the booms and the panel edges. It will be 
assumed that the boom stresses remain fixed at the values calculated by the 
procedures of Sections Al and A2. Thus the additional boom loads are carried by 
increasing the boom areas; the method solves the design problem of finding the 
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proper reinforcement of the booms near cut-outs. This reinforcement is, of course, 
ignored in the basic stressing of sections Al and A2. 
"Tt is convenient to consider separately the symmetrical and anti-symmetrical 
diffusion cases. For the symmetrical case the self-equilibrating stress system at the 
free edge consists of a uniform stress —(¢,,+¢.,)/2 in the panel and equal 
boom loads 
Poo = (Fp + F29) Wt /4 . (A23) 
For the anti-symmetrical case the stress varies linearly from —(c,,—©.,)/2 at 
boom | to (7,,—@2,)/2 at boom 2 and there are equal and opposite boom loads. 
at booms | and 2 respectively. 
At distance z from the cut-out the boom loads due to the symmetrical and 
anti-symmetrical diffusion are defined by P,, and P,, and are given by (see for 
example reference 17), 


i odd i? 
Pe 24 3 L p—zipz/2 | 
Pas ieven Ll 
where 
. : . (A26) 
4-0 
BOOM 4 Pao 300M 4 


= 


1:0 20 0-5 1-0 
z/w z/w 
Fig. 49. 
Diffusion of self-equilibrating symmetrical and anti-symmetrical stress-systems applied at the 
free end of an infinite uniform panel with constant edge stresses. 


In Fig. 49 the ratios P.,/P,, and P,,/Ps. are plotted against z/w for various values 
of uw. Having the increments of boom loads the additional boom areas AB, and 
AB, are obtained from 

Where «, and «, are the boom stresses derived under Sections Al and A2. The 
shear stresses are greatest at the panel edges and are easily obtained from the 
gradients of the boom loads. 
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FURTHER CORRIGENDA TO PART V 
P. 772. The first term on left hand side of equation 301 should read, 


P. 893. 


The last two of equations (374) should read 


tts 


Pc “im 
Pc Pci 
Note that whatever and 
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S. Camm, .Esq., C.B.E., F.R.Ae.S. 
A. V, Cleaver, Esq., A.R.Ae.S. 


Dr. G. P. Douglas, O.B.E., M.C., F.R.Ae.S. 

G. R. Edwards, Esq., M.B.E., B.Sc., A.M.Struct.E., F.R.Ae.S. 

A. G. Elliott, Esq., C.B.E., MS.AE., F.R.AeS. 

W. S. Farren, Esq., C.B.E., F.R.S., F.R.Ae.S. 

Sir A. H. Roy Fedden, M.B.E., D.Sc., M.I.MechE., M.S.A.E., F.R.Ae.S. 
Sir Arthur Gouge, B.Sc., F.R.Ae.S. 

Major F. B. Halford, C.B.E., F.R.Ae.S 


F. 

B. 
Professor A. A. Hall, M.A., Se 
S. Scott Hall, Esq., A.C.G.I., M.Sc., D.I.C., F.R.AeS. 
J. W. F. Housego, Esq., Grad.R.Ae.S. 
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P. G. Masefield, Esq., M.A., F.R.Ae.S. 

W. Tye, Esq., B.Sc., F.R.Ae.S. 

Captain C. F. Uwins, A.F.C., O.B.E., F.R.Ae.S. 

Honorary Librarian—J. E. Hodgson, Esq., Hon.F.R.Ae.S. 

Honorary Treasurer—Captain C. F. Uwins, A.F.C., O.B.E., F.R.Ae.S. 
Solicitor—L. A. Wingfield, Esq., M.C., D.F.C., A.R.Ae.S. 
Secretary—Captain J. Laurence Pritchard, Hon.F.I.Ae.S., Hon.F.R.Ae.S. 


ANGLO-AMERICAN CONFERENCE 

The President, Sir John S. Buchanan, C.B.E., F.R.Ae.S., and the Secretary, 
Captain J. Laurence Pritchard, Hon.F.R.Ae.S., Hon.F.1I.Ae.S., sailed for New York 
on 12th May 1949 to attend the second Anglo-American Conference which is being 
held under the auspices of the Institute of the Aeronautical Sciences from the 23rd 
May until the 4th June 1949. 

Over 60 members of the Society will be attending the Conference. 


ANNUAL GENERAL MEETING 

The Annual General Meeting of the Society was held at 4 Hamilton Place, W.1, 
on 5th May 1949 at which the Annual Report and Balance Sheets of Aerial 
Science Ltd. and Aeronautical Trusts Ltd. for 1948 were approved. The complete 
Report and Balance Sheets will be published in the July Journal. 


FELLOWS OF THE SOCIETY 


The following have been elected Fellows of the Society, as announced at the 
Annual General Meeting on 5th May 1949:— 


Edgar Hilton Atkin. 

James Valentine Connolly. 
Percy George Crabbe. 

Ivan Howard Driggs. 
Arthur Ernest Hagg. 
Robert Ernest Hardingham. 
Arthur Lorne James. 

Jules Louis Albert Jarry. 
William Henry Lindsey. 
John Lloyd. 

Elizabeth Muriel Gregory McGill. 
Basil Gervasse Markham. 
John Knudsen Northrop. 
David Macleish Smith. 
Richard W. Walker. 


THE ROYAL AERONAUTICAL SOCIETY CHARTER SCHOLARSHIP 
The Council have decided that to mark the granting of the Charter of Incorporation 
to the Society a Royal Aeronautical Society Charter Scholarship shall be founded. 
The conditions governing the award of the Scholarship are given below: — 


1. The Scholarship shall be known as the Royal Aeronautical Society Charter 
Scholarship. 

2. The Scholarship shall be awarded annually, or at such other times as may 
be decided by the Council, to assist a student wishing to undertake advanced 
work or study in aeronautics. 

3. A Scholarship or Scholarships up to the value of £300 a year can be awarded 
and these Scholarships, in the first instance, will be awarded for a 
period of one year, although an extension of the Scholarship for a second year 
will be looked upon favourably by the Council if it is clear that the student 
and the work would benefit by such an arrangement. 
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4, Students wishing to apply for the Scholarship shall write to the Secretary of 
the Society not later than Ist June (30th June this year, 1949) stating the 
course they intend to follow, and giving full particulars of their age and the 
results of any examinations taken by them. The Scholarship will only be 
awarded to those who can satisfy the Council that they have the necessary 
qualifications to benefit, and whose proposals for work or further studies are 
acceptable to the Council. Normally, the work should lead to some higher 
degree or post-graduate diploma, although this is not essential. 


5. The Council will select from the applicants those whom it considers suitable, 
and will call such applicants for interview in the offices of the Society before 
such Board as may be decided by the Society. 


6. The decision of the Council on all matters concerned with the award of the 
Scholarship shall be final and binding. 


It is intended that the first award should be made this year. Applications for the 
Scholarship should be made to the Secretary and must be received not later than 
30th June 1949. 


SOCIETY OF BRITISH AIRCRAFT CONSTRUCTORS SCHOLARSHIPS 


The attention of members is drawn to the Society of British Aircraft Constructors 
Scholarships administered by the Society. 

The Scholarships are designed for the assistance of young men who are unable, 
for financial reasons, to obtain training in aeronautical engineering. All holders 
of S.B.A.C. Scholarships are expected to qualify for a technical grade in the Royal 
Aeronautical Society. 

Applicants should be between the ages of 16 and 18 on the Ist September in the 
year of the award. The closing date for scholarship applications for this year is 30th 
June 1949. Application forms may be obtained from the Secretary. 


MEMBER’S NEW APPOINTMENT 


P. G. Masefield, Esq., Fellow, has been appointed Deputy Chief Executive, 
British European Airways Corporation. ’ 


THE AERONAUTICAL QUARTERLY 


The Council have great pleasure in announcing that the first number of The 
Aeronautical Quarterly is now available. Copies may be had from the offices 
of the Society, 4 Hamilton Place, W.1, at 10/3 each post paid to non-members or 
7/9 each to members of the Society. 


The contents of the first number are: — 
Control Reversal Effects on Swept-back Wings .. Pe H. Templeton 
Calculation of Downwash Behind a Supersonic Wing... G. N. Ward 
Estimation of the Effects of a Parameter Change on the 


Roots of Stability Equations .. K. Mitchell 
Flutter of Systems with Many Freedoms... .. W. J. Duncan 
Note on Propeller-Turbine Reduction Methods... wa E. C. Pike 


Determination of the Drag of Jet-Propelled Aircraft in Flight . G. W. Trevelyan 


and D. R. Blundell 
Notes on the Linear Theory of Incompressible Flow Round 
Symmetrical Swept-back Wings at Zero Lift .. és F. Ursell 


GRADUATES’ AND STUDENTS’ SECTION 
Visits 

A visit to Croydon Airport will be made on Saturday, 18th June 1949. Members 
Wishing to take part in this visit should write to the Assistant Hon. Secretary, 


Mr. D. A. Thurgood, Sunnyside Bungalow, Portsmouth Road, Ripley, Surrey, not 
later than ist June 1949, 
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BRANCH NOTICES 

Members of the Society are reminded that if they join the local Branch they 
will not be called upon to pay the Branch Subscription as this sum will be paid 
direct to the Branch by the Society. Society members should get in touch with the 
local Branch Secretary in their area if they wish to take part in Branch activities, 


BROUGH BRANCH 

The Second Annual Flying Display organised by the Brough Branch will be 
held at Brough Aerodrome (kindly lent by Blackburn Aircraft Ltd.) on Saturday, 
18th June 1949. Proceeds will be in aid of the R.A.F. Benevolent Fund and the 
committee are arranging a full programme of events which will again include the 
£100 Cirrus Trophy Competition for Cirrus-engined Auster aircraft. It is hoped 
that Branch officials and others will make a note of the date and if possible respond 
to the invitations which will be forwarded to them in due course. 

The display will begin at 2 p.m. Gates open at 11 a.m. 

Arrival time for visiting aircraft from 11 a.m. 

Admission to public enclosure 2/6, cars and motor cycles 1/-. : 

Further information may be obtained from the Secretary: F. A. Wilkinson, 
A.F.R.Ae.S., Design Office, Blackburn Aircraft Ltd., Brough, E. Yorks. 


JOURNAL BINDING 
The new prices of binding of Journals will be as follows :— 
1948 Volume .. .. 14s. 6d. (including packing and postage) 
Previous Volumes .. ‘16s. 0d. ” 
Cases for 1948 Volume 6s. 0d. ” ” 
Journals should be sent direct to the Lewes Press, Friars Walk, Lewes, Sussex, 
and the remittance to the Secretary at the Offices of the Society. 
Requests for cases, with remittances, should be sent to the Secretary at the Offices 
of the Society. 


CHANGES OF ADDRESS 
To assist in keeping the records of members correct and up to date the Sccretary 
will be glad if all members will notify him as soon as possible of changes of address. 
When notifying changes please give the following particulars:— 
Name (in block letters). 
Grade of membership. 
New address (in block letters). 
Old address. 
Changes of address must be received before the 15th of the month in order to 
be effective for the Journal for the following month. 


NEW MEMBERS 
The following new members have been elected : — 
Associate Fellows 


Basil Twyman Cheverton, Peter John Cooper, Raymond Edwin Cooper (from 
Graduate), John Frederick Foss, John Owen Hounsell, Horace Aubrey Jarrett (from 
Associate), Thomas Hubert Gethin Lewis, Henry Edmund North (from Graduate), 
Geoffrey John Powter (from Associate), Joseph Rawcliffe, Thomas Dorey Sellwood, 
Norman William Wakelin. 


Associates 


John Holt Eaton-Shore, Alexander Sabine Courtenay Lumsden, James Alfred 
Scott Mills, Frank William Proudlove, Frederick George Beaconsfield Randall, 
James John Chalmers Reid (from Student), Denis Maule Roberts, Stanley Speight. 
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Graduates 

Michael Henry Beeching, Brian Robert Carter (from Student), Paul Emile Joseph 
Ghislain Colin, Jal Erach Daboo, Albert Georges Renaud Dechamps (from Student), 
George Frederick Donno, David Dunlop Ewart (from Student), Ernest Alfred Timby, 
Paul Anthony Tunbridge. 


Student 
Eric Cyril Maskell. 


ACKNOWLEDGMENT 


The Council acknowledge with grateful thanks back numbers of the JOURNAL 
from N. B. Garrard. 


ADDITIONS TO THE LIBRARY 


Pamphlets in italics with location reference following in brackets. Books marked 
* and ** may not be taken out on loan. 


A.a.357—Streamline Flow. H. F. P. Purday. Constable. 1949. 

A.a.358—Mechanical Vibrations. W. T. Thomson. Prentice Hall. 1948. 

A.a.359—Introduction to Engineering Problems. R. Q. Brown. Prentice Hall. 
1948. 

A.a.360—Supersonics. R. W. Wood. Brown University, Providence R.I. 1939. 

A.a.361--Evaluation of Residual Stress. K. Heindlhofer. McGraw Hill. 1948. 

A.a.362—Plastic Deformation, Principles and Theories. L. N. Kachanov and 
others. Mapleton House, New York. 1948. 

A.b.82—Hydro and Aerodynamics. Albert Betz. F.I.A.T. 1948. 

BB.b.142—Analysis and Design of Airplane Structures. E. F. Bruhn. Tri-State 
Offset: U.S.A. 1949. 

D.b.281—Airline Record 1948. Roy R. Roadcap. Chicago. 1948. 

E.b.135—-Combustion Engines. A. P. Fraas. McGraw Hill. 1948. 

EE.h.42—Rocket Propulsion Elements. G. P. Sutton. Chapman & Hall. 1949. 

EE.h.43—Control of Airplane Flight Path by Propulsive Jets. H. McKinley 
Conway. Southeastern Research Institute, U.S.A. 1948. 

G.e.C.21—Nickel in Non-ferrous Castings. Mond Nickel Co. 1949. 

H.b.93—Law of Aviation. R. W. Fixel. Michie Co., U.S.A. 1948. 

M.a.34—Principles and Methods of Telemetering. P. A. Bowden and G. M. 
Thynell. Reinhold Publishing Corpn., N.Y. 1948. 

N.a.97—Science and Engineering of Nuclear Power, Vol. II. A. O. Allen and 
others. Addison Wesley Press, N.Y. 1949. 


N.A.C.A. Technical Notes 


1457—Overbalancing in residual-liquidation computation. A. S. Niles. 

1760—N.A.C.A. and office of Naval Research metallurgical investigation of two 
large forged discs of S-590 alloy. J. W. Freeman and H. C. Cross. | 

1765—Office of Naval Research and N.A.C.A. metallurgical investigation of a 
large forged disc of S-816 alloy. H.C. Cross and J]. W. Freeman. 

1789—Investigation of aerodynamic and icing characteristics of recessed fuel-vent 
configurations. R. S. Ruggeri, U. von Glahn and V. G. Rollin. 

1790—Investigation of icing characteristics of typical light-airplane engine 
induction systems. W. D. Coles. 

1800—Laminar mixing of a compressible fluid. D, R. Chapman. 

1899—The method of characteristics for the determination of supersonic flow over 
bodies of revolution at small angles of attack. A. Ferri. 

1811—The longitudinal stability of elastic swept wings at supersonic speed. 
C. W. Frick and R. S. Chubb. 

1812—The application of airfoil studies to helicopter rotor design. F. B. Gustafson. 
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1813—A study of flow changes associated with airfoil section drag rise at super. 
critical speeds. G. E. Nitzburg and S. Crandall. 
1814—Effects of several design variables on turbine-wheel weight. V. L. LeValle 
and M. C. Huppert. 
1815—On compressibility corrections for subsonic flow over bodies of revolution, 
Reissner. 
1816—Triangular wings cambered and twisted to support specified distributions 
of lift at supersonic speeds. B.S. Baldwin, Jr. 
1818—Effect of automatic stabilisation on the lateral oscillatory stability of a 
hypothetical airplane at supersonic speeds. L. Sternfield. 
1819—The response of pressure measuring systems to oscillating pressures. 
I. Taback. 
1820—Sitrength analysis of stiffened thick beam webs. L. R. Levin and C. W. 
Sandlin, Jr. - 
1821—Ambient pressure determination at high altitudes by use of free-molecule 
theory. B. Weiner. 
1822—Elastic and plastic buckling of simply supported metalite type sandwich 
plates in compression. P. Seide and E. Z. Stowell. 
1823—The buckling of parallel simply supported tension and compression members 
connected by elastic deflectional springs. P. Seide and J. F. Eppler. 
1824—Linearised compressible-flow theory for sonic flight speeds. M. A. Heaslet, 
H. Lomax and J. R. Spreiter. 
1825—Compressive buckling of simply supported plates with longitudinal 
stiffeners. P. Seide and M. Stein. 
1826—Linear theory of boundary effects in open wind tunnels with finite jet length. 
S. Katzoff, C. S. Gardner, L. Diesendruck and B. J. Eisenstadt. 
1827—Matrix methods for calculating cantilever-beam deflections. S. U. Benscoter 
and M. L. Gossard. 
1828—-Effect of forebody warp on the hydrodynamic qualities of a hypothetical 
flying boat having a hull length-beam ratio of 15. A. W. Carter and I. Weinstein. 
1829—Data on the compressive strength of T5S-T6 aluminium-alloy flat panels 
with longitudinal extruded z-section stiffeners. W.A. Hickman and N. F. Dow. 
1830—Tension properties of aluminium alloys in the presence of stress-raisers. 
I—Effects of triaxial stress states on the fracturing characteristics of 24S-T 
aluminium alloy. A. W. Dana, E. L. Aul and G. Sachs. 
1833—Notes on the foundation of the theory of small displacements of orthotropic 
shells. F. B. Hildebrand, E. Reissner and G. B. Thomas. 
1837—Elevated-temperature compressive stress-strain data for 24S-T3 aluminium- 
alloy sheet and comparisons with extruded 75S-T6 aluminium alloy. W. M. 
Roberts and G. J. Heimerl. 
1847—Review of emulsified antifoams for aircraft lubricating oils. W. W. Woods 
and J. V. Robinson. 


N.A.C.A. Technical Reports 


825— Analysis of wind-tunnel stability and control tests in terms of flying qualities 
of full-scale airplanes. G. G. Kayten. 

826—A method for determining the camber and twist of a surface to support a 
given distribution of lift, with applications to the load over a swept-back wing. 
D. Cohen. 

829-—Summary of measurements in Langley full-scale tunnel of maximum lift 
coefficients and stalling characteristics of airplanes. H. H. Sweberg and R. C. 
Dingeldein. 

832—A systematic investigation of pressure distributions at high speeds over five 
representative N.A.C.A. low-drag and conventional airfoil sections. D. J. 
Graham, G. E. Nitzberg and R. N. Olsen. : 

834—Effect of compressibility at high subsonic velocities on the lifting force acting 

on an elliptic cylinder. C. Kaplan. 
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843—Jet-boundary ‘and plan-form corrections for partial-span models with © 
reflection plane, end plate, or no end plate, in a closed circular wind tunnel. 


J. C. Sivells and O. J. Deters. 
844—Tank tests to determine the effect on planing-tai hulls of carrying length, 


width, and plan-form taper of afterbody. J. R. Dawson, R. C. Walter and 


E. S. Hay. 
848—The Lagrangian multiplier method of finding upper and lower limits to 


critical stresses of clamped plates. B. Budiansky and Pai C. Hu. 
849—The effect of wall interference upon the aerodynamic characteristics of an 
airfoil spanning a closed-throat circular wind tunnel. W. G. Vincenti and 


D. J. Graham. 
850—A visual and photographic study of cylinder lubrication. M. C. Shaw and 


T. Nussdorfer. 
853—Cylinder-temperature correlation of a single-cylinder liquid-cooled engine. 


B. Pinkel, E. J. Manganiello and E. Bernardo. 

855—Relation between spark-ignition engine knock, detonation waves, and 
autoignition as shown by high-speed photography. C. D. Miller. 

857—Analysis of spark-ignition engine knock as seen in photographs taken at 
200,000 frames per second. C. D. Miller, H. L. Olsen, W. O. Logan, Jr., and 
G. E. Osterstrom. 

860—Analysis of cooling limitations and effect of engine-cooling improvements on 
level-flight cruising performance of four-engine heavy bomber. F. E. Marble, 


M. A. Miller and E. Barton Beil. 
861—Experimental and theoretical studies of surging in continuous-flow com- 


pressors. R. O. Bullock, W. W. Wilcox and J. J. Moses. 


N.A.C.A. Technical Memoranda 
1199—Compression shocks in two-dimensional gas flows. A. Busemann, 
1202—On wind tunnel tests and computations concerning the problem of shrouded 

propellers. W. Kruger. 


A.R.C. Reports and Memoranda 
2140—Flutter at supersonic speeds: derivative coefficients for a thin aerofoil at 


zero incidence. G. Temple and H. A. Jahn. 
2233—Tables for use in the determination of profile drag at high speeds by the 


Pitot traverse method. J. A. Beavan and A. R. Manwell. 
2253—Measurements of the degree of smoothness attained in a laminar-flow wing 
specimen (Short Bros.). R. B. Coles. 
2256—A pplication of thin aerofoil theory to controls having set-back hinge balance, 
with an analysis of wind-tunnel data on aerofotls of finite thickness. H.H. B.M. 


Thomas and M. Lofts. 
2264—Flight tests at high mach number on E28/39 W4041 (single-engined jet- 


propelled aircraft). A. W. Thom, F. Smith and J. Brotherton. 
2341—Body interference on a tractor propeller. The results of 24 ft. tunnel model 


tests. R. J. Monaghan, 


College of Aeronautics 
24_Experiments on an induction type high speed wind tunnel driven by low 
pressure steam. G. M. Lilley and D. W. Holder. 


C.S.I.R. Division of Aeronautics 
ACA-38—Izod, tensile and ey tests on some aircraft steels of se 
manufacture. A. R. Edwards and F. G. Lewis. 


National Research Council of Canada 
MA-206—Comparison of experimental and calculated lift and hinge moment 
parameters for full-span control surfaces. R. J. Templin. 
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MA-207—Tables of the relative density of humid air. G: E. Rickwood. 

MA-213—Wall corrections for a two-dimensional electrical tank. J. Sanders, 

MM-210—IJnterim report on measurement of landing loads on a ski plane. J. P, 
Uffen and A. D. Wood. 

MR-4—Induced precipitation-preliminary experiments on seeding a variety of 
clouds with dry ice. D. Fraser. 

MT-5—Some dimensional considerations.in fluid mechanics. J. G. Laberge, 


J. LAURENCE PRITCHARD, 
Secretary. 
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“Co-operation like 


this 


makes things happen fast!” 


I'd be lost without the 
services of my local 
B.O.A.C. Appointed 
Agent. He’s the chap 
-multiplied 3,500 


times throughout the 
world — who saves 
me hours of preliminary planning. 
Yes, my Agent makes things happen in a 
hurry. And he knows the answers to all 
my riddles about currency restrictions, 
baggage allowances, visas and regula- 
tions. I give him my itinerary and leave 
the details to him. He handles everything 
without fuss or contusion. All I do is pack 
a bag and leave on schedule. 
This same Speedbird efficiency — this 
same concern for your comfort and 


pleasure — applies wherever you may 
fly along 150,000 miles of B.O.A.C. 
routes to forty-two countries on five 
continents. The flight itself is swift and 
sure — comfortable, too. You enjoy com- 
plimentary meals, and there are no extras 
for the countless attentions and courtesies 
that make your journey so much more 
enjoyable. It’s all part 
ot B.O.A.C.’s 30-year- 
old tradition of 
Speedbird service and 


experience. 


GREAT BRITAIN + USA + BERMUDA + CANADA 
MIDDLE EAST - WEST AFRICA + EAST AFRICA 
SOUTH AFRICA + PAKISTAN + INDIA + CEYLON 
AUSTRALIA » NEW ZEALAND - FAR EAST - JAPAN 


B.0.A.C. TAKES GOOD CARE OF YOU 


BRITISH OVERSEAS AIRWAYS CORPORATION IN 


ASSOCIATION WITH Q.E.A., S.A.A., T.E.A.L. 
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100 HOURS 
AT COMBAT RATING 


* 


yy **A 100 hour endurance test 
A under combat conditions, by a 
ag Rolls-Royce Derwent V turbo- 
Vii jet recalls a similar test 


completed by a Merlin piston 
engine in 1945. The Derwent 
V is the first turbo-jet in the 
world to complete such‘a test.” 


Flight 
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